The airborne sound reduction factor in decibels 
of a partition is defined as 


where N, is the energy incident on one side of the 
partition and N, the energy emitted from the other 
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Summary 


Field values of airborne sound insulation of floating joist floors are considerably lower than 
laboratory values. To investigate this difference the indirect transmission past such a floor with 
an average airborne insulation of 43 db has been measured in the field by a method previously 
described. The transmission past a similarly situated floating concrete floor, insulation 49 db, has 
also been measured for comparison. In both cases the value of airborne insulation was found to 
be determined mainly by flanking transmission and it appears that the higher insulation of the 
concrete floor in the same structure is due to the junction of the floor with the walls attenuating 
the flanking energy. There is little indirect insulation with the joist floor since it connects with the 
walls at very few places. 


Sommaire 


Les mesures effectuées sur des planchers flottants & poutrelles donnent des valeurs nettement 
plus basses que celles obtenues en laboratoire. Dans le but d’essayer d’expliquer cette différence, 
on a mesuré la transmission indirecte obtenue avec un plancher dont ]’isolation sonore contre les 
sons aériens était de 43 db. La méthode employée pour ces mesures a été décrite ailleurs. Les résul- 
tats furent comparés avec ceux obtenus pour un plancher flottant en béton massif avec une iso- 
lation sonore de 49 db. On a constaté dans les deux cas que la valeur de l’isolation sonore contre 
les sons aériens était déterminée en grande partie par la transmission latérale. Les valeurs plus 
éleyées obtenues pour le plancher en béton massif dans une méme construction semblent étre 
dues au fait que la jonction du plancher aux parois atténue l’intensité de la transmission latérale. 
Le plancher & poutrelles n’a que peu d’isolation indirecte puisqu’il ne rejoint les parois qu’en cer- 
tains points. 


Zusammenfassung 


Messungen der Luftschalldammung durch schwimmende Fuf8béden auf Balkendecken er- 
gaben in ausgefiihrten Bauten wesentlich niedrigere Werte als bei Laboratoriumsversuchen. Um 
diesen Unterschied aufzukliaren, wurde an einer ausgefitihrten Decke, deren durchschnittliche 
Luftschalldimmung 43 db betrug, die indirekte Ubertragung durch Nebenschliisse untersucht. 
Die angewandte Methode wurde bereits friiher beschrieben. Zum Vergleich wurde die seitliche 
Ubertragung neben einem ahnlich gelegenen schwimmenden FuBboden auf einer Betondecke 
mit einer Dammung von 49 db untersucht. In beiden Fallen ergab sich, daB die Luftschall- 
daimmung hauptsichlich von den seitlichen Nebenschliissen abhangt. Die hGhere Dammung der 
Betondecke unter sonst gleichen Verhaltnissen ist der Verbindung zwischen Decke und Wand 
zuzuschreiben, die die Starke des seitlich vorbeigeleiteten Schalles herabsetzt. Eine Balkendecke, 
die nur an wenigen Stellen mit den Wanden in Verbindung steht, besitzt diese indirekte Schall- 
dammung nur in bescheidenem MaBe. 


1, Introduction R= 20 log p,/p.+10 log S/A 


No. 6 


(2) 


where p, and p, are the sound pressure levels in the 
source and receiving room respectively, S is the 


iti he total ab- 
R=10 log N,/N. 1) area of the partition under test and A t 
0g N/Nz (1) sorption in the receiving room. Under laboratory 


side. When measured in the laboratory, Risgivenby negligible. 


conditions the energy transmitted by indirect paths, 
i.e. other than through the partition, should be 
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When the sound insulation of a partition is meas- 
ured in the field, however, some energy will travel 
to the receiving room via the flanking walls. There 
has previously been some evidence [1] that the 
magnitude of this indirect transmission did not 
become appreciable until the average insulation 
of the partition reached about 50 db. In fact, if the 
laboratory and field measurements of similar con- 
structions are compared it is found that the diffe- 
rences are usually small (compare Aston [2], [3] 
and Parkin [4]). There are, however, considerable 
differences between the laboratory and field results 
in the case of joist floors with floating rafts, as shown 
in Table I. (For both sets of results the frequency 
range of 200 to 2000 c/s has been used as the labora- 
tory results were only available for this range. The 
correction factor is in both cases 10 log S/A to facili- 
tate comparison.) 


Table I 


Measured insulation 


[db] 


Type of floor 


Laboratory Field 


(1) Plain wood joist. Floor- 34 
boards nailed to joists. 
Plasterboard ceiling. 


(2) As (1) but with 50 kg/m? 42 42 
of sand on plasterboard 
between joists. 


(3) As (1) but raft of floorboards 48 | 4] 
secured to battens floating 
on glass silk quilt over the 
joists. 

(4) As (3) but with 50 kg/m? of 57 44 
sand on plasterboard between 
joists. 


33 


It is seen that the difference between the two sets 
of values is significant as soon as the floorboards 
are separated from the joists by a resilient quilt 
and is most marked when pugging is included to 
increase the weight of the floor. 

The discrepancies may be due in part to the dif- 
ference in size between the laboratory specimen 
floors.and the actual floors (about 4.5 m? in the lab- 
oratory compared with about 14 m? in the field). 
This may affect the total stiffness of the floor and 
hence the insulation, but such an effect is unlikely 
to account for the whole of such a large difference. 
In any case such discrepancies are not found for 
other types of partition. The only other reason for 
the discrepancies could be the effect of the indirect 
transmission on the field results. This paper de- 
scribes the measurement of indirect transmission 
past a pugged, floating joist floor and, for compari- 
son, similar measurements on a concrete floor. 
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2. Method 


The measurements were made in the first block} 
of the Building Research Station’s experimental] 
flats at Abbots Langley. This block was divided into} 
two sections, one had floating joist floors and the 
other floating concrete floors. The load-bearing walls} 
were of 28 cm cavity brick, i. e. the floors were car-} 
ried by the inner leaf of brick, 11 em thick. A layer } 
of bitumenised felt was included in all walls at floor} 
level. The joist floor measured was of the fourth} 
type given in Table I and the details of its construc- } 
tion are shown in Fig. 1. It should be noted that | 


j 


Fig. 1. Axionometric projection of part of joist 
floor; 
(1) Skirting-board attached to wall but free from 
floor, 
(2) Floating wooden raft, 
(3) Glass silk quilt, 
(4) Sand pugging, 
(5) Pugging boards, 
(6) Bitumenised felt membrane, 
(7) Open-web steel joists. 


open web steel, not wooden, joists were used in this 


floor, but it makes practically no difference to the 


results obtained and the same discrepancy between 
field and laboratory results occurs for this type as 
for wood joist floors. The floating concrete floor, on 
which measurements were also made, was in a room 
identical to that containing the joist floor and sym- 
metrically placed in relation to the block so that 
all external factors were the same in both cases. The 
main floor, 15 em thick, extended over the whole 
of a flat and in each room the floating screed, which 
was 4 cm thick, rested on a quilt of glass wool turned 
up at the edges to isolate it from the walls. 
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A detailed description of the method for meas- 
uring indirect transmission has been given pre- 
viously [5]. Briefly, in the first stage of the method 
warble tones from loud-speakers are emitted in the 
source room. The normal vibration amplitude of 
all receiving room surfaces below are measured with 
a vibration pick-up. One measuring point to each 
square metre of surface approximately is used. 
Structure-borne sound exciters on the source room 
walls and floor are then used to replace at discrete 
points the energy previously obtained over the whole 
surface from the airborne sound field. These exciters 
are adjusted in position and number, and in the 
energy fed to each, until the vibration amplitudes 
in the receiving room are as nearly as possible the 
same as with airborne excitation. For stage 2 of the 
measurements the exciters on the source room 
floor are removed and each of the remaining ones 
fed with the same energy as in stage 1. This cuts 
out the direct component of transmission together 
with any indirect component that might be present 
in the receiving room walls due to the floor vibra- 
tions; the surface amplitudes in the receiving room 
are again measured. In stage 3 the receiving room 
ceiling amplitudes are kept at the lower level ob- 
tained in stage 2 while those of the remaining sur- 
faces are raised, if necessary by additional exciters, 


_ to the levels obtaining in the first stage. 


The sound pressure levels measured at each stage 
in the receiving room are then used to calculate the 
amount of indirect transmission. In the normal 
measurement of airborne sound insulation the field 
value of the reduction factor, R’, calculated by eq. 
(2) depends both on the directly and indirectly trans- 
mitted sound. If the latter were not present R’ 
would be increased to a value R’+AR, due to 
the direct sound alone, which would correspond 
more closely to the figure obtained in the labo- 
ratory. It is this quantity AR that is calculated 
from the sound pressure levels measured in the 
present method. 

Measurements were made at the three octaves 
centred at 300, 600 and 1200 c/s giving a total fre- 
quency range from 200 to 1600 c/s. The top and 
bottom octaves of the international frequency band 
for sound insulation were not covered because the 
frequency response curve of the vibration pick-up 
had a large peak in the upper octave and because 
at low frequencies the exciters used would not give 
high enough wall amplitudes. 


3. Results 


The graphs of airborne sound insulation as meas- 
ured by the provisional international method are 
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given for both the joist and concrete floors in Fig. 2. 
The average values for the frequency range 100— 
3200 c/s are 43 and 49 db respectively. 
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Fig. 2. (a) Airborne sound insulation of joist floor; average 
(100—3200 ¢/s) 43 db. 
(b) Airborne sound insulation of concrete floor; 
average (100—3200 c/s) 49 db. 
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Fig. 3. Transverse vibration amplitudes of all surfaces of 
the receiving room beneath the floating joist floor 
in the frequency range 400—800 c/s; 
(1) 
(2) x ---x Structure-borne excitation, stage 1, 


(3) 0 o Structure-borne excitation, stage 2 
(0db = I1ms—%). 


Airborne excitation, 


Curve (1) of Fig. 3 shows the wall amplitudes meas- 
ured in the receiving room beneath the joist floor 
when the source room was excited with airborne 
sound at 600-+200 c/s. Reproduction of this pattern 
with structure-borne excitation in stage 1 gave 
curve (2) and removal of the exciters from the floor 
in stage 2 caused the levels to be those shown by 
curve (3). Apart from the ceiling, these levels were 
essentially the same as for stage 1 which made it 
unnecessary to go on to the third stage at this fre- 
quency. Table II gives the sound pressure levels 
appropriate to these curves measured in the re- 
ceiving room, and also the calculated values of AR 
and the amounts of sound energy transmitted in- 
directly. 
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Vig. 4a) and b). Transverse vibration amplitudes of all 
surfaces of the receiving room beneath the floating 
conerete floor in the frequency range 400—800 c/s; 


Airborne excitation, 


(2) x---x Structure-borne excitation, stage 1, 
(3) 0 o Structure-borne excitation, stage 2, 
(4) Composite curve = airborne; 


Curve (1) above, for all surfaces except ceiling 
which is = structure-borne stage 2 from curve 
(3) above, 

(5) x---x Structure-borne excitation, stage 3 
(Qidb==1 anrs—=2)- 


Fig. 4 shows the corresponding set of curves ob- 
tained at 600-+200 c/s for the concrete floor. Cur- 
ves (1), (2) and (3) in Fig. 4a are for: airborne, 
stage 1, and stage 2 excitation respectively as be- 
fore. The full line (4) of Fig. 4b is a composite curve 
obtained from the stage 2 level for the ceiling and 
the air-borne excitation for the other surfaces. 
Curve (5) is the best reproduction of these values 
by stage 3 of the method. The measured sound pres- 
sure levels and values of AR calculated from them 
are again given in Table II. This table is completed 
by the results of measurements on both floors at 
other frequencies, 300-+100, and 1200-400 c/s. 


Table II 


Sound pressure levels Sound 


Fre- {db above 2-10~—4 dyne/cm*] we 
quency 3a | a MER? ee mitted 
: | | indirectly 
Airborne | | ae 
c/s | 
[e/s] excitation 8¢28¢2 Stage2 Stage 3 [%] 


Joist floor 


300+100] 52.0 | 50.0! 46.0 | 47.0 30 
600+200} 50.5 | 50.5 | 48.0 | 48.0 55 
43.5 | 42.5 36.5 | 75 


Concrete floor 


300-+-100} 46.5 46.5 | 39.5 | 45.5 80 
600+200} 48.5 47.5 43.0 | 46.5 60 
1200-400) 41.5 | 41.0 38.0 | 40.0 70 


- 

Note: Sound pressure levels for different frequencies on the 
same floor are not comparative since different exci- 
tation levels were used. 
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4, Discussion 


The accuracy of the method (+-1.5 db in the align- 
ment of individual vibration levels) is such that 
only the order of magnitude of the above results 
can be considered significant. Nevertheless, the figu- 
res show that for the joist floor at mid-frequencies 
the direct and indirect sound are about equal while 
the former is the controlling influence at low fre- 
quencies and the latter is at high frequencies. In the 
case of the concrete floor the indirect is at all fre- 
quencies of greater importance than the direct sound. 
Although the percentage of indirect sound is even 
greater in this instance than with the joist floor the 
absolute value is much less as is obvious from the 
difference in airborne insulation of the two floors. 
This fact is shown clearly in Fig. 5 where the vibra- 
tion amplitudes of the two receiving rooms under 
identical conditions of airborne excitation above 
are compared. It is seen that the wall amplitudes 
beneath the concrete floor are nearly everywhere 
less than those beneath the joist floor. Thus the 
sound radiated from the walls in one room is less 
than that in the other but this conclusion does not 
necessarily hold also for the floor and ceiling for it 
is here that the rooms differ in construction. 
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Fig. 5. Comparison of vibration amplitudes of all surfaces 
of both receiving rooms with identical airborne ex- 
citation in the source rooms above, frequency range 
400—800 c/s; 


(1) 


Amplitudes in room under the floating 
joist floor, 


(2) x ---x Amplitudes in room under the floating 
concrete floor 
(db =I mis—4)- 


The figures of 55 per cent. and 60 per cent. of 
sound transmitted indirectly at 600 ¢/s past the 
joist and concrete floors respectively show that the 
overall insulation between the two rooms would be 
increased by less than 3 db if there were no direct 
sound component. Thus with joist floors carried on 
11 em brick walls there would appear to be a limi- 
ting insulation of about 46 db obtainable at 600 ¢/s 
between adjacent rooms and correspondingly a limit 
of about 53 db with concrete floors on similar walls. 
Obviously an increase in the insulation of the floor 
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itself, to cut out this direct sound, will affect the 
junction of the floor and walls, and may tend to 
reduce the indirect sound also, particularly that 
component travelling from the floor into the recei- 
ving room walls. It is unlikely that the limiting 
value of 46 db would be appreciably increased by 
this means. If the load bearing walls were solid, 
say 23 cm thick, instead of 28 cm cavity construc- 
tion they would not be so easily set into vibration, 
hence would transmit less indirect sound (although 
being thicker they would be better radiators). In 
most practical building construction however only 
two, and possibly only one of the walls would be of 
23 em brickwork, the remaining partition walls 
being thinner. Hence the limiting figures quoted are 
probably of the right order of magnitude for load 
bearing brick constructions. 

All the components of indirect transmission tra- 
velling from source room to receiving room must 
pass the junction of the walls with the floor so it is 
significant to see how this junction differs in the two 
cases. From Fig. | it can be seen that the joist floor 
has practically no connection with the two side walls 
and only connects with the other walls at discrete 
points where the joists enter them. Thus the vertical 
continuity of these walls is only broken at these 
points which constitute a small fraction of the total 
perimeter. (The thin layer of bitumenised felt in the 
walls does form a barrier, but its effect is thought 
to be small and is certainly the same for both cases.) 
With the concrete floor, however, the main 15 cm 
floor extending over the whole area of the flat makes 
a definite junction with all four flanking walls at 
floor level extending over the whole perimeter. 

Consider what happens when structure-borne 
sound energy travels down from the upper room 
by the flanking walls. There is practically no atte- 
nuation when it passes the joist floor because there 
is no discontinuity, hence it causes high vibration 
amplitudes in the receiving room walls below, with 
correspondingly high sound radiation from them. 
In the other case the connection with the concrete 
slab at floor level causes part of the energy to be 
reflected back into the source room walls and part 
to be transmitted past the junction. Some of this 
transmitted energy goes down into the receiving 
room walls and the rest spreads sideways into the 
main floor slab. Cremer [6] has shown that when 
energy in a wall reaches a cross junction with a simi- 
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lar wall, five-eighths of it is reflected and the other 
three-eighths divided equally between the three 
transmitted paths. This type of junction occurs for 
two of the walls, while with the other two the floors 
form a T junction. Assuming that the order of ma- 
gnitude of reflected and transmitted energies is the 
same as for the other type it seems that approxi- 
mately one-quarter of the energy reaching the junc- 
tion arrives at surfaces from which it can be radiated 
into the receiving room. This corresponds to an in- 
sulation of 6db which agrees very well with the 
measured difference in airborne insulation. 


5. Conclusions 


The junction which a concrete floor makes with 
the supporting walls reduces the transmission of 
indirect sound, although more energy is still trans- 
mitted indirectly than directly. No such junction 
exists in the case of joist floors and, although the 
percentage of sound transmitted indirectly is less 
than with concrete floors, the absolute value is con- 
siderably greater. The net result is that a limit is 
set to the insulation of joist floors in load bearing 
brick constructions of about 46 db. This value can 
not be increased unless an effective sound insula- 
ting membrane can be used in the flanking walls at 
floor level. The obvious requirement for structural 
rigidity makes it unlikely that such a membrane 
can be devised. 
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MESURE DES INTERVALLES HARMONIQUES DE LA GAMME 
A PARTIR DE LA SENSATION SUBJECTIVE DE CONSONANCE 


par MEHDI BARKECHLI* 
Université de Téhéran 


Sommaire 

En vue d’établir une relation objective entre la gamme musicale et les caractéristiques physi- 
que des sons, l’auteur a repris, sur les intervalles harmoniques, les travaux antérieurs qui n’avaient 
porté que sur la mesure des intervalles mélodiques. 

Les résultats de cette étude sont les suivants: 

La sensibilité de l’oreille est trés supérieure & ce qu’on aurait pu croire a priori; 

La quarte parait plus «précise» que la quinte et loctave; 

L’ensemble des intervalles harmoniques appartient a la gamme naturelle; 

Un intervalle pris dans les domaines des sons graves (au dessous de 400 c/s) apparait moins 
«précis» que s’il est pris dans la quinte plus aigué. 

Ces expériences, ot la personne de l’exécutant joue évidemment un certain réle, confirment 
dans l’ensemble les résultats de recherches antérieures relatifs  l’existence de deux gammes: l’une 
mélodique concordant avec celle de PYTHAGORE et l’autre harmonique correspondant & la gamme 
d’ARISTOXENE. 


Zusammenfassung 

Mit dem Ziel der Aufstellung einer objektiven Beziehung zwischen der musikalischen Tonleiter 
und den GesetzmaBigkeiten der Akustik hat der Verfasser die friiheren Arbeiten tiber die harmo- 
nischen Intervalle wieder aufgenommen, die bisher nur die Ausmessung von melodischen Inter- 
vallen geliefert hatten. 

Die Ergebnisse lassen sich in den folgenden Punkten zusammenfassen: 

Die Empfindlichkeit des Ohres ist weitaus gréBer, als man das von vornherein erwarten wiirde. 

Die Quarte erscheint ,,praziser‘‘ als die Quinte und die Oktave. 

Die Gesamtheit der harmonischen Intervalle weist eine Verwandtschaft auf zur natiirlichen 
Tonleiter. 

Ein Intervall im Bereich der tiefen Téne (unterhalb 400 Hz) erscheint weniger ,,prizise‘‘ als 
in héheren Lagen. 

Diese Versuche, bei denen die Versuchspersonen offenbar eine gewisse Rolle spielten, bekrifti- 
gen in ihrer Gesamtheit die fritiheren Versuchsergebnisse, die die Existenz von zwei Tonleitern 
vermuten lassen: einer melodischen, die mit der pythagoreischen tibereinstimmt, und einer har- 
monischen, die der von ARISTOXENUS entspricht. 


Summary 


In order to establish a relationship between the musical scale and the physical characteristics 
of sounds, the author resumes, in relation to the harmonic intervals, previous work which was 
only concerned with the measurement of melodic intervals. His conclusions are: 

The ear has an actual sensitivity much greater than one would deduce a priori. 

The fourth sounds more ‘“‘precise” than the fifth or octave. 

Ths assembly of harmonic intervals forms part of the natural scale. 

An interval chosen from among the lower tones (below 400 c/s) seems less “‘precise’’ than one 
taken from the fifth above. 

These experiments, in which the personal equation of the player has a certain part, confirm in 
general the results of earlier researches, viz., that there are really two scales; one the melodic scale 
of PYTHAGORAS, the other the harmonic one of ARISTOXENES. 


Expériences préliminaires 

. Position du probleme 

. Méthode de mesure et appareillage 
. Dépouillement des résultats 
Erreurs 
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1. Introduction 


Le probléme s’est souvent posé d’établir une re- 
lation objective entre les gammes musicales et les 
caractéristiques physiques des sons qui les consti- 
tuent. Il est bien connu, que la notion de hauteur 
subjective, ne coincide pas rigoureusement avec 
celle du logarithme de la fréquence, comme on I’ad- 
met en premiére approximation. Les travaux de 
Syow [1], Srnvens, Vorkmann et Newman [2] 
ont montré que la hauteur dépendait en particulier 
de lintensité, du timbre et de la fréquence elle- 
méme: un méme intervalle paraissant plus grand 
pour les fréquences aigués que pour les fréquences 
graves. D’autres auteurs ont cherché 4 aborder le 
probléme, d’un point de vue plus rapproché de la 
musique et de nombreuses études ont cherché a vé- 
rifier s'il y avait parmi les différents systémes de 
gammes proposés au cours des siécles, un syst¢me 
qui correspondit plus étroitement que les autres 
aux intervalles mélodiques et harmoniques «na- 
turellement» utilisés. — L’objet de la présente étude 
est de fournir les premiers résultats relatifs 4 une 
mesure directe et objective des rapports des fré- 
quences des notes de la gamme présentant les 
consonances principales subjectives. 


2. Travaux antérieurs 


x 


Hetmnoxrz [3] fut le premier 4 trouver une 
justification physiologique a l’existence des nombres 
entiers comme base des intervalles de la gamme. 

Pour donner une vérification expérimentale de 
ces idées, Cornu et Murcaprer [4] effectuérent 
des relevés & laide d’un dispositif mécanique d’en- 
registrement des vibrations sonores, sur un cylindre 
enfumé (phonautographe). Ils ont discerné deux 
types de gammes, l’une mélodique coincidant avec 
celle de Pyruacore et l’autre harmonique déter- 
minée par l’extinction de «battements d’ordre su- 
périeur» entre les deux sons composants. 

Cependant, les exécutants ne se basent pas forcé- 
ment sur l’extinction de ces battements d’ordre 
supérieur pour apprécier la consonance et, d’autre 
part, l’appareillage expérimental utilisé par Cornu 
et Mercapigr, bien qu’ils l’eftissent perfectionné 
dans la suite de leurs essais [5], ne leur permettait 
pas de mesurer avec la précision nécessaire, les fré- 
quences composantes des deux notes sonnant simul- 
tanément: aprés avoir réalisé l'accord, l’exécutant 
reprenait pour la mesure, les deux notes séparément, 
ce qui introduisait un intervalle mélodique. 

La technique ayant fait d’énormes progrés depuis 
le développement de |’électroacoustique, la question 
a été reprise dans ces derniéres années. 


_ Smart [6] et Gruenz [7], étudiant les intervalles 
mélodiques par des procédés électriques, sont arrivés 
a la méme conclusion que Cornu et Mercapier. 

Le probleme des intervalles harmoniques, techni- 
quement plus difficiles, restait entier. 

Van Essrorck et Monrorr [8] ont poussé ces 
études sur un orgue électrique expérimental et 
«orthoclavier» que le premier a réalisé et dans lequel 
Voctave était précisement divisée en 53 commas. 
Ils ont fait jouer des virtuoses sur cet instrument, 
selon les différentes gammes réalisables grace a la 
division en neuf commas de chaque ton (Tableau I) 
et fait voter ensuite un grand nombre d’auditeurs 
sur la valeur esthétique relative des différents sys- 
temes pratiqués mélodiquement et harmonique- 
ment. D’un dépouillement statistique des résultats 
il semble que l’on puisse déduire, d’aprés ces au- 
teurs, que la gamme de PyTHAGoRE serait la plus 
«naturelle», tant au point de vue mélodique qu har- 
monique. 

Indépendamment des résultats de van Espronck 
et Monrort, l’auteur a abordé le probleme de la 
gamme d’un point de vue différent [9]: dans le but 
d’établir des relations entre la gamme persane, base 
de la musique orientale [10], et les gammes euro- 
péennes, nous avons mesuré les intervalles de cette 
gamme en enregistrant des oscillogrammes des mélo- 


Tableau I 


Tableau comparatif des gammes. Intervalles successifs avec 


les unités : comma, savart et cents. 


lcomma = 23,458 ~ 23 cents, 
lsavart = 3,986~ 4 cents. 

do ré mi fa sol la sido 

Rapports 
des fré- 26 ob ott 2 26 2b ott 

Gamme quences 
tempérée commas 8,8 8,8 4,4 8,8 8,8 8,8 4,4 
savarts 50 50 25 50 50 50 25 
cents 200 200 100 200 200 200 100 
do ré mi fa_ sol la sido 
Rapports 9 9 256 9 9 9 256 
des fré- — — ——_ —- — = —- 
Gamme de quences Pee ea Bats 
PYTHAGORE | commas Di Dou teathse BO Se 
savarts BL 51) 23 “bl bY ST 623 
cents 204 204 90 204 204 204 90 
do ré mi fa_ sol la sido 
Rapports 9 10 16 9 10 9 16 
Gaon fee SS > is 8 9 8 16 
@’ARisto. | W°n°es 

XENE commas D7 vee NE Lap eee gee PY ee 
savarts 51 46 28 51 46 51 28 
cents 204 184 112 204 182 204 112 


dies chantées par un grand nombre de sujets. Le 
dépouillement statistique des fréquences relevées 
sur ces oscillogrammes, a montré: 

1. La concordance parfaite dans un tétracorde 
de la gamme iranienne et de la gamme de PytHa- 
GORE en ce qui concerne les intervalles mélodiques 
principaux (do, ré, mi, fa). 

2. Les intervalles mélodiques secondaires que cer- 
tains musicologues avaient qualifiés de tiers et de 
quarts de ton, représentent un systéme de deux 
subdivisions distinctes d’un ton dont le schéma est 
donné en Fig. 1, entrant toutes deux dans le sys- 
téme de PyrHagore: L.L.C et L.C.L1. 


do 1 2 Ble 
! 4 ! 
cabal assests re 
i} 1 
bavi) gett et 
1 i} 
per 2h a she ae 
1 ' 
Se EO arta a | 
do 7 Jeane 
1 1 
ee ey Oe See ee epee Lh a 
do 7 2 ré 
' 1 1 
bade St a eee Orme nak ete 


Fig. 1. Divers intervalles dans la gamme iranienne, con- 
tenus dans un ton majeur (9/8). 


Nous nous sommes proposés de reprendre la ques- 
tion du point de vue de Cornu et Murcapinr, 
mais en évitant les inconvénients signalés, grace a 
VYemploi d’un appareillage moderne, et en nous 
basant essentiellement sur la consonance subjective. 


3. Expériences préliminaires 


Nous avons fait diverses expériences pour con- 
denser les résultats, souvent assez discordants, don- 
nés par différents auteurs et avoir une idée de la 
sensibilité de Voreille dans lappréciation des ac- 
cords. La plupart de ces essais ont été effectués en 
chambre «sans écho» (pour éviter l’influence des 
ondes stationnaires) et en «écoute biauriculaire». 

Les expériences faites avec des oscillateurs ont 
montré qu’on devait distinguer dans l’appréciation 
des accords, entre une oreille acoustiquement sen- 
sible et une oreille musicale. 

Pour l’accord parfait (lunisson) — obligatoire- 
ment mélodique si on veut éviter les battements — 
on peut avec des sources sinusoidales discerner 


1 Par les lettres L et C nous entendons le limma et le 
comma pythagoricien, soit: 


ae 23 90 ts et pee 5 20 t 
~ 236 ou 90 cents e 524288 o ou 20 cents, 
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deux seuils différentiels: le premier correspondant — 
a la perception d’un désaccord, va de 0,5 pour 1000 
a 2 pour 1000 selon la finesse de Voreille du sujet; 
le second correspondant 4 la distinction entre le son 
plus aigu et le son plus grave, qui seule intéresse le 
musicien est déja plus élevé, 1 4 2 pour 1000. 

Pour les autres accords, les erreurs deviennent 
beaucoup plus grandes. L’octave est apprécié avec 
un écart de 3 4 6 pour 1000, et la quinte et la quarte 
de 5 & 50 pour 1000 selon |’éducation musicale de 
Vauditeur. Pour les accords plus complexes, les 
écarts tolérés par l’oreille croissent trés rapidement 
et, & partir de la tierce, un auditeur non musicien, 
méme entrainé, confond les accords voisins les uns 
des autres. 

Cependant un entrainement 4 l’audition des sons 
sinusoidaux, rarement rencontrés dans la pratique, 
permet de réduire appréciablement ces écarts. 

Les expériences effectuées avec les instruments 
de musique, trés riches en harmoniques, donnent 
des résultats trés différents, comme l’ont fait remar- 
quer de nombreux auteurs: c’est la Vorigine des 
discordances entre les résultats fournis par les an- 
ciens expérimentateurs qui opéraient soit avec les 
cordes, soit avec des diapasons, et ceux obtenus par 
les physiciens qui ont repris ces expériences dans 
les vingt derniéres années, comme SHowER, Bip- 
DULPHE [11], Srevens et Davis [12]. 

Avec les cordes tendues du type sonométre, nous 
avons repris quelques uns des essais effectués 4 
Vaide d’oscillateurs et nous avons obtenu des écarts 
beaucoup plus faibles. L’écart tolérable pour recon- 
naitre un type d’accord restant pourtant assez large, 
sauf pour des sujets exceptionnels. Ainsi loctave 
est appréciée avec une tolérance de | pour 1000, 
la quinte et la quarte 1,5 pour 1000 et la tierce avec 
2 pour 1000. 

Dans tout le cours de ces expériences on a opéré 
dans la gamme médium s’étendant de 400 4 2000 ¢/s, 
ou la sensibilité différentielle de Voreille est 4 peu 
prés constante et 4 des niveaux sonores de 45 4 55 db 
environ, pour lesquels la sensibilité auditive est 
maxima. 


4. Position du probléme 


De nombreuses raisons nous ont conduit 4 penser 
qu’une justification de la valeur objective de tel ou 
tel systéme de division de l’octave devait se trouver 
dans l'étude de la consonance subjective des inter- 
valles harmoniques «otes jouées ensemble». On 
peut citer, par exemple, le fait que de toute anti- 
quité les exécutants font l’emploi des diverses con- 
sonances pour accorder leur instrument, méme dans 
une musique monodique comme celle de Vorient, 
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ce qui fait apparaitre la consonance comme la base 
méme de la gamme musicale. 

Les expériences sur la différentiation de la gamme 
-mélodique et de la gamme harmonique, sont de 
deux types: 

1, Le sujet est un musicien qui joue des mélodies 
ou des accords qu’il rend satisfaisants pour lui sur 
- son instrument; un appareil de mesure objective dé- 
termine d’autre part les fréquences correspondantes. 

2. On répete un méme morceau en utilisant suc- 
cessivement les différentes gammes (Pythagore, 
Aristoxéne, tempérée) et on demande aux auditeurs 
d’indiquer leur préférence. 

On fait appel dans les deux cas, 4 deux facultés 
différentes inégalement réparties selon les individus. 
Dans le premier cas, c’est la capacité (d’un violo- 
niste par exemple) de réaliser successivement ou 
conjointement des intervalles donnés par des taton- 
nements dont on est maitre (variation de la position 
du doigt sur la corde). Dans le second cas, c’est la 
capacité de discerner en bloc des phénoménes trés 
subtils que l’on appelle la justesse d’un accord, 
faute de pouvoir en dire plus. Ces deux capacités 
sont trés inégalement réparties chez les individus, 
et, si la premiére est nécessaire 4 la pratique d’un 
instrument de musique 4 sons non fixes, la seconde 
n’est reliée nia la sensibilité acoustique proprement 
dite, ni aux capacités esthétiques du sujet. 

Etant donnée la délicatesse des phénoménes en 
jeu, ’appréciation des auditeurs ne peut étre prise 
comme base dans la mesure objective des inter- 
valles. C’est le travail de musiciens de talent et bien 
entrainés dans ce genre d’expériences, fournissant 
sur leur propre instrument les intervalles mé- 
lodiques et harmoniques soigneusement réalisés, 
dont l’expérimentateur mesure la valeur objective. 


5. Méthode de mesure et appareillage 


Le probléme de la mesure des rapports des fré- 
quences dans l’intervalle harmonique revient a sé- 
parer les deux sons composants joués simultané- 
ment. On ne peut done avoir recours 4 un micro- 
phone qui ne capterait que la résultante des deux 
sons, donnant un oscillogramme pratiquement im- 
possible 4 dépouiller vu la complexité considérable 
des sons musicaux usuels, et d’autre part, comme 
on peut en constater la raison (voir ci-dessus 3), de 
nombreux auteurs ont démontré l’impossibilité 
d’utiliser des sons purs sinusoidaux pour les études 
d’acoustique musicale, le timbre jouant un réle 
essentiel dans la reconnaissance des intervalles. 

Nous avons donc mis au point un dispositif élec- 
trostatique permettant de capter séparément sur 
un violon les vibrations de chacune des cordes ex- 
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citées simultanément par lVarchet. La Fig. 2 donne 
le schéma du principe du dispositif: une mince la- 
melle métallique formant l’électrode était placée 
sous chaque corde et collée sur la touche du violon. 


Hilt 


au Te T C=1000pF 
Moteur Caméra|  |Préampli-| |Préampli- 
ficateur| | ficateur 
Gain=3000 


Commutateur 
électronique a 
trois directions 


Oscillo- 
graphe 


A 
1000c/s 


Fig. 2. Schéma de l’appareillage utilisé pour la mesure des 
intervalles harmoniques produits sur un violon en 
jouant des doubles cordes. 


Les cordes métalliques ou filées étaient reliées 4 la 
masse et les électrodes portées & un potentiel positif 
de 100 V environ 4 l’aide d’une pile. On construisait 
ainsi une sorte de microphone électrostatique et les 
déplacements de la corde provoquaient des varia- 
tions périodiques de capacité, se traduisant par un 
écoulement de charge dq/q=dC/C a travers une ré- 
sistance R intercalée dans le circuit de la pile. Les 
tensions alternatives recueillies étaient transmises 
& l’aide de faibles capacités a la grille de deux pré- 
amplificateurs attaquant un commutateur élec- 
tronique a trois directions et un oscillographe. On 
obtenait ainsi sur l’écran de celui-ci les vibrations 
des deux cordes correspondant aux deux sons com- 
posant l’accord, et on injectait sur la troisiéme trace 
une force électromotrice 4 1000 c/s issue d’un dia- 
pason entretenu restant constant 4 10—® prés et 
permettant de mesurer, s’il en est besoin, la hau- 
teur absolue des sons. 

La capacité d’une corde isolée tendue a quelque 
distance d’une plaque métallique, étant: 

Kel 
imam (4x/d) 
un accroissement de la distance correspond 4 une 
variation de capacité de C a C + dC, telle que: 
Kel 1 4dx 
C4 d0= | =f 
In os 
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La Jongueur d’une électrode fixe ne pouvant guére 
dépasser 5 cm pour ne pas géner le jeu des doigts 
de l’exécutant sur la touche, le systéme était done 
peu sensible et requérait une grande amplification 
de l’ordre de 3000, la tension produite ne dépassait 
guere 0,5 mV. 

La capacité formée par le systéme corde-élec- 
trode variait comme l’inverse du logarithme de la 
distance entre les deux armatures. La différence de 
capacité entre une méme corde et les électrodes 
voisines juxtaposées n’était pas trés grande et l’on 
pouvait craindre une influence de chacune des cor- 
des sur |’électrode de l’autre (Fig. 3), se traduisant 
par une intermodulation de nature électrique de 
chacune des cordes. 


Cordes 
Electrode. Touche 


Fig. 3. Influence électrostatique de chacune des cordes sur 
Vélectrode de autre. 


Pour obvier & cet inconvénient, nous avons été 
conduit 4 modifier la répartition des lignes de force 
du champ électrostatique 4 laide d’une électrode 
auxiliaire reliée 4 la masse et placée entre les deux 
lamelles, comme l’indique la Fig. 4. 


Fig. 4. Elimination de intermodulation de nature électri- 
que des cordes. 


Pour les enregistrements, une caméra 4 déroule- 
ment continu, a vitesse voisine de | m/s, utilisant 
du film de 35mm, était placée devant l’écran de 
Poscillographe dont le balayage était supprimé 
(Fig. 2). 

Les Fig. 5 et 6 donnent les photographies de l’ap- 
pareillage utilisé. 


6. Dépouillement des résultats 


La vitesse du déroulement de la caméra d’en- 
registrement étant loin d’étre constante, surtout 
au démarrage, il était impossible d’admettre que 
léchelle des temps était respectée. On a done eu 
recours, pour la mesure du rapport des fréquences, 
& la comparaison directe des deux oscillogrammes 
dont les figures 7 et 8 donnent des exemples typi- 
ques pour l’unisson et la sixte mineure exécutés par 
l’auteur. 


Fig. 5. Piéce composée de minces lamelles de euivre servant 
délectrodes positives pour les microphones (corde- 
électrode) captant séparément les vibrations des ac- 
cords, avec les électrodes auxiliaires reliées & la masse. 


Fig. 6. Piéce ci-centre (Fig. 5) montée sur le violon de mar- 
que GALIANO servant aux expériences réalisées. 
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Le film 4 dépouiller était placé dans un viseur & 
réticule dont la platine était munie d’un dispositif 
de déplacement micrométrique. 

On prenait sur chacun des oscillogrammes un 
repére facilement reconnaissable, tel qu’une pointe 
A, A’ (Fig. 9). Deux pointes successives A, A,, A’, A’, 
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Erreurs 
Les erreurs résultant de ce mode de dépouille- 
ment proviennent: 
a) de Verreur de visée d’un point reconnaissable 
sur le premier oscillogramme, 
b) des erreurs dues au déplacement du micrométre, 


ei 
LES RRA SO RARE RA RARER BES 


aa 


Meta tetaneeheta stabs ts Uhalptemiingte ane 


Fig. 7. L’unisson la;—la, = 1,00012 = 0,208 cents (écart ~ 0,20). 
1) las, corde libre; 2) la,, 4¢me doigt sur la corde ré (accord réalisé par 
lauteur). En bas les vibrations du diapason entretenu 1000 c/s. L’écart 
est donné par rapport & l’intervalle correspondant dans la gamme har- 


monique. 


Fig. 8. Sixte mineure: ré,—si D zen bas ré,, corde libre- en haut, si Be premier 


doigt sur la corde la (accord réalisé par l’auteur). 


a lorigine, ou B, B,, B’, B’, a lextrémité, repré- 
sentent done la. durée d’une période T pour le 
premier phénoméne et T” pour le second. On comp- 


Extremité 


Origine 


Fig. 9. Mode de dépouillement des oscillogrammes, entre 
A et B, on compte NV vibrations, et, entre A’, et B/’ 
N’ vibrations. 


tait alors sur l’un des phénoménes, le premier par 

exemple, W=300 oscillations complétes (entre A 

et B,), le nombre d’oscillations du deuxiéme phéno- 

méne pendant ce temps était la somme: 

a) de N’ oscillations complétes (entre A’, et B’ dans 
la figure), 

b) de deux fractions d’oscillations 4 l’origine et a 
lextrémité, représentées par les segments PA’, 
et QB’ de valeurs: 

PA’, QB’ 
AVA’, ° BB, 
Le rapport de fréquences est donc: 
PA’, QB 
_ AYA, +N BrBe 


d N 


c) des variations de vitesse du déroulement du film 
pendant une période unique a Vorigine et 4 l’ex- 
trémité, 

d) de l’erreur de lecture du micrométre. 


Il est facile de voir que ce sont les erreurs c et d 
qui sont les plus importantes, mais elles ne portent 
que sur 1/V du diagramme. Ce sont donc des erreurs 
du second ordre qu’on peut estimer négligeables. 

Les expériences ont été effectuées en deux grou- 
pes espacés par quatre mois d’intervalle, pendant 
lequel nous avons perfectionné l’appareillage, en 
particulier le systéme microphonique, pour éviter 
absolument toute influence électrique des deux vi- 
brations lune sur l’autre. L’aspect de l’oscillo- 
gramme obtenu sur les deux canaux correspondant 
aux deux cordes quand l’une d’entre-elles seule- 
ment est attaquée par l’archet, et l’autre maintenue 
par le doigt pour l’empécher de vibrer par réso- 
nance, a permis de voir que la diaphonie tant électri- 
que que mécanique est supprimée. Si on lache la 
corde inutilisée, celle-ci vibre a sa fréquence propre 
par l’effet du couplage mécanique, comme si l’archet 
excitait simultanément les deux cordes; chacune 
de celles-ci vibre sur les notes qui lui sont imposées 
par sa propre longueur. I] n’est pas douteux que ce 
phénoméne joue un rdle essentiel dans le timbre 
des sons produits sur les instruments de musique. 
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Tableau II 


accords et ceux des gammes de PYTHAGORE et d’ARISTOXENE 


Intervalles | | de PYTHAGORE | ad ARISTOXENE | Mesuré cart evel arene 
pe San F ‘ | PYTHAGORE | ARISTOXENE 
I if 


Rapport | Valeursen | Rapport Valeurs en | Rapport | Valeurs en | 
des savarts des savarts des savarts | 
|\fréquences) (et en cents) |fréquences, (et en cents) fréquences (et en cents) | 
| a | 
Fp 0 : OF as ot hn a nee 7 0083 
Pe i . 1,00012 
unisson. 1 0 1 0 0,208 
ré — ré 2 301,030 2 301,030 300,965 
3 4 3 is 1,99970 
octave I 1200,000 1 1200,000 L199,771 
la — mi 3 176,091 3 176,091 176,236 
Ge suc ns 1,50050 
quinte 2 701,954 2 701,954 702,532 
j——1 124,939 124,939 124,913 : 
i cuan i A 1,33326 
quarte 3 498,045 3 498,045 497,942 
‘ 81 5 
fa — la 102,305 a 96,910 96,914 
3 3 64 4 1,25002 
tierce maj. 407,819 386,314 386,329 
$ i, | 
fa — la 32 73,786 ‘6 79,181 78,594 | 
gs oe (as 1,19911 | 
tierce min. 294,134 | 815,640 313,34 
16 — si ee Sov 24n0 jinnies 221,849 221,683 
Se 8 ae | 2 1,66603 
sixte maj. 16 905,866 | 3 884,399 883,698 | 
la — fa 198,725 | 204,120 204,141 
s. 4 ii 8 1,60009 | 
sixte min. 81 792,213 | 5 813,686 813,784 | 
| 


Tableau III 


dans des régions différentes du spectre sonore. 
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Tableau comparatif des intervalles harmoniques mesurés & partir de consonances subjective des 


Valeurs en 
savarts 
(et en cents) 


+0,052 +0,052 
+0,208 +0,208 
—0,065 —0,065 
—0,259 —0,259 
+0,145 + 0,145 
+0,578 +0,578 
—0,026 —0,026 
—0,103 —0,103 
—5,391 +0,004 

| —21,490 +0,015 
+4,809 —-0,586 

| +19,170 —2,336 
—5,561 —0,166 
—22,168 0,661 
+5,419 +0,024 

a Peva | | +0,098 


Tableau montrant le changement de « précision» des intervalles harmoniques quand ils sont joués 


— ee ee 


d’ARISTOXENE 


| YTHAGORE ARISTOXENE 


Valeurs en 


Intervalles | 


de PYTHAGORE 


savarts 


(et en cents) \fréquences 


Rapport | 


des 


savarts 


| (et en cents) 


Valeurs en 


Mesuré 


| des 


| fréquences 


Rapport | Valeurs en 


savarts 
(et en cents) 


Y | 
_ Ecart avec | Ecart avec 
P 


Valeurs en 
savarts 
(et en tents) 


A RR I RR RR 
| } J | 


' Rapport | 
des | 
\fréquences) 
ré — la 3 
evs . 2 
quinte . 2 
la — mi . 3 
Siynid = 
quinte 2 
ré — si 128 
3 3 Pig 
sixtemin. | 81 
| 
Ne deeds sl BAe 0 
3° tk = 
sixte min, 81 


176,091 
701,954 
176,091 
701,954 


198,725 
792,213 
198,725 
792,213 


3 
2 


176,091 
701,954 
176,091 
701,954 


204,120 
| 813,686 
204,120 
813,686 


1,50235 


| 1,50050 


1,59588 


1,60009 | 


{ 
| 
| 
| 


176,771 
704,665 


176,236 
702,532 


203,001 
803,225 | 
204,144 
813,784 


+0,680 +0,680 
+2,771 +2,771 
+0,145 +0,145 
+0,578 +0,578 
+4,276 —1,119 
+17,013 —4,460 
+5,429 +0,024 
+21,571 +0,098 
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Les expériences ont été poursuivies aprés accord 
trés soigné des cordes utilisées, cet accord étant 
vérifié avant chaque expérience. Les cordes sol, ré, 
la, étaient en boyau filées en argent d’excellente 
qualité et la corde mi en acier chromé. La compa- 
raison. des résultats obtenus pour les mémes inter- 
valles, dans les deux groupes d’expériences succes- 
sifs, montre la parfaite stabilité des accords réalisés, 
leurs écarts restant pour tous les intervalles du 


-méme ordre de grandeur. I] va sans dire que les 


doubles cordes ont été toujours jouées sans vibrato. 


7. Résultats 


Les résultats sont condensés dans les Tableaux II 
et III, dans lesquels nous avons donné la valeur des 
différents intervalles mesurés en savart et en cents, 
pour un groupe déterminé d’accords exécutés par 
Vauteur sur son violon de marque GaALtAno, allant 


Fig. 10. Octave rés — ré, = 1,99970 = 1199,771 cents (écart ~ — 0,26). 
1) rés, corde libre; 2) ré,, 36me doigt sur la corde la (accord realisé 
par Vauteur). L’écart est donné par rapport & l’intervalle correspon- 
dant dans la gamme harmonique. 


fmeN tae ht wa Oe wel ook hve Sh yet typed? 


2B OL is Be Sle ee la es 


We eee 


eeeth CM 4 ae aes t 
‘ ak ; : 
Ape NS Sy A, 


Fig. 11. Quarte: mi,—la,—= 1,33326 = 497,942 cents (6écart ~ — 0,1); 
1) mi,, doigt sur la corde ré; 2) la,, corde libre (accord réalise par l’auteur). 


Fig. 12. Tierce majeure: fag 


la, = 1,25002 = 386,329 cents (écart = + 0,02). 


1) fa,, 26me doigt sur la corde ré; 2) lag, corde libre (accord réalise 


par l’auteur). 


Cot ks he ee en ee ee 


Fig. 13. Tierce mineure: fa,—la,—1,1911 = 313,304 cents (écart ~ — 2,34). 
1) fa,, 26me doigt sur la corde ré; 2) la;, corde libre (accord réalisé 


par auteur), 


Fig. 14. Sixte majeure : ré,— si;— 1,66603 = 883,698 cents (écart —= — 0,66). 
1) rés, corde libre; 2) si;, premier doigt sur Ja corde Ja (accord réalisé 


par lauteur), 
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des plus simples aux plus complexes avec leurs 
écarts respectifs par rapport aux intervalles corres- 
pondant dans la gamme de PyTHacorE et dans 
la gamme d’AristoxENnE. Les Fig. 7, 10 4 14 mon- 
trent les photos des enregistrements correspon- 
dants. 

On peut déduire de ces tableaux les conclusions 
suivantes: 

a) La sensibilité de Voreille est trés supérieure & 
ce qu’on aurait pu croire a priori; elle atteint facile- 
ment 1/10000 soit 0,2 cents environ pour les ac- 
cords simples, ce qui est bien supérieur aux préci- 
sions que l’on peut atteindre avec des instruments 
& son purement sinusoidal comme des oscillateurs 
(voir 3). 

b) La quarte, avec l’écart de 0,1 de cents apparait 

“comme plus «précise» que la quinte (0,6 cents) et 

loctave (0,3 cents) ce qui explique qu'elle ait été 
considérée par les anciens comme le cadre mélodi- 
que essentiel et Vintervalle le plus consonant. 
c) L’ensemble des intervalles harmoniques appar- 
tient sans aucun doute a la gamme naturelle, ainsi 
que le montrent les écarts respectifs des intervalles 
mesurés avec ceux des gammes naturelle et pytha- 
goricienne (Tableau IT), les premiers étant toujours 
inférieurs aux seconds pour les intervalles de tierce 
et de sixte majeures et mineures, qui différent no- 
tablement dans ces deux systémes. 

d) Le méme intervalle pris dans la région des gra- 
ves <400 c/s apparait moins «précis» que quand il 
est pris dans une quinte plus aigué. Le Tableau III 
en donne deux exemples pour la quinte juste et la 
sixte mineure. I] en résulte que l’oreille a une sensi- 
bilité 4 la consonance qui n’est pas constante et qui 
demanderait a étre déterminée avec précision: on 
pourrait voir son rapport avec la sensibilité diffé- 
rentielle mesurée par certains auteurs. 


Conclusion 


Dans les conditions un peu spéciales ou nous nous 
sommes placés pour effectuer ces expériences, cette 
étude montre, d’une facon trés nette, que la gamme 
«cnaturellement» utilisée en harmonie et dans l’ap- 
préciation des consonances, est basée sur la gamme 
d’ARISTOXENE et non sur celle de PyTHAGorRE. 
Des expériences antérieures ayant montré que la 
mélodie se basait sur la gamme de PyTHagorg, il 
en résulte une confirmation des conclusions tirées 
par Cornu et Mercapirr, relatives a l’existence 
de deux gammes bien distinctes, l’une pour la mé- 
lodie et l'autre pour ’harmonie, 
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Nous poursuivons actuellement ces expériences 
en cherchant 4 établir ces conclusions sur une base 
plus générale, en faisant appel 4 un nombre suffi- 
sant d’exécutants et de formation musicale diffé- 
rente. 

Cependant nous soulignerons que l’auteur — plus 
familiarisé avec la musique orientale de base Pytha- 
goricienne —a «naturellement», comme nous l’avons 
constaté, réalisé les intervalles de la gamme har- 
monique et que l’assentiment qu’il donne 4 l|’au- 
dition des quatuors 4 cordes bien exécutés par des 
artistes occidentaux, est un indice qu’il n’y a pas 
de différence flagrante entre les accords qu’il réalise 
et ceux qu’il entend. 

On ne peut d’ailleurs manquer d’étre frappé des 
variations observées entre les tolérances des diffé- 
rents intervalles. L’oreille apparait comme un or- 
gane un peu mystérieux, dont les réactions sont 
encore bien mal connues 4 beaucoup d’égards. 

(Regu le 19 Février 1952.) 
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DIE DIFFUSION IN DER RAUMAKUSTIK 


Von W. FURRER und A. LAUBER 


Forschungs- und Versuchsanstalt der Generaldirektion PTT in Bern, Schweiz 


Zusammenfassung 


Es hat sich in den letzten Jahren immer deutlicher gezeigt, daB die klassischen Groen nicht 
geniigen, um das akustische Verhalten eines Raumes hinreichend zu beschreiben. AuBer dem 
Volumen, der Nachhallzeit und der nach geometrischen Uberlegungen bestimmten Raumform 
spielt die Homogenitat des Schallfeldes, die man auch ,,Diffusion‘‘ nennen kann, eine sehr wich- 
tige Rolle. Im vorliegenden Aufsatz wird zunachst gezeigt, da die im stationaren Zustand vor- 
handenen Frequenzgang-Schwankungen, die schon langst bekannt sind, nicht direkt ein MaB fiir 
diese Diffusion darstellen. Dagegen gelingt der empirische Nachweis, dab die ,,mittlere Hohe eines 
Berges dieser Schwankungen ein wichtiger Parameter ist, der die Homogenitaét des Schallfeldes 
unmittelbar beschreibt. Er wird fiir 11 Raume, deren Volumina zwischen 37 m® und 20000 m?® 
liegen, bestimmt, und zwar jeweils fiir die Frequenzgebiete um 375 Hz und 1650 Hz. 


Summary 


It has become more and more evident in the last few years that the classical quantities are not 
sufficient to describe adequately the acoustic behaviour of a room. In addition to the volume, the 
reverberation time, and the shape of the room determined by geometrical considerations, the 
homogeneity of the sound field, which may also be called ‘‘diffusion”’, plays a very important part. 
It is shown in this paper that the frequency irregularities existing in the stationary state, which 
have been known for a long time, are not a direct measure for this diffusion. It is possible, how- 
ever, to prove empirically that the ‘‘mean height of a peak”’ of these irregularities is an important 
parameter directly related to the homogeneity of the sound field. It is determined — for the fre- 
quency regions around 375 c/s and 1650 c/s — for 11 rooms, whose volumes range from 37 m* 
to 20000 m%, 


Sommaire 


Depuis ces derniéres années on se rend de mieux en mieux compte du fait que les grandeurs 
classiques ne suffisent pas pour définir entiérement les qualités acoustiques d’un local. En plus du 
volume, de la durée de réverbération et de la forme du local déterminée d’aprés des considérations 
géométriques, lhomogénéité du champ acoustique, qu’on peut encore appeler < diffusion», joue 
un réle trés important. Le présent exposé montre tout d’abord que les fluctuations de la courbe 
de fréquence qu’on constate 4 l’état stationnaire (d’ailleurs connues depuis longtemps), ne repré- 
sentent pas une mesure pour cette diffusion. En revanche, on peut prouver empiriquement que la 
«hauteur moyenne d’un sommet» de ces fluctuations est un important paramétre qui définit im- 
médiatement ’homogénéité du champ acoustique. Ce paramétre a été déterminé pour 11 locaux 
@un volume compris entre 37 m% et 20000 m*, chaque fois pour les gammes de fréquence voisines 
de 375 et 1650 Hz. 


n einer friuheren Arbeit [1] wurde der 1940 er- 
| reichte Stand der akustischen Behandlung von 
Radio-Studios dargestellt. Er war im wesentlichen 
dadurch charakterisiert, da8 fiir jede Art und GréBe 
eines Studios die optimale Nachhallzeitkurve fest- 
gelegt werden konnte und daf es gelungen war, die 
dafiir notwendige Absorption in einfacher und 
sicherer Weise zu erreichen. Die Raumform wurde 
nur insofern als wichtig betrachtet, als es sich darum 
handelte, grobe Effekte wie z. B. Flatterechos, 
durchhiingende Nachhallgeraden usw. festzustellen 
und zu beheben. Von diesen Betrachtungen abge- 
sehen, wurde aber der Raumform wenig systema- 
tische Beachtung geschenkt; es wurde lediglich 


darnach getrachtet, reflektierende parallele Flichen 
zu vermeiden. 

1951 wurde mit der Fertigstellung des Neubaues 
des Radio-Studios in Bern, sowie Erweiterungs- 
bauten in Basel und Ziirich eine neue Bauetappe 
abgeschlossen, die in raumakustischer Beziehung 
wiederum durch erhebliche Fortschritte charak- 
terisiert ist. Dabei steht aber im Gegensatz zu friiher 
eine systematische Untersuchung des Einflusses der 
Raumform und der diffusen Reflexion im Vorder- 
grund. Im Laufe eines zehnjahrigen Betriebes 1940 
bis 1950 schalten sich namlich deutlich einige Bei- 
spiele von Radio-Studios heraus, die trotz einer 
ganz einwandfreien Nachhall-Frequenzkurve in 


bo 
ou 
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ihrem akustischen Verhalten nicht befriedigten und 
fiir die Regisseure und Studiotechniker Quellen stan- 
diger Probleme bildeten. Das charakteristischste Bei- 
spiel war das Orchesterstudio in Genf (3500 m), das 
bei seiner Eréffnung im Jahre 1940 als hervorragend 
beurteilt wurde, dann aber mehr und mehr Anla}B 
zu Kritik bot. Mit den stetig steigenden Anspriichen 
wurden immer neue Orchester- und Mikrophon- 
aufstellungen versucht, deren Erfolg jedoch unbe- 
friedigend blieb. SchlieBlich wurde 1947 ein Versuch 
mit den in der Zwischenzeit besonders in Amerika 
bekannt gewordenen polyzylindrischen Diffusoren 
gemacht, mit denen die grofen glatten Seitenwande 
beim Orchesterpodium verkleidet wurden. Dieses 
rein empirische Vorgehen war von einem durch- 
schlagenden Erfolg gekroént; das Studio wird seither 
als véllig in Ordnung befunden. Ein anderes ahn- 
liches Beispiel bildete ein Studio in Zurich (480 m?). 

Diese Erfahrungen zeigten klar die Notwendig- 
keit, neben dem Volumen und der Nachhallzeit 
bzw. der Absorption noch einen weiteren Parameter 
einzufiihren. Dieser soll den Homogenitatszustand 
des riumlichen Schallfeldes kennzeichnen, der fre- 
quenzabhangig durch die Art der Reflexion an den 
Begrenzungsflachen bestimmt ist. Es wurden zu- 
nachst ausgedehnte Versuche durchgefthrt, um 
das Interferenzfeld des Raumes im stationaéren Zu- 
stand auszumessen. Solche Messungen sind schon 
aus friitheren Arbeiten bekannt; Fret [2] hat 1935 
bei konstanter Frequenz die Schalldruckverteilung 
in Raumen gemessen, WENTE [3] variierte, eben- 
falls 1935, die von einem Lautsprecher abgestrahlte 
Frequenz langsam bei feststehender MeSstrecke, 
wobei der Frequenzgang des Raumes erhalten wird. 
1950 wird dieses Problem von Boutr und Roop [4] 
theoretisch behandelt ; es gelingt dabei, fiir den paral- 
lelepipedischen Raum mit gleichmaBig verteilter 
Absorption die Frequenzgiinge zu berechnen. 195] 
bentitzten SomERVILLE und Warp [5] fiir die Unter- 
suchung eines Modellraumes (Flachraum) als Kri- 
terium die abgewickelte Linge des Frequenzganges. 

Abb. 1 zeigt einen typischen Ausschnitt eines 
solchen gemessenen Frequenzganges. Als Frequenz- 
gang-Schwankung bezeichnen Botr und Roop: 


v 
= Pmin 


= Pmar ax 
Fy A aid “Ay 
Dabei bedeuten: 
=Pmax Summe der Schalldruckmaxima (,,Berge‘), 
=Pmin Summe der Schalldruckminima (,,Taler“), 


Ay  Frequenzbereich. 


Es schien zunichst, daB diese GréBe Fy ein wich- 
tiges Kriterium fiir das akustische Verhalten eines 
Raumes, ganz besonders eines Radio-Studios, sei, 
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Da ferner, wie bereits erwihnt, eine zuverlassige 
Statistik iber die subjektive Beurteilung einer gro- 
Ben Zahl verschiedener Raiume vorliegt, ist es még- 
lich, nach einer Korrelation zwischen den subjek- © 
tiven Bewertungen und der GréBe Fy zu fragen. 
Dabei miiBte es sich herausstellen, ob Fy ein brauch- 
bares MaB fiir die ,,Diffusion“ darstellt, worunter 
in diesem Falle die-Homogenitit des Schallfeldes 
zu-verstehen wire. 
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Abb. 1. Frequenzgang eines Raumes von 83 m*; 1300--- 
1400 Hz; Ordinatenteilung 5 db. 


Als erster Schritt wurde das MeBverfahren fest- 
gelegt. Das System eines Trichter-Lautsprechers 
diente als quasipunktformige Schallquelle; seine 
verhaltnismaBig groBe mechanische Impedanz ge- 
wahrleistet zudem, da die Strahlung von der aku- 
stischen Belastung praktisch unabhingig wird; 
ferner ist auch der Frequenzgang eines solchen Sy- 
stemes verhaltnismaBig flach. 

Zunachst wurde der Einflu8 des Abstandes Laut- 
sprecher-Mikrophon untersucht, wozu in zwei Raéu- 
men auf bestimmten Me8strecken Fy als Funktion 
des Abstandes gemessen wurde. Der auf das Mikro- 
phon fallende Schall setzt sich aus zwei Anteilen 
zusammen : 

a) aus dem direkten Schall mit der Energiedichte 


E,: 


R= Ee P = Leistung der Schallquelle, 
Oar die d = Abstand von der Schall- 
quelle, 


c = Schallgeschwindigkeit. 


b) aus dem indirekten Schall mit der Energie- 
dichte E, (nach SaBrne): 
E, = 28 ; A=Zai S; =Totale Absorption des 

Raumes. 

In einem Raum entspricht daher jedem Abstand 
d zwischen Mikrophon und Schallquelle ein be- 
stimmtes Verhiltnis des indirekten Schalles zum 
direkten Schall: 
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Abb. 2. Frequenzgang-Schwankungen Fy in Funktion des 


Abstandes Lautsprecher-Mikrophon; 
(1) Raumvolumen 330 m*, Frequenzbereich 375 Hz, 


(2) € 330 m3, " 1650 Hz, 
(3) is 1100 m5, e 375 Hz, 
(4) 44 1100 m%, ns 1650 Hz, 
(5) s 1100 m’, e 375 Hz, 


6) $3 1100 m§, es 1650 Hz. 
(3) und (5), bzw. (4) und (6) wurden an verschiedenen 
Orten gemessen. 


In der Abb. 2 sind die gemessenen Werte von F'y 
als Funktion dieses Verhiltnisses aufgetragen. Es 
geht daraus hervor, da8 fir E,/E,>5 Fy praktisch 
konstant bleibt, so daB fiir alle Messungen von Fy 
d wie folgt gewahlt wurde: 

d=YV01 A. 

Ferner wurde die MeBgenauigkeit der Fy -Mes- 
sungen untersucht; dazu wurden in zwei Raéumen 
je drei MefBstrecken festgelegt, fiir 
die die Messungen jeweils fiinfmal 
wiederholt wurden. Die gréBte Ab- 
weichung zwischen diesen sechs 
MeBreihen betrug -- 3%. Endlich 
wurde der Einfluf des Ortes der_ 
MeBstrecke gepriift, wobei in einem 
Raum von 330 m® zwolf verschie- 
dene Lagen ausgemessen wurden; 
dabei ergab sich eine mittlere Ab- 
weichung vom algebraischen Mittel- 
wert von + 10%. Diese verhiltnis- 
maiBig groBe Streuung macht es 
notwendig, in jedem Raum an ver- 
schiedenen Stellen zu messen; die 
vorliegenden Untersuchungen um- 
fassen jeweils Messungen an sechs 
Orten. 

Fiir die Berechnung der Mittel- 
werte. einzelner Frequenzbereiche 
wurde die von Bort und Roop vor- 
geschlagene Mittelwertsbildung [4] 
iibernommen, wobei mit abnehmendem Gewicht die 
F,-Werte der Nachbarfrequenzen beriicksichtigt 
werden. Durch Verbinden der Mittelpunkte zwischen 
je zwei Punkten der urspriinglichen Kurve wird 
eine neue Kurve erhalten (Abb. 3); durch sechs- 


maliges Wiederholen dieses Verfahrens wird eine 
genugende Glattung des gewohnlich sehr unregel- 
mafigen Kurvenverlaufs erzielt, wobei dies auch 
algebraisch gemacht werden kann: 

Zu sieben aufeinander folgenden Frequenzinter- 
vallen v, bis v, gehdren die Werte Fy,...Fy,. Der 
Mittelwert (fiir Fy,) berechnet sich nun unter Be- 
rucksichtigung der beidseitig benachbarten F',, bis 
Fy, wie folgt: 

Ey 


Fy, + 6Fyy+ 15 Fy, +20Fy4+ 15 Fy; + 6Fyot+Fy, 
GOES 5 alti vu 


Abb. 3. Mittelwertbildung der Frequenzgang-Schwankungen ; 
° o gemessene Werte, 
Q)==-=---- o Mittelwertskurve. 


Abb. 4. Hérspielstudio (83 m*) mit polyzylindrischen Diffusoren. 


Die Gewichtsfaktoren, mit denen die Nachbar- 
frequenzen bewertet werden, sind also die Bino- 
minalkoeffizienten. 

Mit diesem Verfahren wurden zahlreiche Riume 
untersucht, deren Volumina zwischen 37 m* und 


254 


20000 m3 lagen. Aus der groBen Anzahl der Meb- 
ergebnisse werden hier nur zwei Beispiele angefuhrt: 

In einem Hoérspielstudio (83 m’) wurde Fy vor 
und nach dem Anbringen der polyzylindrischen 
Diffusoren gemessen (Abb. 4). Die tiber den ganzen 
Frequenzbereich ausgewertete Grofe Fy fiir beide 
Falle zeigt Abb. 5. Es geht daraus hervor, dai Fy 


its 
be 
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Abb. 5. Frequenzgang-Schwankungen I’, eines Hérspiel- 
studios (83 m*); 
se 22----- ohne polyzylindrische Diffusoren, 

mit polyzylindrischen Diffusoren. 


unter dem EinfluB der eingebrachten Diffusoren wohl 
bei den tiefen und mittleren Frequenzen deutlich 
sinkt, beiden hohen Frequenzen aber kaum. In einem 
Konzertstudio (2230 m’, Abb. 6) hat sich dagegen 
durch den Einbau sehr kraftig profilierter poly- 
zylindrischer Diffusoren die GréBe F, nur unwesent- 
lich verandert (Abb. 7), trotzdem der subjektive 
Erfolg sehr deutlich war. 

Da das Auswerten der GroBe Fy tiber den ganzen 
Frequenzbereich sehr mithsam ist, beschrankte man 
sich auf zwei typische Gebiete. Der Mittelwert des 


Abb. 6. Orchesterstudio (2230 m%) mit polyzylindrischen Diffusoren; vier 
weitere sind im Bilde nicht sichthar, 
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Abb. 7. Frequenzgang-Schwankungen J’, eines Orchester- 
studios (2230 m§); 
aeceeeeee ohne polyzylindrische Diffusoren, 

mit polyzylindrischen Diffusoren. 


einen betrigt 375 Hz, des andern 1650 Hz; dabei 
wird keine Riicksicht auf die RaumgroBe genom- 
men, diese Werte bestimmen sich vielmehr aus der 
Art der Schallquellen, die ja tiberall die gleichen 
sind. Bei der unteren gewihlten Frequenz liegt das 


F, 
3 


[dbHz] 2 


Abb. 8. Frequenzgang-Schwankungen Fy in Funktion des 
Raumvolumens; 
a--0=---- Frequenzbereich 375 Hz, 
Frequenzbereich 1650 Hz. 


Maximum der Energie der mensch- 
lichen Sprache [6], sowie sehr vieler 
Musikinstrumente wie z. B. Kla- 
vier, viele Blasinstrumente, Streich- 
instrumente usw. [7]. Die obere 
Frequenz liegt mit 1650 Hz etwas 
oberhalb der Mitte des ganzen Fre- 
quenzbereiches; dieses Gebiet ist 
sowohl fiir die Verstandlichkeit von 
Sprache als auch fiir viele harmo- 
nische Teilténe der Musikinstru- 
mente wichtig. Vom subjektiven | 
Standpunkt aus lé8t sich auch 
sagen, da der Bereich um 375 Hz 
wichtig fiir die ,,Tonfiille“‘ ist, wah- 
rend der Bereich.um 1650 Hz mab- 
gebend fiir die ,,Klarheit einer 
Ubertragung ist. Frequenzen ober- 
halb etwa 3000 Hz sind in diesem 
Zusammenhang nicht mehr so in- 


W. FURRER und A. LAUBER: DIE DIFFUSION IN DER RAUMAKUSTIK 


teressant, da hier die Richtwirkungen von Schall- 
quelle und Empfanger sowie ihre unmittelbare 
Umgebung eine gréBere Rolle spielen als irgend- 
eine mittlere Raumeigenschaft. 

Abb. 8 zeigt die fiir die beiden erwihnten Frequenz- 
bereiche gemessenen Frequenzgang-Schwankungen. 
Das Bild ist unklar, es ist keine Systematik ersicht- 
lich und insbesondere treten die beiden subjektiv 
als schwierig bekannten Raume: ein Studio von 
480 m? und ein Konzertsaal von 9500 m? (mit 
,,Pfeil’ bezeichnet) nicht besonders hervor. Man 
mu daraus schlieBen, daB die Frequenzgang- 
Schwankungen selbst offenbar nicht direkt ein 
brauchbares raumakustisches Kriterium darstellen. 
Hier setzt nun die folgende Uberlegung ein: Wenn 
das anzustrebende Ziel ein moéglichst homogenes 
Schallfeld ist, dann diirfen die 
Frequenzgang-Schwankungen wohl 
sehr viele kleine ,, Berge“ aufweisen, 
nicht aber nur wenige grofe. Der 
Unterschied zwischen diesen beiden 
Zustinden wird aber durch die 
GroBe Fy nicht erfaBt. Er wird nur 
dann erhalten, wenn wir die mittlere 
Hohe aller ,, Berge“ beriicksichtigen, 
d. h. wenn wir die GroBe Fy durch 
die Anzahl n der ,,Berge“ divi- 
dieren. Auf diese Weise erhalten 
wir die neue GroBe 


Dy=Fy/n, 


worin n die Anzahl der Schalldruck- 
maxima bedeutet. 

Wie die Abb. 9 zeigt, erhalt man 
damit ein vollstandig veriandertes 
Bild, Man sieht zunichst sehr deut- 
lich, daB Dy bei 1650 Hz nur etwa 
halb so groB ist wie bei 375 Hz. Dies 
ist ohne weiteres einleuchtend, da 
es in jedem Raum sehr viel mehr Un- 
ebenheiten gibt, die in der GroBen- 
ordnung der Wellenlange von 20 cm 
(1650 Hz) liegen, als solche der 
Wellenlinge 90cm (375 Hz). Fer- 
ner heben sich nun das subjektiv 
als ,,schwierig‘‘ bekannte Studio von 
480 m3, sowie der ,,schlechte‘‘ Kon- 
zertsaal von 9500 m8 (mit ,, Pfeil‘ be- 
zeichnet) deutlich von den iibrigen 
ab. Der gré8te untersuchte Raum 
(20000m%) war eine Barockkirche 
aus dem 17, Jahrhundert (Abb. 10), 
die eine ausgesprochen diffuse Re- 
flexion aufweist und deren aku- 


stische Verhaltnisse als sehr gut bekannt sind. Weiter 
kann man. auch an den frither von Fret gemessenen 
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Abb. 9. Diffusion Dy in Funktion desjRaumvolumens; 


wareneoe= Frequenzbereich 375 Hz, 
Frequenzbereich 1650 Hz. 


Abb. 10. Barockkirche mit diffuser Reflexion. 


256 


Resultaten zeigen, daB durch die dort untersuchte 
Kassettierung die GroBe Fy (dort in Funktion des 
Abstandes gemessen) zunimmt und erst die GréBe 
Dy den richtigen Einflu8 der Kassetten auf die 
Homogenitit des Schallfeldes erkennen laBt. 

Die so definierte GroBe Dy ist also um so kleiner, 
je homogener das Schallfeld und je groBer die Dif- 
fusion ist; wir méchten sie daher als ein Mai fiir 
die ,,Diffusion‘‘ bezeichnen. Dabei ist es klar, da 
Dy nicht ausschlieBlich von der Form des betreffen- 
den Raumes abhangt, sondern auch noch von an- 
dern Faktoren, so vor allem von der Absorption. 
In einem vollstandig absorbierenden Raum wiirde 
ja Fy und damit auch Dy gegen 0 tendieren. Dieser 
BinfluB konnte noch nicht im einzelnen untersucht 
werden, doch finden sich in der bereits zitierten 
Arbeit von Wente [3] schon Angaben tber den 
EinfluB der Absorption auf die Frequenzgang- 
Schwankungen. Die von WENTE gemessene Kurve 
14Bt sich durch den Ausdruck Fy =ko—18 analytisch 
approximieren. Es geht daraus hervor, daf} der Kin- 
fluB der Absorption mit steigendem mittlerem Ab- 
sorptionskoeffizienten rasch abnimmt. Da alle hier 
untersuchten Raiume einen gréferen mittleren Ab- 
sorptionskoeffizienten als 0,1, in sieben Fallen sogar 
als 0,2 aufweisen, wurde auf eine Beriicksichtigung 
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verzichtet, da sie das Bild nur sehr wenig verandern 
wiirde. Die dargestellte Diffusion Dy ist daher als 
,,Betriebsdiffusion zu betrachten, die neben dem 
Einflu& der Raumform zu einem viel kleineren Teil 
auch noch den Einflu8 der Raumabsorption in sich 
schlieBt. 

Zusammenfassend kann festgestellt werden, daf 
die Diffusion Dy, so wie sie oben definiert ist, ein 
neues und-anscheinend brauchbares Kriterium fir 
die Homogenitat des Schallfeldes und damit fir 
das akustische Verhalten eines Raumes bildet. Sie 
muB8 in akustisch einwandfreien Raumen unterhalb 


der beiden Kurven der Abb. 9 liegen. 
(Eingegangen am 28. April 1952.) 
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Sommaire 


Dans le présent travail, nous nous occuperons des échos des ondes acoustiques produits par la 
stratification de latmosphére en nous fondant sur une méthode de BREMMER et sur une mé- 


thode de SCHELKUNOFF. 


Dans cette étude l’influence des termes qui troublent l’analogie avec les ondes électriques, selon 


le travail précédent [1], apparait clairement. 


Zusammensetzung 


Die vorliegende Arbeit beschaftigt sich mit der Reflexion akustischer Wellen infolge der Schich- 
tung der Atmosphare und stiitzt sich dabei auf eine Methode von BREMMER und eine von SCHEL - 


KUNOFF. 


Der EinfluB der Terme, die die Analogie zu den elektrischen Wellen stéren (siehe die voran- 


gehende Arbeit [1]), kommt klar zum Ausdruck. 


Summary 


This paper deals with the echoes of sound waves produced by a stratified atmosphere, by the 


methods of BREMMER and SCHELKUNOFF. 


In this study the influence of those terms which spoil the analogy with electric waves (cf. the 


preceding paper [1]) becomes clearly apparent. 
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1. Les équations différentielles pour une onde acous- 
tique dans un milieu stratifié par la gravitation 


Dans ce qui suit, nous supposons une atmosphére 
au repos qui est stratifiée le long de l’axe des z posi- 
tifs par suite de la gravitation et dont la pression, 
la densité et la température dépendent également 
de z. 

En négligeant les termes non linéaires, on obtient 
dans ce cas les équations!: 


— P= jone + 8.5, (1) 
dv jo lide, 

pas 8, 2 

dz poy ak Pp ad 2) 


ou g est l’accélération de la pesanteur. 

On obtient ces équations en posant: 

1. Les équations d’Eunmr de l’écoulement des flui- 
des, ot l’on tient compte de la gravitation qui 
agit vers le bas, ce qui se manifeste dans le terme 
qui contient g, 

2. Péquation de la continuité, 

3. ’équation d’état adiabatique. 

Il y a des cas ot l’on peut négliger les termes qui 
ne contiennent pas w dans le membre de droite des 
équations (1) et (2), en particulier quand w est 
suffisamment grand, mais de fagon générale, cela 
nest pas possible. Nous voyons done que l’analogie 
qui existe dans un milieu homogéne se trouve per- 
turbée dans un milieu hétérogéne par ces deux 
termes. 

Cela a pour conséquence que dans un milieu stra- 
tifié d’une maniére continue par la gravitation, les 
ondes acoustiques subissent toujours une réflexion 
interne. Nous n’en donnons pas la preuve ici, le 
lecteur la trouvera dans la référence [1]. 


2. La solution a Vaide d’une méthode W K B, d’aprés 
Bremmer 


Dans le cas des ondes électromagnétiques planes 
qui se propagent dans le plan yz dans une atmo- 
sphére stratifiée ott « est fonction de la coordonnée 
z, Bremer [2], [3], [4] a élaboré une méthode 
particuliérement claire et élégante qui est une va- 
riante de la méthode connue W K B. BremMer 
suppose un champ E polarisé dans le plan yz et une 
seule composante H,. Etant donné l’analogie for- 
melle, mais toutefois incompléte qui existe avec les 
équations qui régissent la propagation des ondes 
électriques, nous allons traiter le probléme acous- 
tique par la méthode de BremmER qui ne nécessite 
que de petites modifications. 


1 DP» Por Por Vi Voir chapitre 4. 


bo 
Or 
«J 


3. Les équations de Maxwell pour les ondes planes 
dans un milieu stratifié — leur résolution d’aprés 
Bremmer 


Pour étre compris par la suite et éviter au lecteur 
de se reporter 4 des ouvrages de références difficiles 
a aborder, nous allons exposer briévement la mé- 
thode de Bremmer. Contrairement 4 BREMMER, nous 
supposons que le champ H est polarisé dans le plan 
d’incidence, le champ E lui étant perpendiculaire 
(composante EF, seulement). 

Dans le systéme de coordonnées de la Fig. 1 les 
équations de MaxweE Lt s’écrivent: 


: Oe Oued 
y Joce(s)Ea= oy Ath ma ? (3) 
oP De 
A jopoule\Hy=—", (4) 
JE 
-Z . z)H, = zt ~ 
ee jou u(2) by (5) 


avec et J! comme fonction du temps, 


il 
“0 ~ An- 9-109 [By], te An LOS? [aim]. 


Introduisons l’angle d’incidence + qui, pour z=0, 
prend Ja valeur ta. On a par suite: 
IEE a AATEC 


D’apres ae on écrit maintenant les trois 
intensités E,, H,, H:, du champ en utilisant le 


=) upu(z ) sin + ( )sinta. (6) 


symbole: 
k=k( (2) =0V eqe(e 2) urote(2 (7) 
% 1 
eee dc ae (8) 
ene Z 


ot Z est ’impédance caractéristique, par suite 


exp (—jky sin 7) 


E,;=— iVubas [x(2) exp (—ifk cost dz)-+ 
Zz 

+79(2)-exp (+j/kcos+dz)], (9) 
0 


H, = exp(—jky sint)- VY cost: 


[x (2): exp (—j [k cost dz) + 
0 


+ y2(z): exp (+j /k cos dz)], 
0 


ish 

H,= exp (—jky sins) l/ ae 

* [y(z): exp (—j/k cos +dz) + 
0 


+Y_(z)- exp Gfk cos + dz) }. (11) 
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Nous remarquons que les équations de MaxwEL1 
ne donnent d’abord que les relations entre les trois 
termes avant les crochets. On peut les multiplier 
par une fonction quelconque; et en introduisant les 
fonctions y,(z) et y2(z) dans les équations de Max- 
WELL on détermine ensuite les amplitudes des on- 
des dirigées vers le haut ou vers le bas sous une 
forme que ]’on peut donner immédiatement. Si l’on 
avait multiplié auparavant par un autre facteur 
f(z) on aurait eu simplement x1,2(z)/f(z) 

En introduisant ceci dans les équations de Max- 
WELL, BREMMER trouve pour y, et 7%, le systéme 
d’équations intégrales: 


Zz Z 
=| ale) ‘exp if cos tdz) os (In y Yeost) dz, 
(12) 


Zz 


X2(2) =x 1(2) -exp( (24 [eossae) 4, 


+ co 


(Iny Y cost) dz 
(13) 


que l’on peut résoudre facilement par approxi- 
mations successives et l’on obtient les réflexions des 
différents ordres; on choisit la premiére approxi- 
mation de facgon telle que pour 


d ames 
q, MVY cost=0, ks 


Ce dernier quotient différentiel représente, comme 
Bremmer l|’a montré, le facteur de réflexion du 
milieu stratifié par unité de longueur. 


4. Les équations acoustiques analogues et le terme 
correctif supplémentaire 


Nous écrivons maintenant les équations différen- 
tielles limitées aux termes linéaires qui régissent la 
propagation des ondes acoustiques dans une at- 
mosphére stratifiée soumise a l’accélération de la 
pesanteur g. Les notations employées sont les 
suivantes : 


Pol? z) = pression de l’atmosphére au repos, 

p (z) = amplitude de la pression acoustique, 

29 (z) = densité de l’atmosphére au repos, 

Vz,y = x,y-composants de la célérité (vitesse 
dune particule d’air), 


exp[—jk(z \y sin +(z)] 


pe {zs (2)-exp [i eos (=) de] +2) exp li/ cos (2)de]| 
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g = accélération de la pesanteur, 

y = rapport des chaleurs spécifiques=¢p/ey—=1,4. — 

Les équations de la propagation du son ont pour — 
base: les équations d’Eurer de ’hydrodynamique, 
Véquation de continuité et léquation d’état adia- 
batique. On a: 


10D Ov, dV; 1 
jop avy ae _ 100.9 tel aa 
ae ae 
op 
Bh aaa ’ oe 
OP ~ eB 
J@ 09 Uz Sprains — ee om (16) 


Poe 

Si nous faisons correspondre maintenant dans 
le schéma suivant les grandeurs électriques de la 
premiére ligne aux grandeurs acoustiques de la 
seconde ligne, nous voyons dans les équations (14) 
et (16) une analogie qui n’est perturbée que par les 
termes soulignés a l’extréme droite. 
Grandeurs électriques 


ey e(z), Ugt(z), Ex, Hz, 
Grandeurs acoustiques 
1 = 
Tht 


Remarquons que cette analogie connue n’est pas 
la seule possible. Si nous négligeons les termes 
soulignés, dans les équations (14) et (16) les relations 
de Bremmer donnent la solution de notre probléme, 
mais elles sont écrites cette fois avec d’autres 
symboles. 

Nous allons étudier maintenant l’influence de ces 
termes; introduisons encore les notations suivantes : 


Yac=1/Zac= 1/V yPo(2)po(2) » (17) 
P02) gin a(z)= 0 1/0). sine 
noes sin t(z) = a sin ta (18) 


ot +(z) et ta sont les angles d’incidence aux altitudes 
z et 0 respectivement. Nous posons d’autre part: 


(19) 


et nous écrivons de fagon analogue a (9), (10), et (11) 


(20) 


vz = exp [—jk(z)y sin +(z) ]/ Yac cos +(2) | —x,(z)- exp [—ifk cos7(z) dz] +7(z)- exp (j/k Cos t(z) dz] | , (21) 
0 0 


— ty= exp [— jk(z)y sinc(z ao Xi (z )-exp [ih cos (2) dz} +22) -exp [j/k cos=(2) de} . 


(22) 


} 
| 


dy, (2) 


dxo(z ) 


Ce systéme a la forme: 


dy. oR ~ 

au = 4y1(2)X1 (2) + 44s (2) Xa (2), (25) 
ae 

at” = 4; 4(2)%4(2) + @22(2) X22) (26) 


que nous allons mettre sous une forme plus pratique 
de fagon qu’il y ait toujours une des fonctions 
cherchées qui figure 4 droite. A cet effet, nous 
écrivons: 


X= PrrX*1, Xo= Pore (27) 
par suite: 
i 
Xa = Pal 1 + Pi 1 = Gy (2) Pi x*1 + G2(%) Po z*2, 
| (28) 
= Pe! Xo + Po X*o’ = Gy1(2) Py ~*1 + Go9(2)Po x*o- 
(29) 


Nous fixons p, et p, de facon a tirer les termes 
soulignés de (28) et (29) 


PX") = 41(2) Pi x*13 Py= G- -exp [fan )ds| » (30) 


Po’ h*2= 4992) P2%*2; Py = C,- exp es (2) dz| . (31) 
Il ressort du systéme (23), (24) que: 
Bem (| dey Se ny 
a= Anl)= 5 (5 ge +E) (32) 
si bien que: 
ie, ae 
rceoe[S [tate] o 
Pi=P2 exp 9 ‘ a dz P pe (33) 
et 
x*)'= x*,° exp (2j {k cos dz) . 
da Rey I _ deo £9 
. E In V Yac cos t+ aes hag £m) , (34) 
= hy-exp(—2i /keosde)- 
a eae es EL 
. la: Iny Yac cos t + 2 ( ” ds of see) > (35) 


Nous voyons done que nous obtenons avec les 
termes supplémentaires un facteur d’amplitude 
purement réel pour nos grandeurs x, et x2. 
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Nous obtenons ensuite en introduisant dans (14), (1é 


Ieyel ° 
asi. = YZ»(z)- exp 23 /k eos de) E — In V Yac cos +(z) + 2 ( si +8)) AGE 


# tr d 
ee = in(#)-exp(—2i/keosda)| 5, In y Yac cos t(z) + = ( 
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5), et (16), le systeme: 
do, 1 (; e+ 8h) 
Po 2 \e dz” por’ 
(23) 
os bs _ de Bee) +%Xo(z): 1 (; deg 520) 
Po dz § pox ae 2X eg de” poy 
(24) 


Ecrivons maintenant ce systéme de facon ana- 
logue a (12), (13): 


L ee 2 CXD (: 2} [hos +dz) - 


=e 


L fl” dey . Be i a 
a a Iny Nae COST -+ pal ae . dz ‘ai Po dz, (36) 


ne ie 1 exp ( (—2j /k cost dz) - 
0 


-| gg IRV Yee cos + : (- eo 


+80) dz. (37) 

Comme Bremuer l’a montré, on peut résoudre 
ce systéme par approximations successives, un 
terme en 7, représente toujours la réflexion pro- 
venant du terme en x, précédent et vice versa. Nous 
ne calculerons que la premiére réflexion qui est due 
a une onde oblique se propageant vers le haut avec 
x40 =1 ot 7, est la premiére approximation 
d’indice 0. Par suite, on a: 


1 PY 1 dey , 8 
4 (0) == a peas : 
‘a oP 2 | ( Po dz fe a) z 


et 


Po dz 


(38) 


= Zz 


LL doa mge : \ 
*,(0)— | exp] — = “Po , 8% Pee rel, 
B bas 1 exp| 9 Il at 4jh cos) d | 


d Ty dep eo : 
: E Iny Yac cost + 9 & Ps + 60) dz (39) 
il do, Bon 
D’ott X2')( ) exp 2 /( - . as gta) d 
; ees dey SP :) 
: _ —4jk cost) dz| - 
/ exp | 2. ( Biel jk co 
+ co 
d a ] (. i deo &Po 
. as Iny Yac cost + Bos Sale +880) dz. (40) 


A partir de 7*,() que Von introduit dans (40), 
on peut calculer ~*,“) et ainsi de suite. Nous n’en- 
trerons pas ici dans le détail des calculs. Nous nous 
contenterons de mentionner l’influence du terme de 
(38) qui fausse l’analogie avec les ondes électriques. 
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5. La méthode de Schelkunoff 


dans le cas de l’incidence normale 


Changeons maintenant (cf. [7], [8], [9]) le sens de 
notre systéme de coordonnées et choisissons z néga- 
tif vers le haut. Nous supposons que l’origine est 
au-dessus de la couche considérée. 

De méme que dans [7] et [8], nous 

, allons étudier la propagation d’une 

onde venant du bas c’est-a-dire de 

y grandes valeurs positives de z et 

traversant la couche de z=h a 

z =0. Nous avions supposé dans [7] 

et [8] que l’onde pour z négatif, 

c’est-A-dire au-dessus de la couche, 

se propageait vers le haut dans un milieu homo- 

géne sans réflexion. Dans le cas des ondes acous- 

tiques, nous supposons que pour z = 9, il y a une 

onde qui sort de la couche et une onde réfléchie 

correspondante dont le facteur de réflexion com- 
plexe est r. 


tad 


2 
Fig. 2 


6. La propagation dans la couche 


Dans le cas de l’incidence normale, les équations 
(14), (15) et (16) se réduisant a: 


dp © gpgls)? 20% 

rs a Pov *P + J@Pov, (41) 

dv ja clon 
recat acca PE saptO 42 
dz Por ey d ) 

nous écrivons, d’aprés ScHELKUNOFF [6]: 

29(2) = Poo + Po(2); (43) 

I I Tl | 
= 44 
PoY @) ae be @) -) 


it 1 
ou —— (z) représente —— en fonction de z. Nous sé- 
PoY¥ PoY 


parons donc les grandeurs ¢» oe en deux parties: 
) 


a) une partie constante correspondant a la valeur 
pour z = 0; 

b) une partie variable supplémentaire. 
Pour plus de clarté, mettons le systéme (41) et 
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(49) 


Veo (Povo = K. 


Considérons comme ScuetKunorr, la partie du 
systéme (45), (46) & gauche du trait vertical comme 
un systéme homogéne, et la partie 4 droite comme 
le terme perturbateur de Vhomogénéité. Nous avons 
pour le systéme homogéne, la solution connue de 
Véquation des télégraphistes dans le cas d’un 
cable homogéne. Si nous résolvons le systéme non 
homogéne avec les termes de perturbation en fai- 
sant varier les constantes, nous obtenons pour p(z) 
et v(z) le systéme suivant d’équations intégrales: 


(50) 


[ple)—io / v(t) Go(t) coshg(s—t) de — 


ty 
x 
I 


-- ok [70 las ()| sin ky (s—t) dt + 


(51) 
= jKo/v(6)m(t) sin hy (e—tydt= 
# [ p(t)m(t) cos k, (z —0) at; 
v2) =[h—- / o(C) a(t) aint hy ae 
Bell bing i 
ees | B® | (6) ]cos het) 4&— 
2 (52) 


— v(t) n(%) cos ko(z —%) dt + 


0 
iL te tae, we 
1K | p(t) m(£) sin k(z—t) a 


ow [p(z)]y et [v(z)]p sont, comme nous l’avons déja 
dit, les solutions du systéme homogéne. Posons: 


(42) sous la forme: [p(z)]o = Jhez + re oho P (53) 
dp _. = 8(Po0+ Po) ,— . 1 jkgz , To ike 
Tak ica aiehie Vela +e ray P48) Web — Oe 
__do ah jo(-) = jo Lg | =e 1 oa A(P99+ Po) ow r, est le facteur de réflexion, c’est-a-dire 
da ps amar o : PoY P Poo + Po dz le rapport complexe de l’onde réfléchie & 
(46) onde transmise au point z = 0. 


Pour plus de commodité, posons: 
& (P00 + Po) 7 
PoY(2) 


1 De A(Po9 + Po) 
Poo + Po dz 


m(z), 


(47) 


=n(z), (48) 


7. La transformation du facteur de réflexion par la 
couche 

Si, au bord supérieur de la couche (z = 9), il 

existe un certain facteur de réflexion r, et que d’aprés 

(55), (56) on ait un facteur de réflexion rj, au bord 
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inférieur de la couche, la formule r;, =f (r, a =) 
0 


représente la transformation du facteur de ré- 
flexion par la couche. Posons ry = 0, nous obtenons 
le facteur de réflexion df a la couche seule. Comme, 
dans le cas considéré, l’onde incidente se propage 


dans la direction — z et londe réfléchie dans la 
direction +z, nous obtenons pour r;, 
P(h) + K(h) v(h) 

ree ra a 2 (55) 
are p(h) Biiath) 

= Veo(h) pal) y (56) 


qui est ’impédance caractéristique pour z = h. 


are oh < *< | a 


0) 


+jsin ky —9)| Kijo(5 (0) — n(¢) —jo(o) K? 


— cos ky(h— ) m(%) | dt + 


+ “jm F 4 | + 
0 


jo © P9(S) 
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Comme notre expression est assez longue, nous 


désignons le numérateur par M, le dénominateur 


par N et nous obtenons: 


r(h) = M/N. (57) 


Nous devons résoudre maintenant les équations 
intégrales (51) et (52), nous nous contenterons de la 
premiére approximation et nous poserons simple- 
ment pour les termes p(z) et v(z) sous le signe somme: 


[p(z)]o et [v(z)]o- 


Nous obtenons aprés avoir groupé les termes 
analogues: 


h - 
+ [| cos ky h— 0B + jo(E (0) K(h) + n(¢ 7 [+ 


pee ee th) 


BIEK, 


-m(¢)|— 


(58) 


Kase” 


: h 
+ “f eo iho’ cosk,(h—t) |- 
0 


ah As 
+ jsink,(h— | Kio(—- (©) + n(t) + jopo(t) ye singles 


-m(°)] — cos ky(h + smc} g 


Nall [1 4 Ae] 4 
+ / okt {00s ry(h—z) [Rt —j o (2 (2) Kem) —mce) A] + 
0 
+ jsinky(h—2) | Kyo (2) —n() + jek) RE m(®)| — 


— cos ky(h + ¢)m(%) dt + 


+ | oe akeh |i 


h 


+ [ oe ikeb cos ky(h — %) [Ri 


ca 


+ jin kyl—2)| Kujo (1. (@)) + + n(0)—joso() 


o (3) Km) +m 


(59) 


+ 


K(h) ae 
Ss Senay 


m()| — cosk,(h + %)m(%) |) : 
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La discussion de ces formules est relativement 
simple si l’on se souvient de celle qui a été faite dans 
le cas d’ondes électriques dans une atmosphére ot 
« varie seul [8]. 

Naturellement, aprés substitution de ¢(¢) par 


Pes (0 , ces deux derniéres formules renferment 
0 


tous les termes que nous avions auparavant. A cela 
s’ajoute un terme analogue avec e,(¢) quin’apparais- 
sait pas parce que nous avions supposé la perméa- 
bilité » constante. D’autre part, nous avons une 
série de termes qui perturbent l’analogie entre les 
ondes électriques et les ondes acoustiques dans un 
milieu stratifié. Si nous supprimons ces derniers, 
nous obtenons, aprés avoir remplacé les grandeurs 
acoustiques par les grandeurs électriques correspon- 
dantes, la propagation d’ondes électriques dans une 
couche ou p. et © varient. 

Nous obtenons directement l’influence des termes 
m(C) et n(¢) qui sont négligés dans la plupart des 
travaux. 

Gomme il ressort de (58), (59) le facteur de ré- 
flexion a la forme: 


A+B 


r= 4D: 


(60) 

Ce qui est nécessaire physiquement; si nous po- 
sons r= 0 nous obtenons la réflexion due a la couche 
seule. 

Nous avons vu plus haut, dans la méthode de 
Bremmer, l’influence des termes qui troublent l’ana- 
logie avec le cas des ondes électriques; dans le cas 
de la méthode de ScurLkunorr les grandeurs addi- 
tives m et n qui reproduisent cette influence se pré- 
sentent dans les intégrales. 

A part cela, la discussion des formules (58) et (59) 
n’est pas difficile mais ennuyeuse. Elle se fait d’une 
maniére tout a fait analogue 4 celle qui est exposée 
en [9] et en [10]. Elle est seulement un peu étendue. 
Pour éviter au lecteur d’avoir & se reporter 4 ces 
références, nous en donnerons ici les lignes princi- 
pales. 

On reconnait aussit6t qu’il existe un cas limite. 


Quand p, (%) et Fe ) sont identiquement nuls 
0 


entre z=0 et z=hA et quand pour z=h la ré- 
sistance caractéristique saute brusquement de K, 4 
K(h), quand en outre r, = 0, il résulte alors de (59) 
et (58), étant donné que toutes les intégrales s’an- 
nulent 
Ko Kil) 
K,+ K(h) 
c’est-d-dire le coefficient de réflexion de la transition 
brusquement discontinue. 


(61) 
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En outre, nous voyons dans les grandeurs m et n 
directement l’influence des membres qui troublent 
V'analogie avec le cas électrique. 

Si nous désignons les parties des formules (58) et 
(59) se trouvant entre parenthéses par K nous avons 
alors des intégrales de la forme: 


h 
 jkot | cosko(h =F &) 
[e ) ne (hk —t) | Bie | 


qui peut se séparer en intégrales: 


h 
[ x-ci806. tik OFE) ge. 
) 


Selon le signe k)¢ dans |’exposant peut s’annuler 
ou se doubler: 
Tl reste par exemple: 


1 jpo(5) 


h 
i Qk joo,(%) ? 
ony [e OF Ky dv. 
0 


La premiére expression représente la surface sous 
la courbe @,(¢) =f(¢) de la couche, la deuxiéme 
expression représente la contribution des éléments 
de cette surface 4 l’influence du chemin parcouru 
par l’onde. On verra ainsi directement dans chaque 
cas se présentant en pratique, comment le facteur 
de réflexion s’obtient & partir de ’hétérogénéité 
existante. (Regu le 17 Janvier, 1952.) 
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DIE AKUSTISCHE STRAHLUNG 
DER RECHTECKIGEN KOLBENMEMBRAN 


Von HEINRICH STENZEL 


Fernmeldetechnisches Zentralamt Darmstadt 


Zusammenfassung 


Es werden allgemeine Formeln abgeleitet, die die Berechnung des Schallfeldes unmittelbar 
vor einer rechteckigen Membran gestatten. Das sich anschlieBende Schallfeld der Um- 
gebung wird durch graphische Integration bestimmt und durch zahlreiche Abbildungen erlau- 
tert. Dann werden allgemeine Formeln zur Berechnung der Strahlungsimpedanz entwickelt 
und die beiden Komponenten fiir verschiedene Seitenverhaltnisse des Rechtecks (a: }) = 1:1, 
a:b = 2:1,a:b6 = 5:1, a:b = 10:1) berechnet und graphisch dargestellt. 


Summary 


General formulae are obtained for calculating the sound field immediately in front of a 
vibrating rectangular membrane. The vicinal sound field is calculated by graphical inte- 
gration and illustrated by many diagrams. Then general formulae for the radiation impedance 
are derived; and finally its two components are calculated and graphically constructed for different 
ratios of the sides of the rectangle (a:b = 1:1,a:b6 = 2:1l,a:b =5:1,a:b= 10:1). 


Sommaire 


On a établi des formules générales qui permettent de calculer le champ sonore juste devant 
une membrane rectangulaire. Le champ sonore adjacent est déterminé par intégration 
graphique et son allure est montrée par de nombreuses figures. On établit ensuite des formules 
générales pour calculer ’impédance de rayonnement, et les deux composantes relatives & 
différents rapports des cétés du rectangle (a:b =1:1,a:b = 2:1,a:b = 5:1, a:b = 10:1) 
ont été calculées et représentées graphiquement. 


Allgemeines 

Wahrend die Untersuchung der akustischen Strah- 
lung einer kreisf6érmigen Kolbenmembran schon 
vor langer Zeit durch Rayteten in seinem klassi- 
schen Buche [1] in Angriff genommen und durch 
zahlreiche Arbeiten anderer Verfasser bis in die 
neueste Zeit [2] wesentlich geférdert wurde, ist die 
akustische Strahlung der rechteckigen Kolben- 
membran bisher wenig behandelt. Zuerst hat Rine- 
GER [3] allgemeine Betrachtungen tiber die recht- 
eckige Kolbenmembran angestellt und in zwei spe- 
ziellen Fallen (Seitenverhiltnis a:b =1:1 und 
a:b =1:7) die beiden Komponenten der Strah- 
lungsimpedanz in Kurvenform angegeben. Doch 
befriedigt das Ergebnis insofern wenig, da die zur 
Berechnung notwendigen Formeln nicht angegeben 
sind. Die spiteren Arbeiten beschriinken sich, abge- 
sehen von der Berechnung des Fernfeldes [4], auf 
die Untersuchung der einen Komponente der 
Strahlungsimpedanz [5], [6] (Strahlungswider- 
stand), die aus der. Formel fiir das Fernfeld abge- 
leitet wird. 

Das Ziel des Folgenden soll es sein, eine einiger- 
maBen vollstindige Ubersicht iiber die Strahlung 


der rechteckigen Kolbenmembran zu geben. Dazu 
sollen bei verschiedenem Seitenverhaltnis des Recht- 
ecks einmal die entsprechenden Schallfelder durch 
die Kurven konstanter Amplitude und zweitens die 
entsprechende Strahlungsimpedanz durch die bei- 
den Komponenten graphisch dargestellt werden. 


I. DAS SCHALLFELD IN DER UMGEBUNG 
DER MEMBRAN 


1. Die Berechnung der Schalldruckkomponenten un- 
mittelbar vor der Membran 

Wir legen eine in starrer Wand schwingende 
rechteckige Kolbenmembran A BC D (Abb. 1) zu- 
grunde, deren Schnelle durch w*=w,-e!! gegeben 
ist. Dabei bedeutet w, die (konstante) Geschwindig- 
keitsamplitude der Membran. Der von der strahlen- 
den Flaiche F in einem beliebigen Aufpunkt M er- 
zeugte Schalldruck setzt sich nach RAYLEIGH aus 
den einzelnen, von den Flichenelementen dF’ her- 
rihrenden, Anteilen zusammen. Wir driicken dies 
durch die Beziehung aus: 


pt=cow: e! (pg + ipm) 


=com-e(°' TX. p (1) 


mit 
: i fei* 
Pat+1pm= Seu fi aad dF, 
i (2) 


Pr Patpm ? tex = Pm/Pa.- 


Abb. 1. 
Zur Berechnung der Schalldruckam- 
x plitude (Projektionspunkt innerhalb 
der Membranflache). 


Die (dimensionslosen) Groen 
1 
ates ee Pma>= ; | Sar 
F 


bezeichnen wir als die (relativen) Schalldruckkom- 
ponenten und p als die (relative) Schalldruckampli- 
tude. 
Ferner bedeuten: 
c die Fortpflanzungsgeschwindigkeit des 
Schalles, 
co die Mediumdichte, 
k=2n/2 die zur Wellenlinge 2 gehérende Wellen- 
zahl, 
el den Zeitfaktor der fortschreitenden Welle, 
r den Abstand des Aufpunktes vom Flachen- 
element, 
Jo, J, die Besselsche Funktion nullter und erster 
Ordnung, 
Y,, Y; die Neumannsche Funktion nullter und 
erster Ordnung, 
H,, H, die Struvesche Funktion nullter und erster 
Ordnung. 


Liegt der Projektionspunkt O des Aufpunktes M 
innerhalb der Membranfliiche (Abb. 1), so zer- 
legen wir die Flache in vier Rechtecke, so daB 
F = OHDI + OIAE + OEBG + OGCH ist. Liegt 
der Projektionspunkt O’ des Aufpunktes M’ auBer- 
halb der Membran (Abb. 2), so zerlegen wir die 
Flache F in vier Rechtecke, so daB F = O’H’D’'l’ 
+ O'T'A’E’ — O'E'B’G’ — O'G’'C’'H’ ist. Der all- 
gemeine Fall setzt sich daher aus den besonderen 
Fallen zusammen, in denen der Aufpunkt senk- 
recht iber einer Ecke des strahlenden Recht- 
ecks liegt. Wir wollen nun zunichst den Schall- 
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druck unmittelbar vor der rechteckigen Kolben- 
membran untersuchen. Dazu geniigt es, sich auf 
den Fall zu beschranken, wo der Aufpunkt mit 
einer Ecke der Membran zusammenfallt. 


Abb. 2. Zur Berechnung der Schalldruck- 
amplitude (Projektionspunkt 
x auBerhalb der Membranflache). 


Die relativen Schalldruckkomponenten pa und pn, 
im Eckpunkt des strahlenden Rechtecks mit den 
Seiten x, y sind nach Gl. af bestimmt durch: 


‘ocik ye+7 
Pat ipm= 9(x, Ly) — Vey ReeRce dédy. (3) 


Setzen wir fiir den Zahler des Integranden die 
Reihenentwicklung ein, so folgt: 


k k3 k® 
o(%¥)= > P— 7.91 Mt pele 
(4) 
1 k2 kA 
+i toy a ee 
mit 


cy 
Pe =| [e+ ean. (m=—1,0,1,2,...) (5) 
0 0 


Im Fall m=2n, d. h. fiir den reellen Teil von 9(, y), 
entwickeln wir den Integranden in Gl. (5) nach dem 
Binom. Dann kénnen wir die Integration ausfithren 
und bekommen: 


P2n (x,y) —- 


(i) e Ln) ) . (6) 


een 
<i (Qnt1)' 3-(2n—1) (2n+1)-1 


Im Fall m=2n—1, d.h. fiir den imaginiren Teil 
von 9(x, y), ist es zweckmiBig, die Integration in 
Gl. (5) getrennt iiber die beiden Dreiecke ADC und 
ABC auszufiihren. Wir finden dann (Abb. 3): 


Pan— 1 (x, x)= 


:, 
[ dy / (G24 ya + | dé je (4 q2y'"*dy. (7) 
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Zur Auswertung von Gl. (7) benutzen wir die Gleichungen: 


x z 
1 1 n—/, 
nf (ath de a(x? + yt (2m — 1) ot fH ai) 
0 0 
x x D0 (64 
(2 n—2) [@ af ney dé = x(x? +2)" ah (2n Es 3) af ab ney ls dé 
0 0 
ue as ret y ab 
= — _ ae sy 
c x q 
4 [ (Bt WP) dé=x(x2+y*) + 37? A (62 + 4?) dé aig B 
i j Aaa Shem Ie RUT 
“4 we Abb. 3. 
2/@ + n2)'/2 de = x(x? + y?)'l at | (624+ nz)! d& Zur Berechnung des Integrals iiber 
é die Flache des Rechtecks. 


fe va np)—!2 dé = log( #+V5 + #) 
7 
0 


2n—1 (2n—1)(2n—3) 


Indem wir die erste, zweite, dritte usw. dieser Gleichungen mit 1, sae iat 4 usw. 
multiplizieren und addieren, erhalten wir: an 2 (2n—2)(2n—4) © 
2 9 i 
on [ean dex | (a8 + aay + ( a (x2 + 42)? /2 4. 
é 
(2n—]) (2n—3)---5-3 , : (2n—1) (2n—38)---5-3 ais 24 4 2 
(2n—2) (2n—4)--4-27 (at ot) sein oes 14.9% log (~ ned ee ®) 
Dann isi ¢ 
- =a 
— n— = n—| 3 
an [an [ten nl ere (Coat a PO Cc a ace 
(2n—1) (2n—2)--- at a\'hoon =e (2n—1) (2n—8)--:5-3  y2"+1 ate tr) 
AF m2) ee ic rales * (2n—2) (2n—4)---4-2 ° 2n41 | LEER ed 


Der zweite Summand in Gl. (7) ergibt sich aus Gl. (9) durch Vertauschung von x mit y. So erhalten wir: 


__ 1-3--+-(2n—1) 1 ontl ees Kat ona x + Vx2+ v2 
Pan 9-475, eral log ( )+y log £ + 


ar y(x? Hr y)'ts n—1, 2n—1 n—2 (2n 1)(2 n—3) 2 Ye 4 A 
<r 2n(2n+1) (x*?+-?) 2+ oF E 9 (x x?+-y?) yt on _9)(an—4) +y?) [x ay +: 
(2n—1) (2 n—3): : 3 2n—2 2n—2 —_ to) 
ee aiee 4-2 ty Ip @—1L2-) 
; | Pe fee 2 
und 9-1 (%,9) = axlog (74) ar +y log (ante ist) 
Mit den Abkiirzungen an 
u=ka, v=kb, w=k ya?+b2 (11) 
bekommen wir dann aus Gl. (3) und Gl. (4) fiir die Schalldruckkomponenten paz und pm im Eck punkt 
eines strahlenden Rechtecks mit den Seiten a und b die Formeln: (12) 


ama aes a |e Wer al] alee (ee Glee (lel 


15° 3.3 61) 7b7? 3-5 5-3 
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u+w Wee! ise oie +w Dp PRR SI al bats ph ; 
2pm —log ( . )\r-3 * ot legal to | tog ("2 hee ee Ret x (13) 
fl _ 13 17, 2 1-3-6 1 a iiytla. 108 Gat ot) 
pe SL re ai [w+ 3 25 66 7i(" HOE 3 ee ee 
1-3-5:7 1 2 2-4 2-4-6 ; 
Miia ie a see 
Aus der Reihenentwicklung der Besselschen Funk- . 
5 : : : v, ajlcos vb, b/eosp 
tion J,(x) ergibt sich: tye sae 
2 Pa+ipm= 1 | ave de + vik de, 


NOB He 
sae tr es ah 


Wed 


Pn at us) 


Damit kénnen wir die erste Zeile in Gl. (13) durch 


v u 
log ay Jo(x) dx + log (’ _ ")/ Jo(x) dx 
0 6 


ersetzen und fiir die Berechnung die fiir u=1, 2 ++: 24 


tabellarisch berechnete Funktion / Jo(x)dx be- 
nutzen [7]. Die fiir beliebiges u, v und w=/u?+v? 
absolut konvergenten Reihen (12) und (13) sind 
dann zur Berechnung von pa und p,, fir kleinere 
Werte von ka und kb geeignet. In der Tabelle I sind 
die entsprechenden Werte fiir ka=5, kb=5 zusam- 
mengestellt. Dabei waren in den Reihen (12) und 
(13) im Héchstfalle (ka=5, kb=5) zehn Glieder zu 
beriicksichtigen, um pa und pm» auf vier Stellen 
hinter dem Komma zu bestimmen. 


Tabelle I 
Pa 
ce 1 2 3 4 5 


+02749, +0,2588 | +0,2256 
+0,4483 | +0,4112 | +0,3423 
+0,4798 | +0,4061 | +0,3136 
+0,4061 | +0,3029 | +0,1956 
+0,3136 | +0,1956 | +0, 0901 


+0,1429 
+0,2405 


1 40,2405 
2 

3 | +0,2749 

4 

5 


+0,4012 
+0,4483 
+0,4112 
+0,3423 


+-0,2588 
+0,2256 
Pm 
a EOE: 2 3 4 5 


+0,1292 
+0,0119 


+0,2237 
+0,2253 


1 +0,2253 
2 

3 +0,1759 

4 

5 


+0,1903 
+0,0898 
+0,0177 
+0,0119 


+0,1759 
+.0,0898 | 
—0,0509 | —0,1375 | —0,1311 
—0,1375 | —0,2177 —0,1830 
—0,1311 | —0,1830 | —0,1356 


+0,1357 
+0,0177 


+0,1357 | 
+0,1292 


Fir groBere Werte ist eine andere Entwicklung 
(nach fallenden Potenzen von w) bequemer, die wir 
nun ableiten wollen. Fihren wir in Gl. (3) Polar- 
koordinaten ein: 

E=pcosy, n=esinp, dF—d&dyn=—ededy, 


so ergibt sich an Hand der Abb. 4 aus Gl. (3): 


(15) 


Pot iPn=5. [ote / go thaleost a gas | ctileostay| 


Wir setzen nach Gl. (11) ka=u. Dann ist: ’ 
(16 
7/2 


vy 
| Picea ats he fs —iu/eos dy we —tufeos Yay) 
0 


Wir fiihren eine neue Variable t== : on y ein und 


setzen : 
7/2 ; 
—iu/cosv Yeh? 
— 7 du a fe ae 
f(u) [e : tVe—1 
i 
Dann ist [8]: 
d : ie T iy 
a Saif rea dt = — 5 [Jo(u)—iYo(u)] 
i 


und daher 


= -f : [— [adorei [Ysera 


fu) (17) 


5 


Lk a 
Abb. 4. Zur Einfiihrung von Polarkoordinaten. 


Fir den zweiten Summanden auf der rechten Seite 


von Gl. (16) erhalten wir nach Einfiihrung von 
t=1/cos ¥ durch partielle Integration: 
7/2 me 
—iut 

—iu/eospa, _ |° ___ dg = 
| iat shiek | #1 
di 1/cos pa 
= @—iu/eos ys (5 —,)—iu e 14! are sin (1/t) dt. 


1/cos Ys 
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Ersetzen wir den arc sin (1/t) durch seine Potenzreihe und benutzen die neue Variable z= ut, so folgt 
wegen ka/cos ),=kya? + b?=w: 


7/2 Be ° 

he nis Seay me Sean Bl omen 
foretanmnM(5 a) ilo f ar gels [So ae). 
Y, 


Die Integrale auf der rechten Seite von Gl. (18) lassen sich durch den Integralsinus si(w) und den 
Integralkosinus ci(w) ausdriicken. Dabei benutzen wir die Beziehungen [9]: 


oo 


aa aes ay si(w) 4 008 20 w w? eeu (—1)" 71?" 1 
ntl (Qn)! yw |2nQn—1) 2n(2n—1)(2n—2)(2n—3) (2n)! + 
w 
__ sin w 1 ww eh. ce ie a a (19) 
wr |2n 2n(2n—1)(2Qn—2) (2n)! 
nee, hee xe ee A cos w iy oD w2 gee lasers of 
So nll) il an 7 Daa 1)ne2) (2n)! fe (20) 
sin w ‘ Ww pe f w? ; pz yo ee 
eer on (a1) ona 1) (2n—2) (Qn—3)* v (2n)! ; 


Fassen wir entsprechende Glieder zusammen, so folgt aus Gl. (17), (18), (19) und (14): 


bi 


[ costujeos.) dy = : ji— | 5sx)dx |—( —s) cos w + U(u, w) (21) 
Bye Ws 0 
u ¢ u 
mit U (uw, w)—— si (w) iL Ty(x) dx + (752 S| aa [ Jo(e) as| aa 
fy 6 
2!cosw 3! sinw ih gs 
=P [Lg fouls 
4!cosw 5! sin w WETS ee RIFT 
hohe we Ue aN sees [a ya usw. 
Ya u 
f sin (u/cos ydg=— 5 [ Y,(«) da — 5 rs ws) sin w— V(u, w) (23) 
‘ F 
mit V(u, w)=— ci(w [a (ae ‘ee w )( - | )ax|— 
3!cosw 2! snw 1 wu : 
a ae e's te [u— 5 at Jo(x) de + (24) 
0 
(oe ee) 1 ee 1-3 u® =| 58 wink aw. 
= u 


at — 2317 2-451 

Es 1aBt sich leicht einsehen, daB die Reihen in Gl. (22) und (24) fir beliebiges u und w (w>u) absolut 

konvergent sind. Betrachten wir namlich fiir hinreichend groBes n das n-te Glied in Gl. (22), so gilt fiir 
den ersten Faktor: 

aay 

(2n—2)! 7 +(2n—1)1 Eel 5 


wn 


268 H. STENZEL: DIE AKUSTISCHE STRAHLUNG DER KOLBENMEMBRAN 


und fiir den zweiten Faktor wegen Gl. (14): 


Nok Teer ba fl (Qn—3)_u2—1 
W'S etch Ocal ae ae eae On sane a 
1-3°° ‘(2n—1) y2n+1 4 ses 


9-4--On  (@n+lyl ~ Qn+1)l" 
Daher ist das n-te Glied in Gl. (22) absolut kleiner als 2w/2n(2n +1) und die Reihe konvergent. Analog 
folgt die Konvergenz fiir die Reihe in Gl. (24). 
Fiir die praktische Rechnung ist jedoch meist eine asymptotische Entwicklung einfacher. Diese 
erhalten wir, indem wir den Integralsinus si(w) und den Integralkosinus ci(w) durch die asymptotischen 
Reihen ersetzen [10]: 


; cosw sin w 2!cosw 3!sinw 4! cos w 5! sin w 25 
8i(w) = w at Nat) eS we we hm (26) 
: cosw sin w 3!cosw 2! sinw 5!cosw 4! sinw 26 
SOLO) == aig oe \ aoe or ltl oe sR aaa (26) 


Dann fallen in Gl. (22) und (24) die Glieder mit dem Faktor / J(x)dx fort und wir erhalten aus Gl. (15): 


2p Sears 4 aimas2 U*(u, w)—U*(v, w) — 5 (1— cos w), (27) 
ty 0 
u v 
. T Tomi Te 
2x Pm = — 9 [Xeax— 5 [Xoe)ax— V*(u,w)—V*(v,w) — 9 sin w (28) 
6 6 
i ! ! si 1 u3 
mit U* (2) = (ee a ee ee [ua gi| to (29) 
cosw sin w 3!cosw 2! sinw 1 us 
¥*(u 10) = ( ww )u—( wt ws ) fu fl (30) 
Die Berechnung der Integrale in Gl. (27) und (29) erfolgt fiir u<16 durch die Formeln [11]: 
[ Su(2)ax = ; u | Ju(u)H’y(u ) + J,(u) Hy(u )}, [st \da— 5 u | Yo(u )H’Q(u) + Y,(u) A a(u)| (31) (82) 
d 


mit Hilfe der tabellierten Funktionen H,(u) und H,(u) [12] und der Beziehung [13] H’)(u)=2/z—H,(u). 
Fir u>16 benutzen wir die Annaherungsformeln [14]: 


Hy Gaya (-, = =) + O(u-%), H, (u)=Y,(u) + = (1 + i) + O(u-4) (33) 


u 


und erhalten daraus (mit einem Fehler kleiner als 5-10—5): 


f 3s(2)ae—1— 2) (1—, 2) I (u) , [ aes Baty) + (1-34) ¥(w). {34) (35) 


Die danach berechneten Werte sind in der Tabelle II angegeben. 


Tabelle II 
Pa Pm 
a ye ea ere SS SE Ore wer ie ee 
ave. 8 Cy PE Te ee is\_— 10 16 20 
8  +0,3818 | 40,3192 | +0,3456 +0,3049 +0,0902  —0,0387 _ —0,0047 +. 0,1025 
10 40,3192 +0,2623 +0,2949 40,2811 10 —0,0387 | —0,1246  —0,1208  —0,0018 
16 | 40,3456  +40,2949 | +40,3072 \ +0,2854 16 —0,0047  —0,1208 —0,0736 | —0,0124 


20 + 0,3049 +0,2811 + 0,2854 +0,2782 20 + 0,1025 —0,0018 —0,0124 + 0,0952 


. Abb. 5. Zur Berechnung der Schalldruckamplitude fiir den 
} Aufpunkt auf der Symmetrieachse. 


1 
2. Die Darstellung der Schalldruckamplitude auf der 
Symmetrieachse 


Wir zerlegen die Fliche der strahlenden recht- 
eckigen Membran in vier kongruente Rechtecke 
(Abb. 5), von denen das eine aus den beiden Drei- 
ecken mit den Flachen F, und F’, zusammengesetzt 
ist. Aus Symmetriegriinden ist dann nach Gl. (2) 
der relative Schalldruck im Aufpunkt M gegeben 
| durch: 


age 
Fi+F, 
' Fuhren wir Polarkoordinaten p, » ein, so ist 
dF, = ededg =rdrdo 
und 
% V-h?+-a2/cos*@ 
eo ikr ok 
- -dF, = | do e  dr= 
re pr 6 
9, 
a Fe er as). 


6 
Analog folgt: 


Q, 
Wists 92( —ikh —ikVk?-+b2/cos?9 
——dF,=<\e a Be do}. 
Tr ik 
F; 0 


Unter Beriicksichtigung von 9, + 9.—7/2 ergibt sich 
dann aus Gl. (36): 


2. 
arr m Or nwaesy eer 
preees m2 fe any atcos9 gs 
0 


, ?: 


ms / eR P+ Bt /coste ge (37) 


T 
0 


H. STENZEL: DIE AKUSTISCHE STRAHLUNG DER KOLBENMEMBRAN 


269 


Dabei sind ¢, und 9, durch tg 9, = b/a und tg9,=a/b 
bestimmt. Die Integrale in Gl. (37) werten wir gra- 
phisch aus, indem wir die Kurven 


cos Y (kh)? +(ka)?/cos? 9, sin (kh)? +(ka)2/cos? 9, 
cos V (kh)? +(kb)2/cos? @ und sin y (kh)? +(kb)2/cos*@ 


in Abhangigkeit von » darstellen und graphisch 
integrieren. Dann konnen wir pa und pm berechnen 
und die Amplitude p=y/p2 + p?, als Funktion von 
kh darstellen. 


20 


16 


0 4 8 12 1) 20 24 28 


Abb. 6. Schalldruckamplitude auf der Symmetrieachse 
2Qa/k= 4; 2b/A= 4, 
———— 2a/A= 6; 2b/A= 6, 
p 2 
Pp =)r2 + Pm? 


TUE sites ee 
—ikh__ soa ald alah nes 
Tw 
0 


m — ® 


04 


0 8 16 4 32 40 48 56 
h/d 


Abb. 7. Schalldruckamplitude auf der Symmetrieachse 


des Quadrats 2a/A = 10, 
wereeeee- des flachengleichen Kreises; 


p=V pat Pm 
1/4 a 
f 2-5 4 —iky h?+a*/cos* 
pa ipg =o th 8] VPRO ay 
0 
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amplitude durch die Werte Zwei und Null gekenn- | 
zeichnet sind [15], liegen die Maximalwerte der qua- if 
dratischen (bzw. rechteckigen) Membran wesent- — 
lich unter Zwei und die Minimalwerte wesentlich ~ 
tiber Null [16] (Abb. 7). 


3. Die Kurven konstanten Schalldrucks in der XY-, 

XZ-und YZ-Ebene 
Wir stellen uns die Aufgabe, die Komponenten | 
des Schalldrucks fiir die drei Ebenen, die XY-, XZ- | 
0 4 8 12 16 20 8B und YZ-Ebene (Abb. 5), in der Umgebung der strah- — 
—— h/n lenden rechteckigen Membran zu berechnen, um © 
Abb. 8. Schalldruckamplitude auf der Symmetrieachse daraus die Kurven konstanten Schalldrucks dar- | 
2a/% = 4; 9, = 0,59, 2b/A = 6, >, = 0,98; zustellen. Nach den friheren Ausfiihrungen geniigt — 
p=Vp24 Ete es, den Aufpunkt senkrecht tiber einer Ecke des | 
1 ¢./. ae Rechtecks anzunehmen. Sind die Seiten des Recht- 
Pa + Pp = sale. 2 i ean ikVh?+ a/cos? do — ecks durch a und b und der Aufpunktsabstand durch 

m5 h gegeben, so folgt nach Gl. (2): 


Qo : 
2 [ —ikVh?+b2/cos? 9 eV patty? 
eae I) do. + ipm=- d ’ 
a Bs Bary ‘| s [ae “VhRe+22+y2 


In den Abb. 6—9 sind die Falle 2a/A=4, 2b/A—4; Fuhren wir die neuen Variabeln u=kx, v=ky ein, 
2a/A=6, 2b/A=6; 2a/A=10, 26/A=10; 2a/A=—4, _ g0 ist 
2b/A=6; 2a/A=6, 2b/A=10 und 2a/A=4, 2b/A=10 Vv kh) gaa 
dargestellt. Dak Pn = 5 | dp [: Ae = tute a 

Im Vergleich zu der kreisformigen Kolbenmem- e+ ue + 02 
bran ist der Verlauf der Schalldruckamplitude der k hy? sg 
quadratischen (bzw. rechteckigen) Kolbenmembran phi ne ree a 
auf der Symmetrieachse grundsitzlich verschieden. V (bh) + -u? + v2 
Wahrend bei der kreisformigen Membran die auf- Fur h=0, d.h. fiir die X Y-Ebene, bestehen nach 
tretenden Extremwerte der relativen Schalldruck- den Ausfiihrungen des ersten Abschnitts keine we- 
sentlichen Schwierigkeiten die Be- — 
rechnung durchzufiihren. Aber fiir ~ 
h+0 sind wir nicht in der Lage 
allgemeine Formeln anzugeben. Wir 
ziehen es deshalb vor, fiir alle drei 
Ebenen graphische Verfahren anzu- 
wenden und diese durch Berech- 
nung der Schalldruckkomponenten 
in der X Y-Ebene zu kontrollieren. 
Dabei sind die neuerdings ver- 
offentlichten Tafeln iiber die ver- 

—h/A allgemeinerten Sinus- und Cosinus- 

Abb. 9. Schalldruckamplitude auf der Symmetrieachse Integralfunktionen [17] von besonderem Nutzen. 
2a/A= 4, ©; = 0,38, 2b/A=10,9.=1,19; Fiihren wir die dort tabellierten Funktionen ein: 


— — ~~ 2a/A=6, 9,=0,54, 2b/A=10, O2= 1,03; 
pa ae z 
p=Vr2+ Pm? sin Yu? + a? 
8(¢,x)= | o-~ —=—de 
yu + az 


91 —— 
‘ pat 2 2 2 
o“ikh_® fg ik) h?+a?/cos pr ae 


(38) 


Pat 1Pm = 


z 
7) a os yz + a2 
2 vara bcos" 9 Ola.) 2 a eee 
open e: e jubtat 


80 folgt aus Gl. (38): 


kb 


1 = 
Pau= af 8 [V (kh)? + v2, ka] dv, 
0 


kb ka 
1 du 
Pn = 2 dv = cen ieee 
™ V (kh)? + u2 + v? 
0 0 


— J, | CW (bh) ke. (40) 


(39) 


0 


: Dabei lift sich das erste Doppelintegral in GI. (40) 
_allgemein auswerten. 

| Es ist: 

hb ha 


| du 
dv 


V(khy+ wet 
Keb 


ka-+y (kh)? + (ka)? + 0? 
ee x V (kh) 


0 


(a+ y(kh)* + (kea)?+ (6) 
= Hog | 1 (hy (by *)* 


kb-+-V/ (kh)? + (ka)? + (k ) 
+kal de 
nie ( V(keh)? + (ea)® 
ka-kb 
—kh arc ‘ oe pay 
khy (kh)? + (ka)? +- (kb)? 


Die zur 0 ag Sera Funktionen 


we kh)? ee 
dv = 
2+ u2 + v2 


“3 ‘ a S[// (kh)? + v2, ka} du, 


2 [% 


= s | C[y (kh)? + v?, ka] dv 


o(ka, kb, kh) = 


(41) 


1—cosy ( kh)? hale v8 


y(ka, kb, kh) = 
V(kh)2+u2 +02 


(42) 


wurden durch graphische Integration bestimmt 
und im Fall kh=0 explizit fiir ka,= 1, 2,3---10 
durch Kurven dargestellt (Abb. 10a—d und 11). 

Danach kénnen die Werte der Schalldruckampli- 
_tuden fiir die im Abstand 2/27 befindlichen Gitter- 
punkte der einzelnen Ebenen mit den entsprechen- 
den Werten der Schalldruckamplituden beziffert 
werden und durch Interpolation die Kurven kon- 
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] 12 2 
Abb. 10. o (kao, kb) = =— | du | du sin Yu* uae : 
ar Vu? + v? 
0 
(a) ftir kay = 0,53 1; 1,53°2; 255. 


(b) fiir kag = 3; 3,5; 4; 4,5; 5. 
(c) fir kay = 6,5; 6; 6,5; 7; 7,5. 
(d) fiir kag = 8; 8,5; 9; 9,5; 10. 


stanter Amplitude eingezeichnet werden. Die so be- 
rechneten Schallfelder der Quadrate mit 2ra)/A= 
2nby/A=3, 6 und 10 und der Rechtecke mit 27a,/A= 
10, 2nby/A=5; 2ray/A=10, 2nb)/A=2,5 und 2na,/A 
=6, 2xb)/A=3 sind in den Abb. 12—17 dargestellt. 
Dabei sind die Seiten durch 2a, und 26, bezeichnet!. 

2 Far die Ausfiihrung der graphischen Integration und 


die Darstellung der Kurven bin ich Herrn Dipl-Phys. 
BROSZE zu besonderem Dank verpflichtet. 
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56 
on 
40 ir 
—— 
32 + , = 
XY 
4 |e 
24 + 
YZ 


‘ WEE 


ol If 


Yo 


Zen 


0 1 2 g 4. 


INS 


kb ka 
1 y seme 2 2 
Abb. 11. (kay, kd) =a /™ | 1—cos/u* + v? 
0 0 


II. DIE BERECHNUNG UND GRAPHISCHE 
DARSTELLUNG DER KOMPONENTEN 
DER STRAHLUNGSIMPEDANZ 


Um die Strahlungsimpedanz der rechteckigen 
Kolbenmembran der Flache F zu finden, haben wir 
den gesamten auf die Fliche F ausgetibten Druck 
P* zu berechnen. Dabei schreiben wir P* in der 
Form: 


P*=Fe-owe”! (Pati Pm) (43) 
und bezeichnen die (dimensionslosen) GroBen P, 
und P,, als die (relativen) Komponenten der Strah- 
lungsimpedanz. Dann ist coPa gleich dem Strah- 
lungswiderstand. Zunichst betrachten wir den, auf 
das im Punkt x, y befindliche Flichenelement dx dy 
ausgeiibten (relativen) Druck “®(x, y)dxdy, der 
(Abb. 18) von den vier Rechtecken 1, 2, 3,4 mit 
den Seiten x, y; a—x, y; a—x, b—y; x, b—y herriihrt. 
Dann ist wegen Gl. (3) und (43): 


a b 
P.+ iP, =p de [ dy 0(x,9) 
und 0 0 
® (x, ¥y) = (2%, oe + 9(a—x, y¥) + 9(a—x, b—y) ate 
+9 (x, b—y). 


Vu? + 0? 


6 7 8 9g 10 
Fae de 


Nun ist aber: 
ab ab 
Ife (x, y)dxdy pee (a—x, y)dady= 
0 


ab 
cil, (a—x, b—y)dxdy= 
0 


ab 
=| {e (x, b—y) dx dy. 
Damit folgt: 


ab 
P,+iPy=p| fe(esdxdy. (44) 
6 6 


Ersetzen wir (x, y) durch die Reihe (4) so folgt: 


ks ks 


k 
P,+iP,= >; %—>-3] 92 + 5.5] Oh — 


AL k2 k4 ; 
+i|) 9-1-3919" 5 4p 
mit 
a b 


4 
= p/ [ols )aedy (n=—1,0,1...) (46) 


Fiir gerades n 14Bt sich die Integration in Gl. (46) 
ohne weiteres durchfiihren. Setzen wir 92,(%, y) aus 
Gl. (6) in (46) ein, so finden wir: 
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‘ouogy-7Z A top ut (q) 


oueqy-Z A top ut (q) "oueqa-Z_X rp ut (q) ‘ouoqa-A X sop ut (8) 
‘ouogn-AN top ut (e) SQ, = wy qrur ‘ouoqm-A NV rep ut (8) fg = MY grur 


UBIqUIOUTUEGIOY USYOSIZeIpeNnb Joep ppes[eyog “FL ‘aqqy UBIQqUIOUIUEGTOY, USYOsTyVIpeNnb sJop ppesTByos “Sl “GqV -[OS[ UOYosTZBApeNb sop ppos[eyos *ZL“AdV 


910 20 70 


©) & 
zi Zig 
OL 
0 
0 


o 


Wo) 


9b \7t\ZLIOF 
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Abb. 15a (oben). 
Schallfeld der rechteckigen 
Kolbenmembran in der 
XY-Ebene; 
kag 10, kbp 9. 


Abb. 15h (Jinks). 
Schallfeld der rechteckigen 
Kolbenmembran in der 
YZ-Ebene; 
kap=10, kbp =5. 


a 
Do, (a, b) = ate SE ot 


(7) gut—2 p2 
+ 3-(2n—1)-2-n 7 
nm a2n—4 pt 


+ §:@n—3)- 3-H T 


N\ p2n 
(n)8 


Dann kénnen wir das allgemeine Glied von P, in Gl. (45) in der Form schreiben: 


D>) 
j 1 
p=0 


1 pe" + 1 
. > ' D2,,(a, b)=— 


onab g2n+12n 
»(2n+1) (3 


8" @n2)! > 


und erhalten schlieBlich, wenn wir b/a=q setzen: 


n= 
n 
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(p + 1) (2p +1) (2n—2p+1) 


l — 2 QPar\2n+2 q XL 
P= DD ( x * (2n+2)1 “20 


Abb. 16a (oben). 
Schallfeld der rechteckigen 
Kolbenmembran in der 
XY-Ebene; 
kay=6, kby)=8. 


Abb. 16b (rechts). 
Schallfeld der rechteckigen 
Kolbenmembran in der 
YZ-Ebene; 
kag == 6,0 kg ——3- 


Abb. 17a (oben). 
Schallfeld der rechteckigen Kolben- 
membran in der X Y-Ebene; 
kag= 10, htbp= 2,5. 


Abb. 17b (rechts). 
Schallfeld der rechteckigen Kolben- 
membran in der YZ-Ebene; 
kag=10, kbp = 2,5. 


an” —2p p2P 


mit 


ay cae : gars Mae 
rot => Pp ) pF p+ NG apFN- 
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y Abb. 18. Tabelle III 
Zur Berechnung der Strahlungsimpedanz. 


0,01548 0,03093 0,07687 0,15056 
0,05706 0,11372 0,27751 0,51011 
0,11288 0,22391 | 0,52895 0,87305 
0,17034 0,33543  0,75399 1,07096 
0,22195 0,43260  0,90789 1,07698 
| 0,26729 0,51406 | 0,98900 | 0,98873 
0,31047 0,58767 | 1,02359 | 0,92521 
0,35556 0,66117 | 1,03904 | 0,93637 
0,40021 0,73580 1,04652 | 0,99124 
0,45160 0,80648 | 1,04218 | 1,03155 
0,53878 0,91499 | 0,98931 0,99694 


SOP PEEPLES 
ocoomronroanono an 


Mit Hilfe der fiir q=0,1; 0,2; 0,5; 1 berechneten 
Funktionen (s. die Tabelle am Schlu8 der Abhand- 
lung) ergab sich fiir P, die folgende Tabelle ITI. 

Um fir ungerades n die Integration in Gl. (46) 
_ auszufihren, verfahren wir ahnlich wie bei der Inte- 


gration fiir 9), nach Gl. (7), indem wir die Inte- 
gration tber das Rechteck durch die Integration 
iiber die beiden Dreiecke ersetzen. Dann erhalten 


wir: 
a b 
b g a ae 
4 4 
®,,_1(4, b) = F dy Pon ay (Bs y) dx + F dx Pan —1(%> y) dy. (49) 
OnE 0 0 
Unter Beriicksichtigung von Gl. (10) ergibt sich dann: 
| ae : 
(2n + 1) vous goat iL Pon—1(% Y) x= f(E y—FO, y) (50) 
6 
mit 2 2 E 21 nd 

© GHP log ZAVEAE) + ye tog (EAVES) 1 (51) 


+ (—1)"y(E2 + 2): [AGy?”— AM y 222 (24 y2) 4 Ay 2m —4 (ee 4 y2)2—_- + + (—1)" An? (2 + y*)"], 
F (0, ¥) = (—1)ry""* 7A) — AD + + (14? | 


und ee *) 
: a —1\r 
Qn+2 mi ) 
ines) 
(ny__ \n—l n—l —1 


n—1).4, = (i) rn 7) 22 (ac be. ae (a) et : (2n—4) 


1 2 2-4-+-(2n—2) 2-4---(2n—2) 
(n) _ meee te ue y 3. 2S" 
(Ont me PS 45-5. an 1) 1 35> Onl): 
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Setzen wir in Gl. (50) fiir — den Wert ay/b ein und integrieren dann von y=0 bis y=b, so erhalten wir 
zur Berechnung der ersten Summanden in Gl. (49), indem wir die Abkiirzung: 


Ans AM ae ee eas (ye. 


benutzen : 
2 an a2rt+2h b+ye+b? 
(2n+1)7- fare iter ie = ams aera) | es } ce 
b2nt+2q (ata? + B2 pens yar +b 
Pee | aaa (—1)" a4! +e 2n+3 


: [ 4) (n) b2n +2 AM h2n (q2-4-b2) + ADH2R—2 (a2 +b2)2—. - - +(—1)"A™6?(a? +b)" ’ 


Zur Berechnung des zweiten Summanden in Gl. (49) haben wir auf der rechten Seite von Gl. (52) a mit by . 
zu vertauschen. Durch Addition ergibt sich dann (mit der Abkiirzung q=b/a): 


4 3-5 ++-(2n—1) 


Pn—1(8,))= gy gata St ee) (53) 
mit a8 
‘1l+yl+@ (I + got") yea 
veon+1 (q) = aarp lon ( (q+V1+@?) ee = 5 log ( aa *#)_(1p “Graig 48 aa (54) 
(—1)"y1+ q? (n) 2n+2)\__ 4”) 2 2 
pa (2n+3)q AX (i a2") AL + o) (La ect 
+ ADL + g-2) (14g —- + (MAM tay. 
Aus es 
l (1h 
P= ZL a! (Qn)! — Don—1 (a, b) 
n=0 
folgt dann schlieBlich: 
aN an\2nt1 1 
Pa OT) ae tymtal Mo ") 
n=0 


Die Berechnung von t2n+1(q) ist wegen der auftretenden Summen: 


(my__ (™ 2 aan ‘) 2-4 as) oul ee (?) 56 
% (p) +3 p }*35\ p } 7° *3-5-- @m—2p +1) \p er 
besonders fiir groBe m und kleine p ziemlich umstindlich. Dabei kénnen wir zur Vereinfachung (oder zur 


Kontrolle) von der Beziehung Gebrauch machen?: 


24> a Be) 


Cp” + Cnn = 1:3.-(Qp +1)-1-3--(m—2p +1)’ (67) 
Zum Beispiel folgt aus Gl. (56): 
Cn =1L, OP =1+ 3 ttset" aay i Sia 
Cam + : Cm m+ 2 (m— ee 2)+-: oe = a 5 
wegen Gl. (57): 14 34 Fett gp ee Le. “ee ; 
m3 (m— +3. 35 2 ie +35. ey l= p91 3-1-3: eH m—3- 


2 Durch den SchluB von n auf n+1 beweist man Gl. (57) zunachst fiir p=0 und dann fiir p = 1, 2... wieder durch 
den Schlu8 von n auf n+1 mit Hilfe der Beziehung: nab = os + om): : 
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Mit Hilfe der fir g=0,1; 0,2; 0,5; 1 berechneten 
Funktionen $2n+1(q) (s. die Tabelle am SchluB der 
Abhandlung) ergab sich fiir Pm die folgende Ta- 
belle IV. 


Tabelle IV 
at/A | b/a=0,1 b/a=0,2 b/a=0,5 bla=1 
0,5 0,10963 0,17713 0,30950 0,43327 
1,0 0,20523 0,32470 0,53467 | 0,66108 
1,5 0,27815 0,42642 0,63115 | 0,61989 
2,0 | 0,33086 0,48411 0,60997 | 0,41268 
2,5 0,36858 0,51315 0,52085 0,20961 
3,0 0,40279 0,52967 0,41830 0,12290 
3,5 0,43389 0,54273 0,33354 0,14239 
4,0 0,46317 0,55136 0,26929 0,18912 
4,5 | 0,48794 0,54930 0,21499 0,19468 
5,0 | 0,50553 0,53402 0,16408 0,16942 


Fir gréBere Werte von az/A wurde P, durch 
graphische Integration ermittelt. 

Der Gesamtverlauf von Pa und Pm ist in der 
Abb. 19 bei verschiedenem Seitenverhaltnis (6/a=q) 
in Abhingigkeit von ka) dargestellt. Will man 


12 7 Tel 


10 Ml Seat 


—-_ kag 


Abb. 19. Die Komponenten der Strahlungsimpedanz der 
rechteckigen Kolbenmembran bei verschiedenem 
Seitenverhaltnis (b/a = q) und des Kreises (—— —). 


die eimzelnen Kurven bei verschiedenem q mitein- 
ander vergleichen, so hat man zu beachten, daB bei 
gleichem ka, die zugehérige Rechteckfliiche F von q 
abhangt. Es ist dann anschaulicher, anstatt ka,als ab- 


hingige Variable die GréBe = yF= ~ V¥q=kayVq 


zu wahlen. Die so fiir gq=1 und q=0,5 erhaltenen 
Kurven sind in der Abb. 20 dargestellt. Man wird 
daraus den Schlu8 ziehen, daB die Strahlungsim- 
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pedanz eines Rechtecks, dessen Seitenverhaltnis 


zwischen 1 und 1/, liegt, angenihert durch die 


Strahlungsimpedanz des flachengleichen Kreises 
ersetzt werden kann. Die in der Formel (47) an- 
gegebene Beziehung wurde bereits frither auf andere 
Weise durch Einfiihrung des Strahlungsfaktors S 
abgeleitet. Dabei besteht die Beziehung: 


(58) 


Indem wir an die weiteren friiher [18], [19] abge- 
leiteten Beziehungen ankntipfen, wollen wir eine 
neue Formel fiir S angeben, die zur Berechnung 
fir groBere ax/A= x und br/A=y geeignet erscheint. 


—~ kag VF 


Abb. 20. Die Komponenten der Strahlungsimpedanz der 
rechteckigen Kolbenmembran bei verschiedenem 
Seitenverhaltnis (b/a = q) und des Kreises (~——). 


Wir gehen aus von der Formel: 


1s 2 Jeo 25 (59) 
S(x, ¥)= 90(*)— 9 gy R1(*) 5-4" BLY Pa") — 
* 2 1 2a 
: sin* x sin z 
0 
und 


1 [2-4---(2n—2) 
E -+-(2n—1) 


Pn—-1 (x)— Qn (x) — 2x2 


—-———— 8,4 2x)| . (61) 
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Mit Hilfe der sich aus Gl. (60) und (61) ergebenden Gleichungen: 


] 1 
Po(%) = Pol(*) + 553 — Oya? 


1 1 2 , O'Sp(2x) 
mil*)= 90l®) + 998 aa2 Tt at 


1 1 2-4 O!S,(2x)  1!8,(2%) 


e(™)=Ol#) + 958 Ont TB! dat det 
1 1 24-6 018,(2x) , 1!8,(2x) , 2!8 (2x) 
a7) e*)chiocs oe ana ae dee rey 


erhalten wir aus Gl. (59): 


s(x,y)= | eol#) + gaa (1 


Ph et l 28 SE OES 
“gi? tea ery aaa a a 


1 
2 
1 2 . O!S8,(2x) 1:3 2° ,/0! 1!8, (2x) 
af sang ed Oly (aa Cs) + Acie 


Fiihren wir noch die Funktion ¢(y) ein: 


2y i 
1 J,(2y) Uy 228 1-342 
vy) = gy f Sela) de It ga gt (62), 
iy 
und beriicksichtigen die Beziehung: 
sin? y 23 y2 Dove 
erty 4! 6! ; 
so ergibt sich schlieBlich: 
1 sin? 
S(2,9)= (2190) + ga [90-5 | + Dm (63) 


mit 


m! Sm (2x) 142? 1:3 25 ++ (2m—1) Q2m+1 y2m 


Am =~ qant3 W)—L+ 3 Gm pa yt t(D ae -2m (2m +2)! ies) 


123 1-3 25 1-3+--(@n—1)  22n41 (2n$1) 2248 
! tens Gee tp ee ae tera L SiS Ae ee a) < Ra 
Wl +5: g og BI tI Tegan On £2)17 <5q -(2n +2) (Qn-+4)! 


y2nt2, 


Daraus folgt, daB fiir das (n+1)te 
Glied A, der Reihe in Gl. (63) die Tahalls wv 
Ungleichung gilt: 


| | Se (2) | (y2/e)041 z=y | 2 | a vill, Renew Ream 10 
| An | < G53 (n +2)! (n +1) (2n +3) o| 0,67240 0,35775 | 0,24866 0,188620,15186 
} : b  0,54539 0,24403 0,16626 0,12631 —0,09915 
und da /S,,(2«)| unter einer endlichen Go: 0,36672 | 0,08730  0,04134 0,02382 0,01506 
Schranke liegt, ist die Reihe An sinty/y? 0,20671, 0,03580 0,00217 0,01529 0,00296 
fiir beliebiges x und beliebiges y Y—sin?y/y? 0,33868 0,208230,16409 0,11101 0,09619 
konvergent. Setzen wir ferner die 1 (v5) 0,04233 , 0,00651| 0,00228, 0,00087,  0,00048 
fiir groBes x und y geltenden Be- ea) ¥y ) 
. . 2 | 
ziehungen [19]: ets a(t) 040905 0,09381 | 0,04362 0,02469 | -0,01554 
- . 
S(x,y) ~ Quay’ EAm |+0,01152 —0,00188 —0,00010 +0,00018 —0,00010 
I S| 042057, 0,09193 0,04352  0,02487  0,01544 
20(*)* or aa Pa| 1,0710 | 0,9364 | 0,9974 | 1,0135 | 0,0828 
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in Gl. (63) ein, so folgt, 
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daB Anhang 
ec 
>» Am fiir groBes x und y gegen Yon 4 1(9) 
m=0 
Null geht und daher vernachlassigt 
werden kann. 0,17643239 0,2867589 0,51066731 0,7433022 
Zur Erlauterung haben wir in 3 0,05077961 0,1049279 0,30178947 0,7821081 
: is 5 0,03361633 0,0690219 0,21015131 0,8023884 
der Tabelle V die Berechnung fir 7 0,02690234 0,0552788 0,1714899 0,8827744 
x=y=2, 4, 6, 8, 10 dargestellt. 9 0,02308814 0,0476225 0,1523822 1,0359522 
‘Wir méchten danach annehmen. 11 0,02055306 0,0425902 0,1419373 1,2886325 
i as prlisag 0,01871396 0,0389721 0,1361248 1,6899751 
et ap oe ? ? ? E) 
da fir x=y=6 (miteinem Fehler 5 0,01730222 0,0362188 0,1332202 2,3240858 
unter 1%) die Annaherung gilt: 17 0,01617499 0,0340397 0,1323949 3,3328102 
Ble oe 65 19 0,01524907 0,0322646 0,1332199 4,9552962 
PT) ‘ (65) 9 0,01447073 0,0307869 0,1354716 7,5978052 
£ ast sin? y 23 0,01380534 0,0295348 0,1390457 11,9554183 
~ G(%) HY) + ay | ¥(¥) — 7 25 0,01322580 0,0284595 0,1439143 19,2263123 
io sh 27 0,1501088 31,4909292 
(x= y=6). 29 0,1577145 52,390132 


Fir noch groBere Werte (x = y= 10) 
konnen wir dann in Gl. (65) ¥(y) 
durch 1/y und ¢(x) durch 


n/2%—1/2x?—sin 2x/4x3 ani let 2-10-1 5-10—1 1 
ee 2 | 3,3666666-10—2 | 6,9333333-10—2 | 2,0833333-10—1 | 6,66666667-10—} 
ersetzen und bekommen die Nihe- 4 | 2,0168666-10—2 | 4,1397333-10—2 | 1,2708333-10—1 | 5,66666667-10— 
rung 6 | 1,4420395-10—2 | 2,9682590-10—2 | 9,3377976-10—2 | 5,52380952-10—1 
Qxy 8 | 1,1231504-10—2 | 2,3218851-10—2 | 7,5750248-10—2 | 5,93650794-10—4 
A Oe) (66) 10 | 9,2034632-10—% | 1,9118309-10—2 | 6,5185195-10—2 | 6,89754690-10—1 
= : “ae 12 | 7,7999418-10—8 | 1,6285246-10—2 | 5,8340311-10—2 | 8,55649906-10—1 
jie 2 (5 2x sin’y 14 | 6,7709502-10—3 | 1,4211111-10—? | 5,3724430-10—* | 1,12229999 
4 mx \ 2x y? 16 | 5,984227 -10—3 | 1,262744 -10—2 | 5,0578530-10—2 | 1,54344367 
(x= y= 10) 18 | 5,36323 -10—3 | 1,137908 -10—2 | 4,8479184-10—2 | 2,20916704 
he Vea 20 | 4,86059 -10—% | 1,03701 -10—2 | 4,717550 -10—2 | 3,26987625 
(Eingegangen am 1. April 1952.) 22 | 44454 -10—3 | 9.5380 -10—3 | 4,65261 -10—2 | 4,9777802 
24 | 4,0967 -10—% | 8,8403 -10—3 | 4,63987 -10—2 | 7,7586618 
Schrifttum 26 | 3,7997 -10—3 | 8,247 10—3 | 4,6774 -10—2 | 12,33652 
[1] Lord RAYLEIGH, The theory of 28 | 3,544 -10—3 | 7,737 -10—3 | 4,761 10—2 | 19,95174 
sound evel L. Macmillan, London ©)29 -| 3.321 -10—3 | 7,29 -10—3 | 4,89 -10—2 | 32,7398 
1877. Pel) Sy -10—3 | 6,9 -10—3 | 5,06 *10—2 | 54,278 
[2] PACHNER, J., J. acoust. Soc. Amer. $2. 49,0 “10—* | 6,6 “10—* | 5,8 ore 91,431 
21 [1949], 617; 23 [1951], 185. 38 pp habtdisntice 
[3] RIEGGER, H., Wiss. Veréff. Siemens 40 4 58. 
3 [1924], 67. 49 796 
[4] STENZEL, H., Elektr. Nachr.-Te ch. 44 1392 


6 [1927], 247. 

[5] STENZEL, H., Ann. Phys. Lpz. 7 [1930], 953. 

[6] MOLLOY, C. T., J. acoust. Soc. Amer. 4 [1948], 395. 

[7] JAEGER, C., J. Math. Phys. 27, [1948], 211. 

[8] WATSON, G. N., A treatise on the theory of Bessel 
functions. University Press Cambridge 1922, S. 180. 

(9] KEIKITIRO TANI, Tables of si(x) and ci(y), Naval 
Experimental and Research Establishment, Meguro, 
Tokyo 1931, 8. 6 u. 7. 

{10] KEIKITIRO TANI, [9], S. 2. 

[11] WATSON, G. N., [8], S. 752. 

[12] WATSON, G. N., [8], S. 666—697. 

[13] WATSON, G. N., [8], S. 329. 

[14] WATSON, G. N., [8], S. 333. 

[15] BACKHAUS, H. und TRENDELENBURG, F., Z. techn. 
Phys. 7 [1926], 630. 

[16] ScHocH, A., Akust. Z. 6 [1941], 318. 

{17] The Annals of the Computation Laboratory of Har- 
vard University 18, [1949]. Harvard University Press 
Cambridge, Mass. 

[18] STENZEL, H., Ann. Phys. Lpz. 7 [1930], 955 und 956. 
eae hier die Abkiirzung S,, (22)=y mad, 13), (22) 
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Nachtrag bei der Korrektur 


Herr R. L. Prirronarp, New York, mit dem ich 
iiber die Ergebnisse der vorliegenden Arbeit korre- 
spondierte, machte mich freundlicherweise auf die 
beiden im Jahre 1951 erschienenen Arbeiten auf- 
merksam : 


1. Nomura, Y. and Arpa, Y., On the radiation 
impedance of a rectangular plate with an infi- 
nitely large fixed baffle. Sci. Rep. Res. Inst. T6- 
hoku Univ., Series B, 1, 2 [1951], 337—347. 


2. Cueraty, D. N., O Soprotivlenii Priamougol- 
noi Plastiny Kolebliushsheisia v Vyreze Ploskoi 
Stenki (Concerning the resistance of a rectangu- 
lar membrane vibrating in an opening in a plane 
wall). Prikladnaia Matematika i Mekhanika 
(Moskau) 15 [1951], 439444. 
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Ich bin Herrn PrircHarp be- 
sonders fiir die Ubersendung der 
beiden Photokopien dankbar. 


AuBerdem erschien 1952 eine 


kurze amerikanische Arbeit (Let- 


ters to the Editor): 


3. Swenson, G. W., jr. and Joun- 
son, W. E., Radiation impe- 
dance of a rigid square piston — 


in 


an infinite baffle. J. acoust. 


Soc. Amer. 24 [1952], 84. 


Ie 


Die japanische Arbeit be- 
nutzt die Integration in der 
komplexen Ebene 
duensatz). Der fiir Pa ab- 
geleitete Ausdruck ergibt 
sich dabei in derselben Ge- 
stalt wie beiuns, wahrend fur 
P» ein wesentlich umstiénd- 
licherer Ausdruck erscheint. 
Doch 148t sich auch hier 
die Ubereinstimmung mit 
unserem P,, durch Reihen- 
entwicklung nachweisen. 


.Die russische Arbeit 


bringt das Ergebnis in Form 
eines bestimmten Integra- 
les. Daraus lassen sich un- 
sere Funktionen yo,(q) und 
bon41(q) mit Hilfe des Tay- 
lorschen Satzes berechnen. 


.In der amerikanischen 


Arbeit treten unendliche 
Reihen auf, deren Koeffi- 
zienten bestimmte Integrale 
sind. Diese Integrale lassen 
sich ohne Schwierigkeiten 
auswerten und daraus die 


Xon(q) und ve,41(q) bestim- 
men. 


0,1000000000 


0,0336666667 | 


0,0201686667 
0,0144203952 
0,0112315041 
0,0092034632 
0,0077999418 
0,006770950 
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0,200000000 | 
0,069333333 | 
0,041397333 | 
0,029682590 | 
0,023218851 | 
0,019118309—— 
0,016285246 | 
0,014211111 | 


| 


q= 5 


0,5000000000 
0,2083333333 
0,1270833333 
0,0933779762 | 
0,0757502480 | 
0,0651851946 | 
0,058340311 
0,053724430 
0,050578530 
0,04847918 
0,04717550 
0,0465261 
0,0463987 
0,0467740 
0,0476097 
0,0488957 
0,0506058 
0,0528560 
0,0555823 
0,0588617 
0,0627522 
0,0673255 


q = 0,5 


0,5106651 
0,3017888 | 
0,2101510 | 
0,1714894 | 
0,1523820 
0,1419371 
0,1361245 
0,1332199 
0,1323945 
0,1332194 
0,1354709 
0,1390448 
0,1439135 
0,1501068 
0,1577002 
0,1668104 
0,1775951 
0,1902534 
0,2050312 
0,2222241 
0,2421871 


0,005984227 0,01262744 
0,00536323 | 0,01137908 
0,00486059 ~—«*0,0103701 
0,0044454 | 0,0095380 
0,0040967 | ~=—-0,0088403 
(Resi- 0,0037997 0,0082471 
0,0035437 0,0077367 
0,0033208 0,0072933 
32 0,0031249 0,0069045 
34 0,0029514 0,0065612 
36 | 0,0027967 0,0063047 
38 | 0,0026579 0,0059825 
40 0,0025327 | 0,0057367 
42 0,0024191 | 0,0055147 
i 
Y3n+1(9) 
2n-+1| q= 0.1 | q= 0,2 
1 0,1764324 |  0,2867576 
3 0,0507797 | 0,1049273 
5 0,0336164 | 0,0690215 
7 0,0269024 |  0,0552785 
9 0,0230881 0,0476223 
11 | 0,0205530 0,0425899 
13 | 0,0187138 0,0389719 
15 0,0173022 0,0362186 
17 0,0161754 0,0340396 
19 | 0,0152489 0,0322646 
21 | 0,0144705 0,0307869 
23 0,0138046 0,0295349 
25 | 0,0132268 0,0284598 
27 | 0,0127193 0,0275261 
29 | 0,0122692 0,0267073 
31 0,0118665 0,0259838 
33 | 0,0115034 0,0253401 
35 | 0,0111741 0,0247642 | 
37 0,0108736 0,0242466 | 
39 0,0105980 | 0,0237795 | 
41 0,0103444 | 0,0233564 | 


1,0000000000 
0,6666666667 
0,5666666667 
0,5523809523 
0,5936507936 
0,6897546897 
0,8556499056 
1,122299922 


1,54344367 

2,20916704 

3,26987625 
4,9777802 
7,7586618 
12,3365194 
19,9511444 
32,7397988 
54,2782568 
91,4309050 
155,1518126 
265,5850180 
458,1855533 
796,0459560 


q=1 


0,7433024 
0,7821082 
0,8023885 
0,8827745 
1,0359517 
1,2886330 
1,6899748 
2,3241120 
3,3328105 
4,9552959 
7,5978080 
11,9554080 
19,2263122 
31,4909851 
52,3859481 
88,3071373 
150,5683845 
259,2891737 
450,4240962 
788,5187210 
1389,9353787 


Bemerkenswerterweise stimmen die so aus 2 bzw. 3 berechneten Funktionen 73,(q) und 43,,,;(q) yoll- 
stiindig iiberein, weichen aber in der Gestalt von unserem 7p,(q) und 4o,+(q) nicht unwesentlich ab. 
Und zwar ergibt sich: 


XEn(4) 


(2n +2) ¥S,44(q) = log (q+ VIF @) +q"42 log ( 


2q_|(1+q°)" 


ey ey ee | [21+ ( 


2n(1+q?)?— 


2n+1) (2n—1) 


2[L+qets—(1+ qr? 


(2n+3) (2n+4)q — ’ 


eae 
q 


2n (2n—2)--+2 


2n+1 


Ltt] +++ a1) Qn—l)--- 


3 cl+ge| + 


2.4..-2n 


2n+3  3-5---(2n+1) 


1 = 2 
a [1+9q2e+8 — (1-+-4?)"t+*:]+qy1 He fie *] + 3 (1+-q*) (0-4 q"*] ee 


2.4--+ (2n—2) 


+ 3.5..-(2n—1) 


(1 + e102} 


; 
| 
| 


TH. VOGEL: CONDITIONS AUX LIMITES DE L’EQUATION DU SON 281 


Wir haben diese Formeln zur 


Tabelle IV* 
Kontrollrechnung benutzt und 


ar/r | ‘a = 0,1 | la = 0,2 | b/a = 0,5 a=1 

die genaueren Werte in den Ta- aig fi | si i 
bellen fiir y3,(q) und 43, ,1(q) zu- 0,5 | 0,10969 0,17712 0,30950 0,43327 
sammengestellt. Mit Hilfe dieser 1,0 | 0,20510 0,32471 0,53467 0,66133 
2m , 1,5 |  0,27836 0,42639 | 0,63116 0,61996 
Tabellen wurden einige Unrich- 2,0 0,33054 0,48416 0,60997 0,41245 
tigkeiten der Tabelle IV fir P,, 2,5 0,36906 | 0,51304 0,52088 0,21042 
beseitigt (vgl. Tabelle IV*). 3,0 | 0,40205 0,52965 0,41831 0,12322 
3,5 0,43374 0,54276 | 0,33355 0,14265 
4,0 | 0,46343 0,55129 0,26947 0,18837 
4,5 | 0,48799 0,54945 0,21506 0,19302 
5,0 0,50541 0,53385 0,16415 0,14982 

| 


SUR LES CONDITIONS AUX LIMITES DE L’EQUATION DU SON 


Par THEODORE VOGEL 
Centre National de la Recherche Scientifique, C. R. 8. I. M., Marseille 


Sommaire 


Aprés avoir indiqué ce qu’une théorie ondulatoire de l’acoustique des salles doit conserver 
des résultats de SABINE, l’auteur discute la notion d’impédance normale spécifique, montre par 
quelle suite de raisonnements elle s’introduit, et les conséquences qu’on peut en tirer dans un cas 
simple de propagation. Les expériences faites dans ce cas ne semblent pas vérifier ’hypothése, 
laquelle demande donc a étre utilisée avec précaution. On termine en suggérant des modifications 
possibles & la condition aux limites, aboutissant aux notions d’impédances spécifiques normale 
et tangentielle. 


Zusammenfassung 


Zunichst wird darauf hingewiesen, daB jede Wellentheorie der Raumakustik mit den SABINE- 
schen MefBergebnissen vereinbar sein muB. Es wird gezeigt, durch welchen Gedankengang sich der 
Begriff des normalen spezifischen Wandwiderstandes einfiihren 1a8t und welche Folgerungen sich 
daraus fiir einen einfachen Fall der Wellenausbreitung ergeben. Die Experimente scheinen in 
diesem Fall jedoch nicht mit der Hypothese tibereinzustimmen, so daB bei Anwendung der Theorie 
Vorsicht geboten ist. Es wird daher vorgeschlagen, die Theorie dadurch zu erweitern, daB man 
neben dem normalen einen tangentialen Widerstand einfiihrt. 


Summary 


A wave theory of room acoustics should account for those of SABINE’s results which are sub- 
stantiated by facts. These are discussed, and it is shown how the boundary condition conducive 
to the notion of specific normal impedance is most naturally introduced. However, the conse- 
quences of that hypothesis do not seem in good accord with experiments, in one simple case at 
least; caution should therefore be used when applying this boundary condition. A more involved 
condition is tentatively suggested, taking into account two specific impedances: one normal and 
one tangential. 


1. Introduction: 
Remarques sur la théorie de W. C. Sabine 


On se propose de discuter dans ce qui suit quel- 
ques aspects du probléme de la propagation du son 
dans un espace clos; et il sera sans doute bon de 
commencer par définir briévement le point de vue 
ou l’on se place. On croit que le réle du physicien 
consiste 4 choisir ses fondations, qui seront des faits 
expérimentaux avérés; 4 construire sur ces fonda- 

_ tions un édifice rationnel, en suivant des régles com- 
munes sinon a toute la physique théorique, du moins 


a tous les édifices voisins; et 4 aboutir 4 la prévision 
de faits observables, qui seront effectivement avérés 
4 leur tour. Cohérence théorique et possibilité d’ap- 
plications expérimentales sont done deux condi- 
tions 4 remplir simultanément pour faire oeuvre de 
physicien. 

Si l’on examine de ce point de vue l’héritage lé- 
gué, en matiére d’acoustique des salles, par W. C. 
SasrnE [1] et ses continuateurs, on trouve, d’une 
part, un ensemble de faits avérés, pouvant se ré- 
sumer ainsi: importance primordiale, pour la qua- 
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lité de l’audition, de la facon dont le son s’éteint 
dans la salle aprés arrét du fonctionnement des 
sources; possibilité d’exprimer approximativement 
cette loi de décroissance par un nombre, dit «temps 
de réverbération», lequel devra, pour chaque type 
de salle, étre compris dans un intervalle optimum. 
L’autre part de l’héritage consiste en des énoncés 
tendant a expliquer les constatations précédentes, 
c’est-a-dire & les rendre vraisemblables a priori. Ces 
énoncés conduisent 4 une formule qui exprime le 
temps de réverbération en fonction d’une grandeur 
d’ailleurs mal définie, le «coefficient d’absorption» 
des parois. Enfin, on postule que pour chaque élé- 
ment de paroi, ce coefficient est spécifique du maté- 
riau constituant. 

Cette loi de Sapinr semble étre, depuis assez 
longtemps déja, l’objet d’un malentendu entre 
physiciens et techniciens: ces derniers croient qu’il 
en existe une justification théorique cohérente (celle 
de JAEGER ou ses variantes), dont les premiers con_ 
naissent toute la faiblesse; mais par contre, on voit 
des théoriciens essayer de mieux justifier la loi, ou 
d’en raffiner ]’6noncé, dans la croyance que l’en- 
semble des mesures de coefficients d’absorption et 
des conclusions qu’on en tire donne entiére satis- 
faction aux techniciens des salles d’audition. Or, 
voici ce qu’écrit l’un des plus compétents de ceux-ci: 

«Il ne semble pas douteux que les calculs et les 
recommendations de SaBinz ne conduisent a éviter 
la plupart des défauts du premier ordre que peut 
présenter l’acoustique d’une salle; mais l’expéri- 
ence semble bien montrer que des défauts du second 
ordre, encore importants, ne puissent étre corrigés, 
ni par une application plus méticuleuse, ni par des 
perfectionnements de la théorie de Sapinz. Ce que 
nous disons 14 ne peut guére étre mis sous forme 
mathématique précise; mais c’est une conviction 
que l’on acquiert au contact quotidien des pro- 
blémes réels» (J. Morr) [2]. 

Il semble donc que le mieux soit d’abandonner 
carrément toute référence 4 la théorie géométrique 
de SAaBINE, en ne conservant que la partie avérée 
de son héritage; il s’agit donc de construire une 
théorie de la propagation du son dans les espaces 
clos, qui soit cohérente avec le reste de la mécanique 
des fluides (puisqu’on a affaire au mouvement de 
lair sous l’action de forces), et qui arrive 4 rendre 
compte de l’importance de la réverbération, et a la 
préciser. Le probléme ainsi posé est évidemment un 
probléme aux limites, c’est-a-dire qu’il faut adopter 
une loi du mouvement pour les molécules d’air, 
exprimée par une équation aux dérivées partielles, 
et une condition aux limites, que devra satisfaire 
la solution en tout point des parois: cette condition 
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exprimera de son cété une propriété physique in 
trinséque des parois. 

La question a été traitée sous cet angle, & la suite 
de Lord RaviEreH qui lui consacra quelques pages ~ | 
rapides mais fondamentales, par F. H. VAN DEN | 
DuneeEn d’une part [3], par Morsz, Borr et leurs 
nombreux collaborateurs de l’autre [4]. Des résul- 
tats théoriques de-grande valeur ont été obtenus, 
et_le trésor-des faits expérimentaux a été enrichi 
dé nombreuses mesures; mais une mesure n’a de 
signification que par rapport 4 une vue théorique 
des phénoménes, et celles dont il s’agit reposent 
sur la notion d’impédance acoustique, suggérée par 
les travaux théoriques cités: et peut-étre l’ardeur 
des expérimentateurs les a-t-elle emportés fort en 
avant sur un terrain qui n’était pas entiérement 
assuré. Telles quelles, leurs mesures fournissent des 
données dont il n’est plus permis de ne pas tenir 
compte; mais on croit utile de jeter cependant un 
regard en arriére, pour essayer de mieux se rendre 
compte de la validité des hypothéses adoptées aussi- 
tot qu’avancées. 


2. L’équation de d’Alembert 


et ses conditions aux limites 


Z 


En ce qui concerne l’équation différentielle du 
mouvement, il faut d’abord se demander s’il est 
légitime de la linéariser en supposant les mouve- 
ments infiniment petits. Or cette petitesse doit 
s’entendre par rapport aux plus petites dimensions 
intervenant dans la description du domaine d’inté- 
gration; les amplitudes de vibration varieront prati- 
quement entre 5-10—* cm pour un son 4 peine au- 
dible (2:10-4 baryes de pression maximum) 4 la 
fréquence la plus basse, et 0,06 cm pour un son trés 
intense (2: 10— baryes, ou 140 db) & une centaine 
de Hz: de sorte que des détails de l’ordre du cen- 
tiéme de micron dans le premier cas, de l’ordre du 
cm dans le second, seront «grands». L’hypothése 
de la petitesse des mouvements sera done toujours 
acceptable. 

Peut-on dés lors utiliser l’équation classique de 
p’ALEMBERT, éventuellement complétée par un 
terme dissipatif linéaire en 39/dt? On observera 
simplement, 4 cet égard, que méme écrite en termes 
de pression, l’équation de pD’ALEMBERT suppose 
lexistence d’un potentiel des vitesses, et que cette 
existence n’est assurée que lorsque le théoréme de 
LAGRANGE s’applique, c’est-d-dire lorsque les accé- 
lérations restent constamment bornées: ce qui ex- 
clut les problémes 4 ondes de choc, coups de cla- 
quettes, etc. Une expérimentation au moyen de tels 
dispositifs ne pourra done étre légitimement inter- 
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prétée par référence, fit-elle implicite, 4 une théorie 
fondée sur l’équation de p’ALEMBERT. 

Le choix de la condition aux limites souléve des 
questions beaucoup plus délicates. 

I] faut d’abord remarquer que l’on serait 4 peu 
prés désarmé, tant du point de vue expérimental 
que du point de vue théorique, si le systéme diffé- 
rentiel n’admettait pas de solutions séparables quant 
au temps, c’est-a-dire de la forme p= ©®(x,, x,, x3): T(t), 
des modes et des pulsations propres; et ensuite, que 
les difficultés théoriques resteront considérables si 
ces modes propres ne forment pas une base ortho- 
gonale (au moins en quelque sens élargi), permet- 
tant de développer en série la solution de tout pro- 
bléme non homogéne. 

En ce qui concerne le premier point, on démontre 
 aisément le lemme suivant: La condition néces- 
saire et suffisante pour qu’un systéme dif- 
férentiel composé de l’équation 


1 2 
A.9 — ge 19 hae =0 


et d’une condition aux limites soit sépa- 
rable par rapport au temps est que la con- 
dition aux limites soit homogéne par rap- 
portagetases dérivées partielles par rap- 
port au temps. Si l’on pose 


eae 
A,® met atu ine 
— = ¢ = ’ 
> . onst =~ 
Je facteur temporel devient T = exp t, avec 4 
racine de |’équation 


R24 he2A—22—0. 


Aprés division par une puissance convenable de 
exp dt, on trouvera bien un systéme différentiel 
purement spatial en ©; mais il convient d’observer, 
d’abord, que l’existence ni l’unicité de la solution 
ne sont assurées en général (lorsque la condition 
aux limites n’est pas linéaire); et que lorsque le 
probléme est linéaire, c’est-d-dire que la condition 
aux limites est de la forme 


Lanny 8x? =0, 


le systéme en ® dépend du paramétre X% aussi bien 
par son équation différentielle que par sa condition 
aux limites 

Lann "92D = 0 


(sauf cas exceptionnels), 4 étant une fonction de >. 
Cette particularité, qui explique certains énoncés 
_incorrects que l’on trouve dans les travaux récents 
(comme quoi les différents modes propres corres- 
pondraient 4 des conditions aux limites différentes 
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— ce qui est contradictoire avec la notion méme 
de mode propre), a une incidence importante sur 
Vorthogonalité des modes: en effet, si ® et ©’ en 
sont deux, correspondant aux valeurs propres di- 
stinctes 1 et 2’, le théoréme de GruEn s’écrit ici 


2/2 — 22 key a0’ 
ais ea ee 
c on on 
D s 


Il y aura orthogonalité stricto sensu dans D si la 
condition aux limites est de la forme 


+ «D=0 (« constante absolue), 
et dans ce cas seulement; c’est un résultat bien 
connu; mais il est bon d’observer (ce qui ne semble 
pas avoir été signalé jusqu’ici) que le probléme est 
encore maniable si, dans la condition aux limites 
ci-dessus, « est de la forme u— «,(x) A? (% constante 
absolue, «,(x) éventuellement fonction des coordon- 
nées): dans ce cas, en effet, la condition 
Re 2 

ane [Peace ay ahaa) ds] = 0 
exprime une orthogonalité des modes propres dans 
D +S, relativement 4 une distribution de STIELT- 
ges de densité du=dr dans D et du=a,(x)dS 
sur S. 

On voit par ce qui précéde combien les calculs 
seraient facilités si l’on pouvait adopter une hypo- 
thése qui conduirait 4 la condition aux limites 


dD/In+aD=—0. Or cette hypothése s’accorderait 
bien avec les idées que l’on peut avoir sur les deux 


cas extrémes: frontiére parfaitement rigide (aucune 
particule ne peut la traverser, done condition ex- 
primant que la composante de la vitesse normale a 
la frontiére est nulle, d9/In=0 ou I@/dn=—0), et 
surface fluide libre (pression constante, soit nulle, 
sur la surface en question, done condition pe: 39/3t=0 
ou encore ®=0); de plus, elle jouirait de tout le 
prestige que lui confére l’analogie thermocinétique, 
puisque 
sé + x0 =0 

on 
est la condition aux limites classique pour le cas 
du passage de chaleur partiellement empéché. 

Pour formuler ’hypothése en question, il faut 
lexprimer en termes d’observables, vitesse des 
particules et pression: or, en multipliant les deux 
membres par le facteur temporel exp it, on voit que 
exp at- g®/an =d9/an est la composante normale 
de la vitesse des particules au voisinage des parois, 
et ®-exp 2 est proportionnel a la pression sur la 
paroi, p-39/8 = prA®-exp 2; de sorte que l’hypo- 
thése s’énonce: le rapport de la pression de 
lair sur la paroi a la vitesse normale est 
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indépendant de la répartition spatiale de 
état vibratoire. Autrement dit, ce rapport est, 
en tout point de la paroi, spécifique de celle-ci. Ce 
rapport a d’ailleurs les dimensions d’une impédance, 
et on l’appellera impédance spécifique nor- 
male de la paroi. 


3. Discussion de la notion d’impédance 
spécifique normale 


Il faut maintenant se poser deux questions: 
1. Phypothése faite permet-elle d’utiliser les faits 
relatifs A la réverbération, et comment? 2. les con- 
séquences qu’on peut en déduire sont-elles toutes 
vérifiées par l’expérience ? 

Pour ce qui est du premier point, on voit immé- 
diatement que le temps de réverbération d’un mode 
propre est lié au décrément, c’est-a-dire a la partie 
réelle de la valeur propre qui correspond au mode. 
Les considérations qui suivent s’inspirent, avec 
quelques modifications, du beau travail de VAN 
DEN DUNGEN. 

Soit ® un mode propre correspondant a la valeur 
propre ; en désignant par des astérisques les gran- 
deurs imaginaires conjuguées, on aura les systémes 
différentiels 


2 2\* 
aso - ©=0(D) |a.0*— a o*—0 (D) 
et 
aD pA eae isp pa* 
Pe Pao vigy ae on)“ ®) 


et par application du théoréme de GREEN aux fonc- 
tions ® et O*, 


* 
[ (@s,0*—0*A,0) ds— [(e5, ot oe) do 
D s 


on 0 


2 ¥*__ 32 r* 
soit CO [eota— ef (e- 
D NS) 


On peut d’ailleurs supposer normés les modes pro- 

pres, qui ne sont définis qu’é un facteur d’amplitude 

prés; de sorte qu’il est licite de poser /®®*dr=1 
D 


*) OO* do. 
3 


sans rien perdre en généralité. Si l’on met en évi- 
dence les parties réelles et imaginaires de A et de 3, 
soit 


A=S+ie; 3=—r-+is, 
on aura done 
4iwSd 
oe =3 | (54—,) o*do— 
pe 3 3 
$ 


me O0* do= 


sOO* do [ee 
rts? r2+ 52 


= 2i 


i Aca 
i 
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d’ot pe 
nate et se 
5 Seale at: 
ra, 2 sDO*dc 
oc? f r?+s2, 
iS) 


Les deux cas particuliers les plus simples sont | 
ceux ot 3 est le méme en tout point de la frontiére, | 
et celui ot de plus sa partie imaginaire est nulle; 
on a alors respectivement 
its ce rel 
Qeo(r? +s) + 
en désignant par I Vintégrale [D®*do, qui peut 

8 ’ 


2] 
et 


rg + pe*s I 2 


d’ailleurs recevoir la signification physique de mo- — 
yenne des modules de pression sur la paroi au 
temps 0, 4 un facteur de proportionnalité prés. | 

On voit que méme dans le cas trés simple ou . 
V’impédance normale est = constante réelle, sa — 
connaissance ne suffit pas & déterminer 5; il faut 
encore connaitre la répartition moyenne des mo-— 
dules de pression sur la paroi, c’est-d-dire avoir 
résolu préalablement le probléme de l’état vibra- 
toire de lenceinte. 

Mais peut-étre a-t-on été trop difficile, et trop 
précis; peut-étre faut-il donner & la notion de temps 
de réverbération un contenu & la fois plus large et 
plus vague ? I] est connu que l’équation des ondes 
admet des solutions fondamentales planes, de sorte 
qu’il est loisible de ne considérer dans tous les pro- 
blémes que des superpositions de telles composantes, 
d’incidences diverses. Admettons que le temps de 
réverbération moyen qui caractérise la qualité 
dune enceinte pour une pulsation donnée soit, 
comme le veut SaBINE, inversement proportionnel 
4 un coefficient d’absorption moyen, ce dernier s’en- 
tendant pris comme la moyenne des coefficients 
d’absorption aux incidences diverses, et celles-ci 
étant uniformément distribuées: si l’on suppose la 
paroi plane et indéfinie, ’hypothése d’une impé- 
dance normale spécifique permet facilement le cal- 
cul et la mesure de ce coefficient sous l’incidence *: 
si x0 est l’équation de la paroi, on aura pour 
systéme différentiel 

20 dD w? a®D 

ox? T By2 Ta Rm (355 


d’ot la solution 


a 0, 


®= exp(ioy sin $/c) [exp (iwx cos $/¢) + 
+r-exp(—iw~x cos $/e)] 
avec 
(1—r)3 cos } — (l+r)pe=0 
soit 
3h cos }— pe 
~ 3 cos 9+ pe” 


Posons 
3—=Zee-expia, Z et « réels; 
il vient 
rré— R2 — 1 +2? cos? 3—2Z cos « cos > 


1+ Z2 cos?9-+ 2Z cos « cos 9° 


Le coefficient d’absorption moyen serait alors 
T/2 
a=2/a(%) cos? sind d9;  a(9)=1—R?. 
0 
On remarquera que de la mesure de a(%) sous deux 
incidences distinctes, d’ailleurs quelconques, on 
peut déduire Pimpédance normale. Dans le cas 
théorique de Vabsorbant type «gdteau de miel» 


étudié par Lord RayLeicu, «=7/4; on peut poser 


Zy2=1/M (0<MS1) 


et alors 


LI 2M? 

Cette expression passe par un maximum d’environ 
0,83 pour M 0,4, et @ est alors sensiblement égal 
a (55°). C’est pour cette raison que, désirant mesurer 
le coefficient d’absorption sous incidence donnée 
pour avoir des résultats ayant une signification 
physique précise, nous avons adopté depuis 1945 
Vincidence de 55° au laboratoire d’acoustique du 
Centre de Recherches Scientifiques, Industrielles 
et Maritimes de Marseille. En réalité, il faut bien 
dire que le choix de cette incidence n’a pas une im- 
‘portance trés grande, du fait que les mesures mon- 
trent que le cas de RayLercu «=7/4 est assez rare- 
ment rencontré dans l'étude des matériaux in- 
sonores courants. Toutefois, la mesure basée sur le 
principe ci-dessus présente un autre intérét, qui 
importe davantage ici: c’est de permettre une ex- 
périence cruciale vérifiant la validité de Vhypothése 
d’une impédance normale spécifique: en effet, lors- 
qu’on fait tomber, en chambre sans écho de grandes 
dimensions, une onde plane sur un assez grand pan- 
neau plan incliné & 9, il se forme un systeme d’ondes 
stationnaires, et l’on peut aisément observer les 
modules maxima et minima de pression au voisinage 
du panneau; c’est en somme I’expérience du tuyau 
acoustique de E. T. Paris, faite dans de meilleures 
conditions (pas d’effets de paroi, peu d’effets de 
bord). 

Que |’on forme alors l’expression 


oo |Pmas| [Pmin| : +R tz Se 1 +R? 
3 [Pmin| (0 jaa a LR ya Re? 
elle vaut 
1+ Z? cos? 9 


i Z cos « cos 9’ 


et si Z et « sont vraiment indépendants de la distri- 
bution ®, done de 9, on peut les éliminer entre les 
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Or 


expressions de m correspondant & trois incidences 
distinctes; pour préciser, il doit exister entre les 
trois valeurs observées de m la relation 


mM, COS 3, — My COS 4, 
Mo COS % — m, cos a 


cos? $,— cos? 9, 
COS? 9y— cos? 9, 


En particulier, si les angles choisis sont 0°, 45°, 60°, 
cette relation prend la forme extrémement simple 


Ol, = 


w= Onn 


De nombreuses vérifications expérimentales ont 
été faites au Laboratoire du C. R. 8.1. M. A titre 
d’exemple, le tableau ci-aprés reproduit les résul- 
tats obtenus avec un échantillon de Tropal: 


Tableau I 


fréquence [Hz] 


S NWN] 


58 | 0,58 | (valeur théorique: 1) 


Les mesures ont été recommencées bien des fois, 
avec une dispersion relativement faible, et il est 
peu probable que Verreur commise sur chacune 
delles puisse atteindre 10°%, de sorte que lerreur 
sur uw. n’est sans doute pas supérieure a 20°% en tout 
état de cause. On voit que le désaccord entre les 
conséquences de Vhypothése et les faits observés 
est trop considérable pour ne pas étre inquiétant. 


4, Peut-on suggérer d’autres conditions aux limites? 

On peut remarquer que jusqu ici l’hypothese de 
Vimpédance normale spécifique a fait Vobjet de 
vérifications assez indirectes; c’est ainsi qu’en l’ap- 
pliquant au probléme de la propagation du son 
dans une enceinte rectangulaire, on obtient des 
résultats trés satisfaisants (voir en particulier C. M. 
Harris [6]). Or ce en quoi ce probleme semble 
différer de l’expérience ci-dessus, c’est l’incidence 
des modes propres sur la paroi, car cette incidence 
est quasi nulle dans la salle parallélépipédique, si 
Von admet que l’absorption n’a pas trop troublé les 
modes. L’hypothése discutée serait-elle donc valable 
pour les ondes planes presque normales, et ne le 
serait-elle plus aux fortes incidences? Ceci sug- 
gérerait alors une modification de la condition aux 
limites, faisant intervenir des termes se réduisant 
au gradient normal pour l’incidence nulle: par ex- 
emple, on pourrait ajouter un terme de gradient 
tangentiel, ou remplacer le gradient normal par le 
gradient total. 

Si l’on essaye cette derniére forme dans le cas de 
lexpérience décrite, il faudra remplacer 39/dx par 
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1/cos $39/dx, d’ou 
(l—r)4—(1+r)pc=0 


et m serait alors une constante indépendante de 4. 
Le Tableau II indique les écarts relatifs des valeurs 
précédemment données par rapport 4 leur moyenne 
arithmétique. 

Tableau II 


fréquence [Hz] | 


™ = (my + my + m,)/3 2,19 
My/m 1,07 
m4/m 0,99 
Ma/™m 0,94 


Dans le cas le plus défavorable, l’écart n’est que 
de 10%. 

La suggestion semblerait donc assez séduisante ; 
mais on ne peut songer sérieusement 4 une condi- 
tion telle que 


C) 
“Ia —po- is = 
SPATE) ea? 


(A,: premier paramétre différentiel de Lams), non 
seulement & cause de sa non-linéarité (difficulté qui 
ne serait peut-étre pas dirimante, car on observera 
qu’en vertu de lemme du n° 2 la condition permet 
des solutions séparables par rapport au temps, ce 
qui est l’essentiel), mais encore parce qu’elle con- 
duit & supposer les absorptions sous incidence ra- 
sante aussi fortes que celles sous incidence normale, 
ce que l’expérience contredit. 

Il reste possible que la condition la plus accep- 
table soit celle suggérée en premier lieu, et qui est 
en quelque sorte intermédiaire entre les deux que 
Yon a discutées. Cette condition pourrait s’écrire, 
par exemple, 

Be rade <0 
ou d9/ds désigne la dérivée tangentielle de 9 sur la 
parol: comme 


do \2 do \2 
so= (GY (0) 


on arrive 4 une telle expression en linéarisant la 
précédente, moyennant I’hypothése supplémen- 
taire, que des analogies optiques ne paraissent pas 
rendre choquante a priori, d’un rapport spéci- 
fique de grad, & grad, dans la couche absorbante. 
Le probléme aux limites serait alors ramené a la 
eatégorie des problémes a dérivée oblique donnée, 
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qui ont été étudiés dans la théorie du potentiel, 
notamment par Boutteanp et Grraup [7] et plug 
récemment par Ficnera [8]. 

Si l’on reprend le cas de l’expérience sous inci- 
dence $, on trouvera 4 


An COS 9+ Bs sin 4—e¢ % 
an cos $— as sin (cue ec ; 


dou, 81-3, — Z,,- “exp ia, et 3.= 
1+ Z2 cos? $+ Z? sin? §—2 Z, cos a, sin 4 


™= 7 “cos v,, 608 9—4 ZZ, COB (aq — a) Sim GOSS 


7 7 
ZL, CXp tee, ; 


qui donne naturellement les mémes résultats que_ , 
précédemment pour l’incidence normale, et une~ 
absorption nulle pour l’incidence rasante, comme il 
convient. On pourra encore former, en principe, | : 
une relation entre les m observés sous cing inci- | 
dences distinctes, mais les calculs effectifs sont — 
pénibles. 

Ce ne sont la, bien entendu, que des suggestions, 
qui demanderaient a étre examinées de plus prés, 
tant A la lumiére de la théorie que des expériences. 
Du moins espére-t-on avoir montré que la notion 
d’impédance normale spécifique doit étre utilisée 
avec quelques précautions, et que son domaine de 
validité est encore loin d’avoir été nettement déli- 
mité. 

L’auteur est heureux d’avoir ici occasion de 
remercier I. StRoumza du Laboratoire d’Acous- 
tique duC. R. 8. I. M., quia effectué avec beaucoup 
de soin les mesures dont quelques-unes ont été 
citées dans ce travail. 

(Regu le 15 Aodt 1952.) 
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ON ACOUSTIC DIFFRACTION THROUGH AN APERTURE 
IN A PLANE SCREEN 


by JOHN W. MILES 


University of California, Los Angeles 
Summary 


The scattered wave consequent on the incidence of a monochromatic disturbance on a plane 
screen perforated by an aperture is constructed with the aid of generalized Fourier transforms. 
After formulating the integral equation defining this scattered wave, the total transmission cross 
section is cast in variational form (of the SCHWINGER type) and the result compared with the 
KIRCHHOFF approximation. A similar result is given for scattering by a plane disc. Application 
is made to a circular aperture. 


Sommaire 


On construit, au moyen de la transformation de Fourier généralisée, l’onde diffusée résultant 
de Vincidence d’une perturbation monochromatique sur un écran plan qui est percé d’un orifice. 
Apres avoir établi l’équation intégrale qui définit cette onde diffusée, on calcule la section efficace 
totale de transmission sous la forme d’une variation (du type «SCHWINGER») et on compare le 
résultat obtenu avec l’approximation de KIRCHHOFF. On donne un résultat analogue, relatif a la 
diffusion par un disque plan. Application est faite & un orifice circulaire. 


Zusammenfassung 


Die Beugungswelle, die von einer monochromatischen Stérung an einer Offnung im ebenen 
Schirm hervorgerufen wird, lat sich mit Hilfe der verallgemeinerten Fourierschen Transforma- 
tion darstellen. Nach der Formulierung der Integralgleichung dieser Beugungswelle wird der 
totale Durchlassigkeitsquerschnitt in Variationsform gebracht (nach SCHWINGER), und das Er- 
gebnis mit der KIRCHHOFFSCHEN Naherung verglichen. Ein ahnliches Resultat wird fiir die Beu- 
gung an der ebenen Scheibe angegeben. Als Beispiel wird der Fall einer kreisrunden Offnung im 
ebenen. Schirm behandelt. 


been superfluous after the appearance of the defini- 
tive papers of Levine and Scuwincer [7]. How- 
ever, it seems worthwhile to present the formula- 
tion of the problem via generalized Fourier trans- 
forms, since the evaluation of the variational ex- 
pression for the ‘aperture impedance” often is 
facilitated by the introduction of these transforms. 
(cf. [8]. The application of transform methods to 
the electromagnetic problem has been treated in 
an earlier paper [9], but here also the stronger use 
of the variational principle by LevinE and ScHwINGER 
[10] is to be noted.) Moreover, some additional 
points of interest will be brought in. 


{ 1. Introduction 


We consider here the related questions of diffrac- 
tion through an aperture in a plane screen and 
around a plane disc. The low frequency (wave- 

length large compared to aperture or disc dimen- 
: sions) analysis of this classical problem was given 
originally by Rayierex [1], [2], while at high fre- 
quencies (wave-length small compared to aperture 
or dise dimensions, a situation unlikely to be re- 
alized in audio-acoustics) KircHHOFF'’s approxima- 
tion is valid [2]. Exact solutions are available for 
a half plane [3], a rectilinear slit or strip [4], and 
a circular aperture or disc [5], the complexity of 
which (increasing in the order mentioned) estab- 
lishes the desirability of an approximate method 


2, Basic equations 


to fill the gap between the limiting approximations- 
For this purpose, a suitable vehicle is furnished by 
ScuwinceEr’s variational principle [6], which has 
been applied explicitly to both the circular disc [7], 
and circular aperture [8]. 

The present paper is a condensation of a more 
extensive investigation undertaken in 1947, the 
publication of which (in its entirety) would have 


On the linearized acoustic theory, the vector velo- 
city (q) and perturbation pressure (p) resulting from 
a monochromatic disturbance of pulsatance o may 
be expressed in terms of a dimensionless potential 
(®) according to 


q=c Re[exp(jat)V®], 
P=— ec? Re[jk exp(jot)®] 


(2.1) 
(2.2) 
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where ¢ is the sonic velocity, oe, the mass density of 
the undisturbed medium, k the (dimensionless) wave 
number, viz. (2.3) 


in terms of a characteristic length J, to which all 
subsequent linear measures are referred, and © is 
a solution of the Helmholtz equation 


V0 + k?2O=0. — 


k=al/e 


(2.4) 


In the following, for analytical convenience, we may 
assume k to have a small, negative imaginary part, 
but in the end results it is implicitly real’. 

In general, we pose ® as a function of the position 
vector r, but since we shall focus our attention pri- 
marily on the plane z=0, it is convenient to intro- 
duce the vector e for the value of r there, such that 


r=p+z. (2.5) 


According to the usual custom, z denotes the scalar 
value of z etc. 


3. Boundary conditions 


We consider a disturbance ®,(r) that is incident 
on the plane z=0 from the left (z <0). This plane is 
divided into two (open) regions, o and 7, corre- 
sponding to an aperture and perfectly reflecting 
screen, respectively. The boundary conditions ap- 
propriate to these two regions are continuity of both 
pressure (®) and velocity (V®) across co and van- 
ishing normal velocity (®,) on +. Singularities (in 
V®, not ®) may exist at the boundary (edge) be- 
tween o and +t, but they must be integrable. 

In the absence of the aperture the incident wave 
®,(r) would give rise to the reflected wave ®;(r—2z) 
in z<0, while the field in z>0 would be null. The 
presence of the aperture leads to a transmitted wave 
(r) in z>0, and we therefore are led to write 


®_(r) = 0,(e +2) + O(e—z)— (3.1a) 
®, (r)=0 ( +2) (3.1b) 
where the subscripts —-and + imply z <0 and z> 0, 
respectively. It is evident that O_ and ©; have 
the same normal derivative at z=0, so that it re- 
mains to satisfy the boundary conditions of vani- 


shing normal velocity on the screen and continuity 
of pressure across the aperture, viz. 


®(p—z) > 


0(e)=0, inc; (3.2) 
®(c) =9,(c), ep inc. (3.3) 


In addition, the appropriate (SomMERFELD) con- 
ditions at infinity must be satisfied. For our purpose, 


1 We use the time dependence exp (+jat), whereas LE- 
VINE and SCHWINGER [6], [7] and BOUWKAMP [5] use 
exp (—iwt). For comparison i——j. 


J.W. MILES: ON ACOUSTIC DIFFRACTION THROUGH AN APERTURE 


_ for either diffraction through an aperture or around © 


it suffices to require ® to vanish at infinity when k 
has a negative imaginary part. : 

The problem of diffraction through an aperture o | 
now may be posed as: find a solution to (2.4) satis- — 
fying (2) and (3) and vanishing at infinity. We re- i 
mark that, while this formulation is equally valid » 


a disc, it is appropriate primarily to the former case, — 
a preferable formulation for the latter case being — 
set forth in section 6 below. (But if k is very large © 
the preference tends to be reversed; ef. [6].) 


4, Aperture integral equation 


The classical approach to the problem of deter- 
mining the solution to (2.4) in z>0 for preseribed 7 
values of either ® or ®, on z= 0 is due to RAYLEIGH, | 
who first determined the appropriate Green’s func- | 
tion [11]. Physically interpreted, Raviuten’s result — 
represents the disturbance at any point in space in © 
terms of a source distribution over the plane, cor- © 
responding to the normal velocity there. 

Alternatively, a solution to (2.4) may be exhibited 
in the form 


r= Jk (P98) FO) YO 2): (4.1) 

zt - exp [—j(k2—22) * 2] dS(a) | 

provided that (A, 0) satisfies the two dimensional 
equation (A being a spectral vector) 


vey +a2d=0. (4.2) 


(1) represents the solution as a spectral superposi- — 
tion of the elementary solutions 


H(A, 9) exp [— j(h? — 2") * z] 
with the amplitude distribution 
k(k2 — 22) '* F(x) dS (2). 

On the assumption that the imaginary part of 
k is negative, the branch cuts for the radical] 
are taken from +k to +k + joo in order to satisfy 
the requirement of vanishing © at infinity, or in the 
limit where k is real the 4 path of integration is to 
be indented under —-k and over +k, whence those 
waves for |A| >k are evanescent. We emphasize 
that, in general, > itself may comprise two (scalar) 
parameters (see, e. g., appendix). 

If we introduce the modified aperture velocity, f, 
in accordance with 


®,(r)|, = 
fF) 


Now, in virtue ie the Sturm-Liouville properties 
of (2), the 4(4, 2) can be determined to form a com- 


o=— Af (e) (4.3) 


(1) yields (2: e )dS(a). (4.4) 
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plete (infinite) set and may be normalized in such 
a way that the relation inverse to (4) is 


F(a) ={fos 


Wa.) dSie) (4.5) 


where the bar indicates the complex conjugate. 


(4) and (5) constitute a generalized Fourier trans- 


form pair, examples of which are given in the ap- 
pendix, and F(A) is appropriately designated as the 


“transform” of f(g) 
If we now iavuke the boundary conditions (3.2) 
and (3.3) in (1), we have 


-[k(—2)"* FG) vG, ¢)aSG)— 


Di(oe ), ° ino (4.6a) 
JFQ)¥(, e)dS(2)=0, — pins (4.6b) 
25 

which is a dual integral equation [12] for the deter- 
mination of F(A). 

The integral equation for f(9) is obtained by sub- 
stituting (5) in (6a) and reversing the order of inte- 
gration, whence 


Jes(e, e)f(e')48'= ie), eine; (4.7) 
Bo(°s 6) —/R(B—2) "405 0) 4005 AS (2). (48) 


x 
The original diffraction problem is now reduced to 
the solution of either (6) or (7), after which the total 
fields may be obtained from (3.1), together with (1). 
In the subsequent discussion we shall assume the 


incident wave to be plane and write 
®;(r)= exp(— jkn;:r) (4.9) 
where n,; is a unit vector normal to the wavefront 
(cf. (7.1) infra). 
5. Variational aperture impedance 


The (timewise average) power transmitted through 
the aperture is given by 
=4$ kp c3l*oz 


where o, is the transmission cross section. As a mea- 
sure of the complex power transmitted, it is con- 
venient to introduce the aperture impedance 


(5.1) 


Zg=Ro+jXg=o—/exp(jkni-p) f(e)dS (5.2) 
o 


the real (R,) and imaginary (X,) parts of which re- 
present the propagated and evanescent waves, re- 
spectively, the normalisation being such that 

o,=oR,. (5.3) 
Z, may be cast in variational form? by multiplying 


both sides of (4.7) by f (e) dS, integrating over the 
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aperture and dividing by the square of the magni- 
tude of the resulting, right hand side, whence 


Zs =Yo= 65+ jBo= 
—c/ | fe. fasas’ | 
oo 


(5.4) 
Jexp(jkn, p) fas 


This expression is stationary with respect to first 


order variations of fle) about the true solution to 
(4.7); moreover, it is invariant under a scale trans- 
formation of f. 

In order to exhibit the variational expression in 
terms of the transforms, we substitute the kernel g, 
from (4.8) and, in the denominator, f from (4.4), 
whence (5.5) 

) 


ee FOyeds(a)/|/Fotas(a)) 
ON 
where oi is the transform of ,(¢). 

It is evident from (5) and the prescribed phase 
of the radical that both the real and imaginary parts 
of Y, are positive definite, whence, in virtue of the 
variational principle, the true values are absolute 
minima. We remark that this does not establish the 
real part of Z, (i. e., the transmitted power), albeit 
stationary, aS a minimum. 

It is interesting to observe that, insofar as the 
aperture velocity distribution is assumed 
to be independent of the direction of inci- 
dence (up to a constant multiplier), the transmis- 
sion cross section for oblique incidence is simply 
related to that for normal incidence according to 


= | exp (jkn;-) f(e) as//f(e) ds|*-o,(n,). (5.6) 


This result is valuable primarily for small k sin 9; 
(0; being the polar angle of incidence), since for 
large k sin @; it would be quite wide of the mark to 
assume the aperture field to be independent of the 
direction of incidence. 

The Kircunorr approximation leads to the fol- 
lowing expressions for the aperture impedance 


(ef. [7]) 


Zox =o cost;f fexpl—ikn;-(¢—e)]- (6.1) 
oo 
*8a(P, o’)dS dS’, 
Zox=o- cos; [k (k?—22) (5.8) 
S, 
Ss 2 ate 
+fexp Gino} 40% sb.asf'as 


2 SCHWINGER and LEVINE [7] demonstrate that the 
differential scattering cross-section, as well as the total 
cross-section, is stationary. 
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The result obtained by posing f=f; in the variatio- 
nal expression (4) is designated as the ,,variational 
KircHHorr’’ approximation, and comparison of 
the result so obtained yields the curious relation 


Yon = Zax 800 9; (5.9) 
The comparisons between o,, and o;7° in the 


limits of high and low frequencies have been given 
by Levine and Scuwineur [7]. 


6. Reflection from a disc 


The formulation of the disc problem follows more 
or less by analogy, and we are led to introduce (as 
a measure of the complex, reflected power) the disc 
admittance 


(6.1) 
Y,;=G6,+jBr=7-!cos @; [ exp(—jhn,- ¢) ®(¢) ds 
Az 


where + is the disc area. The reflection cross-section 
then is given by (6.2) 

The corresponding variational expressions are 
Z,=Y, = (6.3) 
= <sec20; | [ Og, ds ds’ / | fexp (jkn; -) dS i 
aU eee Tv (6.4) 
—|2 — — —|2 
o*(a}/'ds(a)/|/Fo*ds(a)| . 

S 
mn 


Gen rae 


Zq= 2 fk (222) 
Px 


In this case, the real and imaginary parts of Z,, the 
reciprocal of Y,, are positive definite and absolute 
minima. 

If it is assumed that O=©; = exp (—jkn,-p) 
thereby violating the boundary condition of conti- 
nuity of pressure at the edge, it is readily shown 
that the resulting integral for Z, obtained from (3) 
is divergent, whence the ordinary and variational 
Kircunorr reflection cross-sections are found to 
be infinite and null, respectively. However, a suit- 
able approximation joining the correct low and 
high frequency limits can be developed by sepa- 
rating out the low frequency limit (k—0) and set- 
ting ®=9,; in the remainder, whence 


kZ.(k)= + sec?@; [ pg )dS dS’ /| fo, as]? + 
TT Tv 


+71 sect@, /[(k2@—22)"* —jjay]- (6.5) 
Sy 


—> —> 


-| fexp(ikn,-e) ¥(3, 9) aS|" a8 (3), 


Bole’) = limkg,=j /irpe eu Pras) (66) 
> rt 


where ®, (p) is the static pressure distribution. 
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In actually evaluating the second term in (5) it © 
is frequently convenient to note that it is directly 
related to the integral of (5.8) with respect to k, | 
provided that k is held constant where it occurs in | 
the incident wave. / 

Eq. (5) is not,-of course, a consistent approxi- 
mation, since it poses a different © at k = 0 than 
for finite k, and it is, therefore, no longer a varia- 
tional approximation in the sense that the error 
in Z, is of the order of the square of the error in the 
assumed field. Nevertheless, it has the merits of 
giving the correct limits at both k = 0 and k = ~ 
and furnishing a continuous and rationally derived 
approximation at intermediate points. (Cf. the 
treatment of electromagnetic diffraction by a cir- 
cular disc or aperture in [9].) 

In concluding this section, we remark that Ln- 
vinE and ScHWINGER [7] start with the dise for- 
mulation for an aperture to obtain a second varia- 
tional formulation (finite c) that is superior to that 
of section 5 at high frequencies. Similarly, the ana- 
logue of this second principle for reflection from a 
disc could be obtained by starting with the formu- 
lation of section 3. 


7. Circular aperture 


The diffraction of a normally incident plane wave 
through a circular aperture has been treated pre- 
viously by variational methods in [8], while the case 
of a circular disc has been treated quite thoroughly 
by Scuwinesr and Levine ([7], after appropriate 
changes in terminology, since they treat the equi- 
valent problem of diffraction through an aperture 
subject to the boundary condition ® = 0 on the 
screen). We shall consider here only the case of 
oblique incidence on a circular aperture, approxi- 
mating the aperture admittance on the basis of the 
static aperture field. 

In polar coordinates (p, e), a plane wave incident 
at polar angle ©; and azimuthal angle 9; is speci- 
fied by (7.1) 

®;= exp { —jk [p sin 9; cos (p—9;) + z cos @;]}. 


The (3, °) are given in the appendix; assuming f (ec) 
to be independent of ¢ only the term n = 0 enters, 
and substitution in (5.4) yields 


oo 1 2x 
Y,=tk il (k2—22) J [ Io(Ae) f () ede dol 


< ae (7.2) 
‘| / [exp [—ike sin ©; cos (9—e,)]f (e) edede|*ada 
00 


where the characteristic length is taken as the ra- 


dius of the aperture. Carrying out the angular inte- 
grations reduces this to 
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Yo=tk/(e—xy 
r (7.3) 


1 1 
-| [ Joe) F (0) ¢ de | f Jo(ke sin ©) f (e) e del?ada. 
0 0 


The static aperture velocity is specified by 
fley=(—e) 
Substituting this in (3), the radial integrations can 


be effected via W 12.11(1) and 3.4(3) (W refers to 
Warson’s treatise [13]) with the result 


Y,o=$k [k sin 9;/sin (k sin 9;)]?- 


(7.4) 


(7.5) 


fx (2 92)" sin? a da. 

ny 
To carry out the A integration we introduce the 
change of variable »=k, differentiate with respect 
to k and identify the real and imaginary parts of the 
result via W 10.4(1) and 6.14(3), respectively, whence 


¥e= 3 [k sin ©;/sin (k sin @,)}?- 
2h 


[ [Hy (*) + jJo(x)] dx 


(7.6) 


where H, is Struve’s function. 

The result (6) may be expected to furnish a good 
approximation for small k. (At normal incidence [8] 
it is quite satisfactory up to k=1.) Expanding in 
powers of k, and substituting in (5.3), we find [14] 


= = S [1+ (0.03916 —}sin?@,) k2 + O(k4)] (7.7) 
which is exact to the approximation indicated, in 
virtue of the variational principle together with the 
use of the aperture field that is exact at k=0. (In 
[8] the number 0.03916 was given as 0.03917 due 
to a small numerical error, the first number being 
correct, as determined by Bouwxkamp’s analysis. 
The fact that the variational approximation must 
be exact to this order was pointed out by H. Lr- 
VINE [15].) 
8. Appendix 


In Cartesian coordinates (x, y), (4.2) reads 
Vest wy tNY=0, p=ixtjy (Al) 
and the appropriately normalized solutions are 


4 (A, e)= (2n)— exp (j°9), (A 2) 
A=iptjy, aS(A)=dyudy (A 3) 
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and (4.4) and (4.5) reduce to the well known, double 
Fourier transform pair. 
In polar coordinates (pe, ©), (4.2) reads 


2 (Pb .)o+ Yoo + (Ae)? 4= 0 (A 4) 


and the required solutions and corresponding trans- 
forms are 


> => 


(2, e)= (2m) "J, (Ae) exp(jng), (A 5) 
flee=)>, [FameGejrdr,, (a6) 
n=—oQ0 
oo 27¢ Ses 
F(a,n)= f f f(e,e) ¥ (a, e)e dede. (A 7) 


In this case the integration over the % spectrum 
must be interpreted as including the summation 
over the discrete spectrum (n) associated with the 
angular (9) distribution of f(e, e), in addition to 
integration over the spectrum associated with 
the radial variable (p). The transform pair of (A 6) 
and (A 7) is, of course, the well known Fourier- 


Bessel representation. 
(Received 26 th February, 1952.) 
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LETTER TO THE EDITOR 


LE DALE Re bOs THE. EDITOR. 


GAP-LENGTH FORMULA IN MAGNETIC RECORDING 


by W. K. WESTMIJZE 
Philips Research Laboratories, N. V. Philips’ Gloeilampenfabrieken, Eindhoven, Netherlands 


In the optical reproduction of a sinusoidal record- 
ing (density as well as amplitude modulation) the 
loss produced by the finite width d of a homogene- 
ously illuminated slit on the reproduction of a re- 
corded wavelength % is given by the well-known 
formula 


which is found by a simple integration. For large 
wavelengths the output is proportional to d. 

This formula is sometimes transferred to the case 
of the loss produced by the finite length / of the gap 
of a magnetic reproducing head. The processes how- 
ever being quite different, this is certainly not 
correct. For this case LtUpuckx [1] derived under 
simplifying conditions the analogous formula 


sin (rl/A) 
G (rl/a) = eran (1) 
As is to be expected on physical grounds G is inde- 

pendent of | for large wavelengths. 
This however is to be con- 
b sidered only as a first ap- 
proximation. A more de- 
tailed analysis shows that 
for a tape of unit relative 
permeability and a repro- 
duction head with pole 
pieces of infinite permeabi- 
lity and of the shape illus- 
trated in Fig. 1, the gap 
loss is given by the formula 


Fig. 1. 
Cross-section of pole tips 
(a) and tape (b) of the 
type of head discussed. 


1/2 


S(xl/) = ah x sin . (ton x—x + 3) dx + 
T r Y. 
0 


1 
Leeeerch x 21 1 1—x 
+ A sin “fe o|5 (2+ 2h 155) 


0 


S represents a damped oscillating function of 
analogous behaviour as G. The most important 
differences are 


1) the first zero occurs for 10.88 and not as in 
(BE) forta 2c 

2) for large values of the argument, S may be 
approximated by 


sin (wl/% + 7/6) 


S (1/2) =0.699 Gal) 


(2) 


and therefore the amplitude of the successive 
extremes decreases less rapidly than for G. 


Hence the gap length as calculated from the wave- 
length of the first zero with 1= \ gives a value for the 
“effective gap length” that is 10 % larger than the 
actual gap. 

It has been proposed to calibrate the recorded 
flux as a function of wavelength with the aid of wide- 
gap measurements [2]. Experimentally it was found 
that in order to obtain results consistent with 
narrow-gap measurements a decrease of 2 dB per 
octave had to be applied ifformula (1) was used. This 
is in accordance with our calculation that predicts a 
rise of S with respect to G proportional to (rl/A)'>. 
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Fig. 2. Voltage E across reproduction coil (arbitrary units) 
against f-l/v; f—frequency, gap length, v=tape 
speed. ; 


In Fig. 2 the reproduced signal for a tape with 
constant amplitude of the magnetisation as caleu- 
lated with eq. (2) is plotted against the dimension- 
less quantity f -l/v(—1/,), where f is the frequency 
and v the tape speed. 

A more detailed review will be published shortly 
in the Philips Research Reports. 
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